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Abstract 


This  experiment  studied  the  effect  of  streamwise  vortices  on  unsteady  turbulent 
boundary-layer  separation.  The  objectives  were  to  document  the  flow  field,  to 
characterize  the  time  response  of  the  boundary  layer,  and  to  understand  the  actual 
mechanisms  by  which  the  streamwise  vortices  modify  boundary-layer  behavior. 

A  new  configuration  for  non-obtrusive  three-component  Laser  Doppler  Anemom- 
etry  (LDA)  determined  the  phase  averaged  velocity  and  Reynolds  stress  compo¬ 
nents,  in  an  unsteady  water  tunnel,  at  a  momentum  thickness  Reynolds  number  of 
1840.  The  streamwise  vortices  were  created  by  three  pairs  of  half-delta  wing  vortex 
generators,  while  the  boundary-layer  separation  was  controlled  through  impulsively 
initiated  opposite-wall  suction,  which  created  a  strong  adverse  pressure  gradient. 

The  time  response  of  the  freestream  velocity  demonstrates  that  convection  is 
the  primary  mechanism  by  which  vortex  generators  modify  the  response  of  the 
boundary  layer.  Modeling  the  response  with  two  time  scales  proved  to  be  a  useful 
technique.  There  is  an  initial  fast  response  throughout  the  boundary  layer  which 
is  unaffected  by  the  presence  of  vortex  generators,  followed  by  a  slow  or  convective 
response,  the  magnitude  of  which  is  substantially  modified  by  the  presence  of  the 
vortex  generators. 

In  the  presence  of  an  adverse  pressure  gradient,  the  strength  of  the  embedded 
vortex  decreased  subtantially  after  several  convective  time  periods.  In  addition,  the 
spanwise  gradient  of  (u)  appeared  to  produce  streamwise  vorticity  of  the  opposite 
sign  by  straining  the  spanwise  vorticity  of  the  boundary  layer.  All  the  signifi¬ 
cant  transport  of  spanwise  vorticity  appeared  to  occur  in  the  viscous  region  of  the 
boundary  layer.  While  the  regions  of  turbulent  kinetic  energy  production  were  mod¬ 
ified  by  the  presence  of  streamwise  vorticity,  this  did  not  appear  to  be  a  dominant 
mechanism  in  the  interaction  between  the  streamwise  vorticity  and  the  separating 
boundary  layer.  These  new  insights  may  be  helpful  to  persons  attempting  to  model 
this  very  complex  unsteady  turbulent  flow. 
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1.  Introduction 


1.1  Background 

The  purpose  of  this  program  was  to  understand  the  effect  of  vortex  generators 
on  an  unsteady,  separating,  turbulent  boundary  layer.  Vortex  generators  have  been 
in  use  on  aircraft  since  the  mid-1940’s  and  it  is  well  established  that  they  help  delay 
flow  separation.  Recently,  a  number  of  researchers  have  studied  vortex  generators 
in  steady  flow  to  understand  the  nature  of  vorticiiy  transport  and  the  changes  in  the 
turbulence  structure  resulting  from  streamwise  vorticity  embedded  in  the  boundary 
layer.  While  steady  flow  is  the  logical  first  step  in  understanding  the  effects  of  vortex 
generators,  much  of  the  interest  in  these  flow  control  devices  relates  to  applications 
in  unsteady  flow  regimes,  such  as  supermaneuverable  aircraft.  The  next  logical 
step,  then,  is  to  understand  the  differences  between  how  vortex  generators  behave 
in  steady  flow  and  how  they  behave  in  unsteady  flow.  This  program  was  intended 
to  be  a  first  step  in  that  direction. 

While  vortex  generators  are  presently  used  primarily  as  passive  control  devices, 
future  applications  may  require  their  use  as  active  control  devices.  During  the 
design  of  an  active  control  system,  engineers  need  to  know  the  step  response  of  the 
various  components  of  the  system.  In  this  program,  the  portion  of  the  system  that 
is  to  be  controlled  is  the  boundary-layer  on  a  flat  surface,  and  the  disturbance  to  the 
system  is  a  step  change  from  zero  to  adverse  pressure  gradient.  The  response  to  this 
disturbance  must  be  known  for  the  system  both  with  and  without  vortex  generators. 
The  results  from  this  program  provide  this  basic  step  response  information  as  well 
as  some  insight  into  the  mechanisms  by  which  vortex  generators  modify  the  nature 
of  an  unsteady,  separating,  turbulent  boundary  layer. 

All  the  data  for  this  program  were  obtained  in  a  computer-controlled,  unsteady 
water  tunnel,  using  a  compact  three-component  Laser  Doppler  Anemometry  (LDA) 
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system  developed  to  meet  the  constraints  of  this  facility.  An  array  of  six  vor¬ 
tex  generators  generated  streamwise  vortices  in  a  steady  region  of  the  flow.  As 
these  vortices  convected  downstream,  an  adverse  pressure  gradient  was  impulsively 
started.  Phase-averages  of  the  three-component  LDA  measurements  documented 
the  development  of  the  mean  velocities  and  all  six  of  the  normal  and  shear  Reynolds 
stresses. 

1.2  Program  Objectives 

The  specific  objectives  of  this  program  were  as  follows: 

(1)  Describe  how  vortex  generators  change  the  nature  of  the  step  response  of  the 
boundary  layer  to  an  unsteady  adverse  pressure  gradient. 

(2)  Understand  the  mechanisms  by  which  streamwise  vortices  modify  the  behav¬ 
ior  of  a  boundary  layer  subjected  to  an  unsteady  adverse  pressure  gradient. 

(3)  Develop  and  test  a  compact  three-component  LDA  system  that  is  compatible 
with  the  constraints  of  the  unsteady  water  tunnel  used  in  this  program. 

(4)  Provide  a  set  of  profiles  of  velocity  and  Reynolds  stresses  taken  under  un¬ 
steady  flow  conditions  to  supplement  the  results  already  published  for  steady  flow. 

Chapters  2  and  3  will  present  the  governing  equations  and  previous  work  pub¬ 
lished  in  the  literature.  Chapter  4  will  describe  the  experimental  facility  and  the 
LDA  configuration,  while  chapter  5  will  discuss  the  validation  of  the  facility  and 
LDA  systems.  Chapter  6  will  outline  the  experimental  plan,  while  chapter  7  will 
present  the  measured  velocity  profiles  and  Reynolds  stresses.  Finally,  chapter  8  will 
analyze  the  results  in  terms  of  time  response  and  mechanisms.  This  report  will  then 
conclude  with  a  concise  summary  in  chapter  9. 
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2.  The  Governing  Equations 


The  governing  equations  can  be  expressed  in  many  different  forms  using  a  variety 
of  notations  and  symbology,  and  in  fact,  they  should  be  expressed  in  the  form  that 
is  most  useful  to  the  particular  study  at  hand.  This  chapter  provides  a  brief  review 
of  the  relevant  governing  equations,  stated  in  terms  that  are  most  relevant  to  this 
program. 

2.1  Alternate  formulations  of  the  governing  equations 

For  steady  turbulent  flows,  the  standard  practice  is  to  identify  each  variable  as 
having  a  time-averaged  mean  component  and  an  instantaneous  fluctuating  compo¬ 
nent.  That  is, 

u  =  U  +  u'y  (2.1) 

where  u  is  the  instantaneous  velocity  composed  of  U,  the  time-averaged  component, 
and  u',  the  fluctuating  component. 

In  unsteady  flows,  the  time-averaged  value  may  or  may  not  have  significance, 
depending  on  the  frequency  of  the  unsteadiness  and  the  response  time  of  the  flow. 
What  is  generally  of  more  interest  is  the  “average”  value  that  would  be  seen  at  a 
particular  point  or  phase  of  the  unsteady  event.  If  the  same  event  is  repeated  many 
times,  then  a  phase-averaged  value  can  be  defined  at  each  phase.  That  is, 

u  =  (u)  +  u\  (2.2) 

where  u  is  the  instantaneous  velocity  composed  of  (u),  the  phase-averaged  velocity 
and  u',  the  fluctuating  component.  If  the  event  being  described  consists  of  an 
initial  condition  followed  by  a  disturbance  and  a  response  time,  then  this  double 
decomposition  is  probably  the  appropriate  description. 
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If,  on  the  other  hand,  the  event  being  described  is  actually  cyclic  in  nature,  then 
a  triple  decomposition  is  appropriate. 

u  —  U  +  u  +  u'  (2.3) 

In  equation  2.3,  u  is  the  deterministic  component  of  the  instantaneous  velocity. 
There  is  clearly  a  connection  between  this  component  and  the  phase-averaged  ve¬ 
locity,  which  is 

(u)  =  U  +  u.  (2.4) 

In  practice,  u  is  found  by  first  obtaining  (u)  and  then  subtracting  the  mean  velocity, 

U. 

In  this  experiment,  the  pressure  gradient  was  suddenly  changed  and  the  flow 
was  allowed  to  fully  adjust  to  that  change.  Under  these  conditions,  neither  the 
time-averaged  values  nor  the  deterministic  values  have  much  meaning  by  them¬ 
selves.  Therefore,  the  quantities  measured  and  calculated  will  be  expressed  as 
phase-averaged  quantities.  For  comparitive  purposes  and  to  establish  a  baseline 
for  symbology,  the  phase-averaged  equations  are  presented  in  section  2.2,  and  the 
deterministic  equations  are  summarized  in  section  2.3. 

Cartesian  Tensor  Notation.  Since  many  of  the  following  equations  are  lengthy 
in  their  most  compact  form,  cartesian  tensor  notation  will  frequently  be  used 
throughout  this  report.  The  conventional  orthogonal  axes  x,  y,  and  z  will  be  in¬ 
dicated  as  *i,  X2,  and  13,  and  an  arbitrary  direction  will  be  indicated  as  ir  In  a 
similar  fashion,  the  conventional  velocities  u,  v,  and  w,  will  be  indicated  as  uj,  u 2, 
and  U3.  An  overdot,  u,  denotes  the  partial  derivative  with  respect  to  time,  and  a 
comma,  ut,y,  denotes  a  partial  derivative  with  respect  to  the  direction  Xj.  Also,  in 
this  notation,  any  indices  that  are  repeated  in  a  given  term,  such  as  ut,*,  indicates 
a  summation  of  that  term  over  all  three  directions.  When  a  summation  is  not  in¬ 
tended,  Greek  indices  will  be  used,  such  as  ua,a.  This  symbology  will  become  clear 
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in  the  context  of  the  equations  that  follow.  Finally,  the  isotropic  tensors  and 
ty jk  will  occasionally  be  used.  They  are  defined  as  follows. 


*ii  = 

/  1  if  »  =  J 

\  0  otherwise 

(2.5) 

^={o' 

if  ijk  belongs  to  123123 
if  ijk  belongs  to  321321 
otherwise 

(2-6) 

2.2  The  Phase- Averaged  Governing  Equations 

The  phase-averaged  equations  take  on  the  same  form  as  the  conventional  time- 
averaged  turbulence  equations,  except  that  the  time-average  is  replaced  in  all  places 
by  the  phase-average.  Since  the  development  of  these  equations  is  found  in  many 
textbooks,  they  are  presented  only  in  summary  fashion  below.  In  all  cases,  the 
fluid  is  assumed  to  be  incompressible ,  a  very  good  assumption  for  this  experiment 
in  water. 

Continuity 


(uy)»y  ~  0)  —  ® 


(2.7) 


Momentum 


(«*•>  +<UyUy),y  =  ^  (p)  ,t  +u(ux),JJ  -(u'u'),y 


Evolution  of  Turbulent  Kinetic  Energy 


=  P  —  D  — 


(2.8) 


(2.9) 
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In  equation  2.7,  P  is  the  production,  D  is  the  dissipation,  and  Jy ,y  is  the  flux  of 
turbulent  kinetic  energy  in  the  jth  direction. 


Evolution  of  the  Reynolds  Stress 

+  (<«*> (utu;))  1*  =  pij  +  Tij  -  Ay  -  Ji]k,k  •  (2.10) 

In  equation  (2.8),  in  addition  to  the  production,  dissipation,  and  flux  named  previ¬ 
ously,  Tij  is  the  pressure  strain. 

Pij  =  -(*'*«;•)  («»)>*:  -(ujfctt  •)(“>>»* 

t*'j  =  “(p'(uiv+ttyn)> 

Ay  =  2u{u'i,ku,j,k) 

Jijk  =  “((Ay)*.-*  +  (AJ)fyk)  +  (ukuiu'j)  ~  vWiu']),k  • 


Evolution  of  Vorticity 

In  this  study,  evolution  of  both  the  spanwise  vorticity  and  the  streamwise  vor¬ 
ticity  are  of  interest.  The  spanwise  vorticity  equation  characterizes  the  vorticity 
embedded  in  the  boundary  layer,  while  the  streamwise  vorticity  equation  charac¬ 
terizes  the  vorticity  introduced  into  the  flow  by  the  vortex  generators.  The  following 
equation  applies  to  both  streamwise  and  spanwise  vorticity. 


6 


K)  +((«jK«i»,>  =  (w>X«i)y  +‘,(ui)<ii  +(u>“J>y )  -  (“y“J)>j  • 


(2.11) 


An  alternate  formulation  does  not  transform  the  turbulent  velocity  products  into 
turbulent  vorticity  terms.  It  is  useful  because  it  expresses  the  equation  in  terms  of 
Reynolds  stress  terms,  which  can  be  directly  measured  in  an  experiment. 

(«*)  +«**y)<"»))»y  =  (wy)(«»).y  +*'("*).«  -eijk(u'qu'k)’qj  ■  (2-12) 

Equation  2.11  is  particularly  helpful  in  that  it  lends  insight  into  the  roles  the 
Reynolds  stress  terms  play  in  equation  2.12.  The  Reynolds  stress  terms  are  a 
combination  of  turbulent  convection  and  turbulent  vortex  stretching,  and  will  in 
general  be  referred  to  as  turbulent  transport.  The  vorticity  equation  will  be  ex¬ 
panded  term  by  term  and  discussed  at  length  for  both  streamwise  and  spanwise 
vorticity  in  chapter  8. 

2.3  The  Deterministic  Governing  Equations 

The  key  deterministic  equations  governing  incompressible,  constant  property  flow 
are  presented  below.  Following  the  philosophy  expressed  in  section  2.2,  since  these 
equations  have  been  developed  in  various  forms  by  a  number  of  sources,  including 
Brereton  and  Reynolds  (1987),  the  resulting  equations  will  not  be  re-derived  but 
are  presented  in  a  summary  fashion  below. 

Continuity 


Uj  ij  —  0,  Uj  ,j  0, 


“>•>  =  0 


(2.13) 


Momentum 


ut  +  (Uj&i  +  UjUi  +  u'u'  +  UiUj),j  =  —  ~pn  +!/«*•  ,yj 


(2.14) 
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where 


pij  =  ~  UJ buifi*»fc  ~  uluSfirk  -u'ku'jUi’k  -u'ku'iUvk 


Tij  =  ~p  P'K.J  +“y«  ) 


Dij  =  2l/u'i,ku'j,k 


Jijk  =  -(p'u'Ajfc  +  P>ui^jk)  +  “fc“iuy  +  «jb“i«y  — 


Evolution  of  Vorticity 

«.*  +  =  (n;ut,y  +wjfcf/t-,y  ) 

-(uynt-  +  uywt  +  u'u/|),y  +(wyut  +  wju{),y  H-^u>t,yy  . 


(2.18) 
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Mathematical  theories  from  the  happy  hunting 
grounds  of  pure  mathematics  are  found  suitable 
to  describe  the  airflow  produced  by  aircraft  ... 
—  Theodore  von  Karman,  1954 
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3.  Literature  Review 


Since  this  program  brings  together  two  technologies,  the  study  of  streamwise 
vortices  in  steady  flows  and  the  study  of  unsteady  separating  turbulent  flows,  the 
first  two  sections  of  this  review  will  be  divided  in  the  same  manner.  In  the  third 
section  of  this  review,  selected  papers  will  be  discussed  which  consider  the  combined 
subjects  of  vorticity  in  unsteady  flows.  Because  there  is  a  great  deal  of  literature 
in  each  of  the  first  two  fields,  only  the  articles  which  relate  most  directly  to  this 
program  will  be  discussed  here.  For  a  broader  review  of  each  separate  subject 
area,  Pauley  and  Eaton  (1988)  reviewed  streamwise  vortices  and  Reynolds  and 
Carr  (1985)  reviewed  unsteady,  driven,  separated  flows. 

3.1  Streamwise  Vorticity  in  Steady  Flows 

Vortex  generators  have  been  in  use  on  aircraft  since  the  mid-1940’s.  Schubauer 
and  Spangenberg  (1960)  attribute  the  first  vortex  generator  to  H.  D.  Taylor  in  1946. 
Much  of  the  early  work  was  geared  toward  specific  design  applications  rather  than 
fundamental  studies  of  boundary  layer  interaction.  In  addition,  since  many  of  the 
results  were  of  a  proprietary  nature,  they  were  only  published  in  corporate  reports 
or  memos.  This  section  will  begin  with  one  of  the  earliest  works  published  in  an 
engineering  journal. 

Schubauer  and  Spangenberg  (1960)  studied  forced  mixing  to  “increase  the  pres¬ 
sure  recovery  that  a  boundary  layer  will  withstand  before  separation  occurs.”  They 
evaluated  various  mixing  devices  in  air,  including  plows,  scoops,  twist-interchangers, 
ramps,  tapered  fins,  domes,  shielded  sinks,  and  vortex  generators.  The  vortex  gen¬ 
erators  were  of  a  trapezoidal  shape  recommended  by  Taylor  in  1948,  and  were 
arranged  in  pairs,  producing  counter-rotating  vortices.  They  learned  that  “mixing 
on  a  coarse  scale  by  relatively  large,  widely  spaced  devices”  was  far  more  effective 
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them  fine  scale  mixing,  and  that  under  these  conditions  multiple  rows  were  less  ef¬ 
fective  than  a  “single  row  of  devices  properly  spaced  and  properly  stationed.”  They 
concluded  that  while  the  mixing  devices  differ  considerably  in  their  effectiveness 
and  drag  penalty,  they  all  affect  the  mechanics  of  flow  in  about  the  same  way. 

Pearcey  (1961)  studied  the  tradeoffs  between  various  vortex  generator  configura¬ 
tions.  He  considered  co-rotating  pairs,  counter-rotating  pairs,  and  alternate  pairs  of 
co-rotating  and  counter-rotating  which  give  the  so-called  biplane  effect.  He  noted 
that  the  key  design  issues  are  the  strength  and  the  disposition  of  the  individual 
vortices  at  the  location  where  separation  is  to  occur.  The  effectiveness  of  each  vor¬ 
tex  generator  was  governed  by  the  initial  spacing,  height,  strength  and  direction 
of  rotation.  In  addition,  the  local  effects  of  the  generator  blades  on  the  boundary 
layer  itself  were  important.  Co-rotating  vortices  were  effective  when  their  initial 
spacing  was  greater  than  three  times  the  height.  Counter-rotating  systems  of  vor¬ 
tices  had  the  characteristic  that  they  moved  away  from  the  surface,  which  limited 
their  effectiveness  to  generally  20  generator  heights.  Pearcey  recommended  dimen¬ 
sional  ratios  for  optimum  performance,  allowing  for  trade-offs  between  maximum 
effectiveness,  range  of  effectiveness,  and  minimum  drag.  As  a  baseline,  he  recom¬ 
mended  co-rotating  pairs  because  of  their  greater  range  of  effectiveness  and  their 
insensitivity  in  performance  to  blade  size,  incidence,  and  spacing. 

Perkins  (1970)  evaluated  the  formation  of  streamwise  vorticity  in  terms  of  the 
steady,  incompressible,  turbulent  vorticity  equation.  He  discussed  secondary  cur¬ 
rents  of  the  first  kind  and  second  kind  as  originally  categorized  by  Prandtl  in  1952. 
Secondary  flows  of  the  first  kind  are  caused  by  vortex  stretching  due  to  skewing  the 
mean  flow,  whereas  secondary  currents  of  the  second  kind  are  caused  by  turbulent 
transport  and  turbulent  vortex  stretching.  Perkins  systematically  evaluated  the  rel¬ 
ative  importance  of  each  of  the  terms  in  this  vorticity  equation.  While  his  treatise 
emphasized  the  natural  evolution  of  streamwise  vorticity,  his  conclusions  also  apply 
to  the  growth  or  decay  of  vortices  that  are  intentionally  introduced  into  a  turbulent 
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boundary  layer.  The  two  terms  that  he  found  were  of  major  import ance  were  the 
turbulent  transport  terms,  -^fcW2  —  w'2)  and  (jp  —  jpjv'w1. 

Mehta,  Shabaka,  Shibl,  and  Bradshaw  (1983)  were  one  of  the  first  to  report  re¬ 
sults  for  longitudinal  vortices  embedded  in  turbulent  boundary  layers.  They  looked 
at  an  isolated  vortex,  a  vortex  pair  with  “common  flow  down”,  and  a  vortex  pair 
with  “common  flow  up”.  They  introduced  the  vorticity  in  the  settling  chamber 
of  a  wind  tunnel,  then  measured  mean  velocities  and  the  Reynolds  stresses  after 
the  flow  went  through  a  9  to  1  two-dimensional  contraction.  They  concluded  that 
the  longitudinal  vorticity  drastically  altered  the  turbulence  structure.  The  velocity 
contours  were  qualitatively  as  would  be  expected  due  to  effects  of  passive  convec¬ 
tion  by  the  secondary  flow.  Triple  products  suggested  that  turbulent  transport  of 
Reynold’s  stresses  occurred  by  gradient  diffusion,  quantification  of  which  they  felt 
would  be  very  difficult. 

Shabaka,  Mehta  and  Bradshaw  (1985)  reported  in  greater  detail  on  their  results 
from  a  single  vortex  embedded  in  a  turbulent  boundary  layer.  They  presented 
contours  of  the  longitudinal  vorticity,  the  Reynolds  stress  terms,  the  eddy  viscosity, 
and  the  turbulent  transport  of  turbulent  kinetic  energy  and  shear  stress.  They 
observed  a  region  of  negative  shear  stress  on  the  downwash  side  of  the  vortex  and 
regions  of  vorticity  of  opposite  sign  surrounding  the  main  vortex.  There  were  large 
variations  in  the  shear  correlation  coefficient  and  the  stress-energy  ratio,  and  the 
correlation  coefficients  for  v'w'  and  u'w'  were  close  in  value  to  u'v1.  Finally,  they 
concluded  that  because  contours  of  the  eddy  viscosity  were  far  more  complex  than 
related  contours  of  the  turbulent  transport  inferred  from  the  triple  products,  that 
the  concept  of  a  transport  velocity  is  more  plausible  than  the  concept  of  an  eddy 
viscosity.  This  implies  that  stress  transport  models  may  be  needed  to  successfully 
compute  such  flows. 

Mehta  (1984)  studied  the  interaction  of  a  vortex  with  a  mixing  layer.  He  observed 
that  the  vortex  induced  strong  secondary  motions  and  to  some  degree  redistributed 
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the  turbulence,  in  particular  u'w1.  In  general,  there  was  nc  active  interaction  be¬ 
tween  the  vortex  and  the  mixing  layer. 

Mehta  (1985)  reported  the  effect  of  a  longitudinal  vortex  on  a  separated  boundary 
layer.  He  made  two-component  LDA  measurements  on  an  axisymmetric  model  at 
subsonic  velocities  in  a  transonic  wind  tunnel  and  generally  saw  consistent  trends 
with  those  reported  for  the  mixing  layer.  The  vortex  delayed  separation  on  the 
downwash  side  of  the  vortex  but  caused  earlier  separation  on  the  upwash  side.  Of 
particular  interest  was  the  large  diffusion  of  u'v1  away  from  the  surface  as  indi¬ 
cated  by  the  triple  correlation  u'v ,2,  and  the  large  production  of  u'v',  indicated  by 
v,2(dU/dy).  He  suggested  that  the  vortex  first  generates  the  normal  stresses,  t/2 
and  u/72,  which  then  produce  u'v'  by  interaction  with  the  mean  strain  rates.  He 
also  observed  a  constant  value  of  the  shear  stress  correlation  coefficient  well  into 
the  boundary  layer. 

Mehta  (1986)  reported  on  the  same  experiment  performed  at  transonic  speeds. 
In  contrast  to  his  subsonic  results,  the  vortex  reduced  the  extent  of  the  separated 
region  and  caused  reattachment  to  occur  sooner.  Also  in  contrast  to  the  subsonic 
results,  the  vortex  core  had  much  lower  values  of  v72  and  there  was  no  increase  in 
u'v12  on  the  upwash  side  of  the  vortex. 

Westphal,  Eaton,  and  Pauley  (1985)  investigated  the  development  of  a  fairly  weak 
streamwise  vortex,  embedded  within  a  turbulent  boundary  layer,  in  the  presence  of 
a  moderate  pressure  gradient.  Their  study  was  done  in  air.  They  concluded  that 
the  vortex  could  be  characterized  quite  well  using  the  contours  of  the  streamwise 
vorticity.  Second,  they  observed  substantial  core  growth  relative  to  that  expected 
for  a  free  vortex.  While  the  accelerated  core  growth  was  anticipated  due  to  the 
adverse  pressure  gradient,  the  growth  was  more  than  could  be  accounted  for  by 
Batchelor’s  (1967)  simple  theory.  As  the  core  grew,  due  to  interaction  with  the 
wall,  the  core  took  on  an  increasingly  elliptical  shape.  In  contrast  to  claims  by 
other  authors,  the  vortex  did  not  meander  as  it  traveled  downstream.  Although 
image  effects  caused  the  vortex  to  stray  from  the  center  line  of  the  plate,  it  did  so 
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hi  o.  consistent  «viid  predictable  manner.  Distortion  of  the  Reynolds  stresses  were 
also  more  pronounced  in  the  presence  of  the  adverse  pressure  gradient.  Finally, 
they  observed  an  apparent  thickening  and  thinning  of  the  boundary  layer  on  the 
respective  upwash  and  downwash  sides  of  the  vortex. 

Eibeck  and  Eaton  (1985)  looked  at  the  effects  of  embedding  a  single  longitudinal 
vortex  in  an  otherwise  two-dimensional  turbulent  boundary  layer.  They  measured 
three  components  of  the  mean  velocity  using  a  4-hole  pressure  probe,  a  3-hole  probe 
and  a  single  hot  wire.  They  measured  shear  stress  using  a  2-component  surface  fence 
gauge.  In  the  downwash  region,  the  vortex  created  a  locally  thinner  boundary  layer 
with  subsequently  higher  skin  friction  and  heat  transfer.  In  the  upwash  region, 
it  created  a  locally  thicker  boundary  layer  with  subsequently  lower  skin  friction 
and  lower  heat  transfer.  While  there  were  3-dimensional  distortions  in  the  outer 
flow  (y+  values  greater  than  200),  the  flow  near  the  wall  (y+  values  less  than  200) 
remained  essentially  2-dimensional  as  was  evidenced  by  the  log-linear  profiles  which 
were  similar  for  the  entire  flow  field. 

Liandrat,  Aupoix,  and  Cousteix  (1985)  computed  the  effect  of  longitudinal  vor¬ 
tices  embedded  in  a  turbulent  boundary  layer.  They  found  that  simple  models 
based  on  the  Boussinesq  hypothesis  worked  well  for  a  single  vortex  but  failed  in 
part  for  the  more  complex  case  of  a  vortex  pair.  By  integrating  the  Reynolds  stress 
transport  equations  along  streamlines,  and  with  the  help  of  a  simple  diffusion  model, 
they  obtained  better  results.  Their  results  emphasized  the  need  for  good  data  in 
building  models  for  the  transport  of  the  Reynolds  stress. 

Cutler  and  Bradshaw  (1986)  studied  the  interaction  between  a  strong  longitudinal 
vortex  and  a  turbulent  boundary  layer  in  air.  A  full  delta  wing  placed  in  front  of 
a  flat  plate  allowed  the  trailing  vortices  to  pass  onto  the  plate  while  the  wake 
passed  under  the  plate.  They  measured  the  mean  velocities  and  pressure,  all  of  the 
Reynolds  stresses,  and  all  of  the  triple  products.  In  the  core  of  the  vortex  there  was 
very  little  mixing.  They  felt  that  the  high  axial  velocity  in  the  core,  in  combination 
with  vortex  wandering,  could  account  for  all  the  observed  shear  stress,  u'v^,  and 
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most  of  the  observed  normal  stresses,  u/2,  t^2,  and  v(?.  As  the  vortices  interacted 
with  the  boundary  layer,  secondary  vortices  of  opposite  sign  first  formed  along  the 
surface,  then  broke  away  and  were  entrained  into  the  primary  vortex.  They  felt 
that  this  breaking  away  process  caused  the  observed  large  levels  of  u'vj1,  u'2,  and 
^2. 

Mehta  and  Cantwell  (1987)  studied  the  mean  flow  and  turbulence  structure  of 
a  half-delta  wing  in  air  in  the  absence  of  interactions  with  other  surfaces  or  shear 
layers.  They  felt  this  was  needed  as  a  baseline  for  comparison  with  other  interaction 
studies  where  a  similar  type  vortex  was  used.  They  concluded  that  the  vortex  was 
fully  developed  in  a  distance  of  3  vortex  generator  heights,  because  the  vorticity 
and  circulation  had  reached  constant  values,  and  because  of  the  coincident  location 
of  the  maximum  vorticity  and  maximum  velocity  deficit.  They  measured  and  re¬ 
ported  contours  for  all  the  Reynolds  stresses  except  v'w'.  The  turbulence  was  nearly 
isotropic  as  indicated  by  nearly  equal  values  and  time  histories  of  the  three  normal 
Reynolds  stresses.  While  contours  of  the  secondary  shear  stress  were  complicated  in 
form,  the  authors  felt  they  were  similar  in  nature  to  the  results  reported  by  others. 

We8tphal  and  Mehta  (1987)  looked  at  the  effect  of  lateral  oscillation  of  a  stream- 
wise  vortex,  when  embedded  in  a  turbulent  boundary  layer.  Their  primary  objective 
was  to  see  if  the  flattening  in  the  vortex  core,  seen  in  previous  work,  could  be  at¬ 
tributed  to  vortex  meander.  They  observed  that  a  lateral  (spanwise)  oscillation  of 
the  vortex  generator  caused  a  flattening  of  the  vorticity  contours  when  compared 
with  vorticity  contours  for  a  stationary  generator.  They  also  saw  some  changes  in 
the  Reynolds  stresses,  particularly  v'2  and  u'w1 .  While  they  concluded  that  vortex 
meander  will  cause  flattening  in  the  vortex  core,  they  could  not  say  whether  or  not 
this  was  the  cause  of  flattening  in  previous  studies. 

Mehta  and  Bradshaw  (1988)  continued  their  series  of  studies  in  a  detailed  report 
on  a  vortex  pair  with  “common  flow  upward”.  With  this  configuration,  the  pair 
of  vortices  tended  to  rise  to  a  height  above  the  surface  equal  to  approximately 


twice  the  boundary-layer  thickness.  While  the  distribution  of  u'u'  was  indicative 
of  “convection  of  pre-existing  boundary-layer  turbulence  by  the  vortices”,  v'v'  and 
w'w'  showed  symmetrical  maxima  in  the  core  regions.  In  general,  turbulent  kinetic 
energy  transport  appeared  to  occur  from  the  high  energy  core  region  to  surrounding 
regions.  Modification  of  the  Reynolds  shear  stresses  appeared  to  occur  by  rotation 
of  the  fluid  around  the  vortex  core.  In  general,  they  concluded  that  the  vorticies 
modify  but  do  not  destroy  the  anisotropy  of  the  turbulence. 

Pauley  and  Eaton  (1988)  investigated  the  development  of  longitudinal  vortex 
pairs  embedded  in  a  turbulent  boundary  layer.  Following  their  earlier  work  with 
Westphal  et  al.,  they  sought  to  characterize  how  the  history  and  interaction  of  vortex 
pairs  varied  with  pair  configuration,  spacing,  and  strength.  Configurations  studied 
were  counter-rotating  with  “common  flow  up”  and  “common  flow  down”  and  co¬ 
rotating.  Their  results  confirmed  earlier  predictions  by  Pearcey  that  interaction 
with  the  boundary  layer  would  result  in  rapid  diffusion  of  vorticity  accompanied 
by  rapid  vortex  growth.  While  Westphal  et  al.  found  this  was  true  for  a  single 
embedded  vortex,  Pauley  and  Eaton  found  an  even  stronger  effect  when  vortex 
pairs  were  placed  close  enough  to  influence  each  other.  This  effect  was  the  most 
pronounced  for  counter-rotating  vortices  with  “common  flow  up” .  In  contrast  to  the 
peak  vorticity,  the  total  circulation  of  t.hp  primary  vortex  (obtained  by  integration) 
decreased  much  more  slowly.  The  primary  mechanism  for  this  was  skin  friction  at 
the  wall,  which  was  the  strongest  in  the  “common  flow  down”  case  where  the  vortices 
were  the  closest  to  the  surface.  One  final  observation  regarding  the  “common  flow 
down”  configuration  was  that  varying  the  strength  of  the  vortices  did  not  change  ihe 
amount  of  boundary  layer  thinning  but  rather  the  width  of  the  thinning  region,  due 
to  compensating  effects  from  the  stronger  image  vortices.  With  respect  to  turbulent 
structure,  they  discovered  that  a  velocity  deficit  in  the  vortex  core  was  associated 
with  an  increase  in  turbulent  kinetic  energy  and  that  the  upwash  region  exhibits 
high  turbulent  kinetic  energy  production.  Finally,  they  found  they  could  classify 
the  type  of  strain  by  plotting  the  data  in  terms  of  the  invariants  of  the  anisotropic 
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Reynolds  stress  tensor.  Various  Reynolds  stress  models  were  tested  against  their 
data. 

3.2  Unsteady,  Separating  Turbulent  Flows 

Perhaps  the  appropriate  starting  point  for  unsteady  flow  is  with  the  two  classic 
problems  posed  by  Stokes  (1851).  While  the  analysis  was  for  laminar  flow,  the 
analytical  solution  for  each  problem  is  straight-forward  and  each  solution  yields  a 
characteristic  length  which  is  very  useful  in  parameterization  of  unsteady,  turbulent 
flows.  Stokes’  first  problem  was  a  flat  plate  started  abruptly  from  rest,  for  which 
the  solution  is: 


The  equation  above  defines  a  characteristic  thickness,  t*3  =  2 \fui.  Stokes’  second 
problem  was  a  flat  plate  subjected  to  a  sinusoidal  freestream  oscillation.  The  solu¬ 
tion  is: 

w0=Mi)c  <3-2> 

where  the  characteristic  thickness,  ls  =  (2 rz/H)1/2,  becomes  very  small  at  high 
frequencies. 

Karisson  (1959)  performed  one  of  the  first  experimental  studies  of  an  unsteady 
turbulent  boundary  layer.  He  used  a  wind  tunnel  with  a  sinusoidally  driven  shutter 
system  to  produce  a  sinusoidal  variation  superimposed  on  a  mean  freestream  veloc¬ 
ity.  He  characterized  the  velocity  response  in  terms  of  a  sum  of  the  mean,  the  first 
harmonic,  and  a  residual  containing  higher  harmonics  and  the  chaotic  portion  of 
the  turbulence.  He  observed  no  systematic  variation  of  the  mean  velocity  profiles 
with  amplitude  or  frequency,  even  for  amplitudes  as  high  as  34%.  In  the  inner 
portion  of  the  boundary  layer,  the  velocity  led  the  freestream  in  phase  by  as  much 
as  35°,  while  at  low  frequencies  there  was  a  lag  in  phase  in  the  outer  portion  of  the 
boundary  layer. 
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Cousteix,  Desopper,  and  Houdeville  (1977)  studied  the  structure  and  develop¬ 
ment  of  a  turbulent  boundary  layer  in  an  oscillatory  flow.  They  measured  the  mean 
velocity,  u'2,  and  uV  for  comparison  with  the  prediction  from  two  different  tur¬ 
bulence  models,  a  mixing-length  model,  and  a  transport  equation  model.  While 
their  tests  were  limited  to  low  Strouhal  numbers  (0.0127  based  on  boundary  layer 
thickness),  they  concluded  that  the  general  behavior  of  the  boundary  layer  and 
the  structure  of  the  turbulence  were  not  affected  by  the  flow  oscillations.  Further¬ 
more,  they  concluded  that,  within  the  limits  of  low  Strouhal  number  flows,  both 
the  mixing-length  model  and  the  k  —  e  transport  equation  model  did  a  good  job  of 
predicting  the  behavior  of  the  turbulent  boundary  layer. 

Cousteix,  Javelle,  and  Houdeville  (1981)  extended  the  work  described  above  to  a 
broader  range  of  Strouhal  numbers.  They  studied  a  Strouhal  range  from  1.5  to  18, 
where  the  Strouhal  number  was  based  on  the  distance  from  the  virtual  origin.  They 
observed  a  periodicity  in  the  mean  velocity  that  appeared  to  be  the  consequence 
of  a  travelling  wave  with  a  convective  velocity  of  0.75Uoo.  A  similar  periodicity  in 
the  phase  averaged  turbulence  intensity  inferred  a  convective  velocity  of  0.851/qo. 
The  shear  stress  correlation  coefficient  maintained  a  constant  value  of  0.45  across 
the  boundary  layer,  consistent  with  observations  from  steady  flow.  Across  their 
range  of  Strouhal  numbers,  they  concluded  that  the  mean  flow  was  not  affected 
by  the  imposed  unsteadiness.  However,  for  Strouhal  numbers  greater  than  5,  the 
inner  region  law-of-the-wall  did  not  hold  and  did  not  have  a  consistent  logarithmic 
overlap  area  with  the  outer  region  of  the  boundary  layer. 

Simpson,  Chew,  and  Shivaprasad  (1981)  characterized  the  nature  of  the  unsteady 
separating  turbulent  boundary  layer  in  a  series  of  reports.  (Also  see  Shiloh  et  al). 
Their  experimental  work  was  done  in  an  air  tunnel  using  a  programmable  rotating 
blade  damper  to  create  a  periodic  freestream  unsteadiness.  In  a  sequence  of  three 
reports,  they  documented  all  the  mean  velocities,  all  the  Reynolds  stresses,  and  the 
flatness  and  skewness  factors  for  a  separating  boundary  layer. 
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Jayaraman,  Parikh  and  Reynolds  (1982)  studied  the  effects  of  a  steady,  uniform 
upstream  velocity  with  a  sinusoidally  varying  downstream  velocity.  They  observed 
that  while  the  mean  flow  was  relatively  unaffected  by  the  unsteadiness,  the  periodic 
velocity  field  was  strongly  affected  and  the  law  of  the  wall  could  not  be  applied 
to  phase-averaged  profiles.  They  concluded  that  the  Strouhal  number,  based  on 
the  streamwise  distance  in  the  unsteady  region,  was  the  dominant  parameter  in 
describing  the  effects  of  freestream  oscillation.  In  general,  at  low  Strouhal  numbers, 
the  phase  of  the  boundary  layer  tended  to  lag  the  freestream,  with  the  inner  region 
of  the  boundary  layer  leading  the  outer  region.  At  high  Strouhal  numbers,  the  phase 
variation  was  confined  to  a  small  “Stokes  layer”  near  the  surface,  consistent  with 
the  classic  solution  to  Stokes  Second  Problem.  Finally,  the  integral  parameters, 
displacement  thickness  and  shape  factor,  showed  periodic  variations.  While  the 
amplitude  of  oscillation  decreased  with  increasing  frequency,  at  high  frequencies 
the  phase  of  the  displacement  thickness  was  nearly  180°  out  with  respect  to  the  free 
stream. 

Simpson,  Shivaprasad,  and  Chew  (1983)  continued  their  study  of  the  separating 
turbulent  boundary  layer.  They  presented  time-averaged  and  phase-conditioned 
profiles  of  the  streamwise  velocity,  the  back-flow  velocity,  and  the  Reynolds  stress 
components  u'2,  v'2,  and  —u'v'.  They  stressed  the  importance  of  characterizing 
the  flow  in  terms  of  the  intermittency  of  forward  motion,  7^,  and  they  presented  a 
number  of  correlations  between  7^  and  the  phase  conditioned  shape  factor,  H.  The 
same  correlations  held  for  unsteady  flow  as  for  steady  flow  except  that  the  unsteady 
cycle  produced  a  hysteresis  loop  which  increased  in  size  with  increasing  frequency. 
Finally,  they  observed  that  the  near-wall  velocity  led  the  pressure  gradient  by  as 
much  as  135°. 

Binder  and  Kueny  (1981)  and  then  Binder  et  al.  (1985)  studied  large  amplitude 
oscillations  in  the  wall  region  of  a  turbulent  channel  flow.  Their  experiments  in  a 
water  facility  varied  the  Reynolds  number  by  a  factor  of  2,  the  frequency  by  a  factor 
of  15,  and  the  amplitude  from  5%  to  70%  of  the  centerline  velocity.  This  included 
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flow  regimes  where  reverse  flow  occurred  at  the  surface  for  portions  of  each  cycle. 
They  concluded  that  as  long  as  the  frequency  was  sufficiently  high  such  that  the 
the  Stokes  layer  was  within  the  viscous  sublayer,  then  the  phase  and  magnitude  of 
the  velocity  and  wall  stress  followed  the  Stokes  solution.  This  restraint  required  the 
Stokes  thickness,  i3 ,  to  be  such  that 
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In  this  region,  the  near  wall  velocity  fluctuations  ranged  from  0  to  over  the 
distance  of  the  Stokes  thickness,  and  the  phase  led  the  center  line  phase,  approaching 
45°  at  the  wall.  For  intermediate  values,  10  <  If  <  20,  the  behavior  began  to 
deviate  from  the  Stokes  solution,  and  for  If  >  20,  the  behavior  was  more  quasi¬ 
steady  in  nature.  They  also  concluded  that  their  data  confirmed  their  criteria  for 
predicting  reverse  flow  (1981), 
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Finally,  even  though  in  some  cases  the  deterministic  fluctuations,  u'2,  were  25  times 
larger  than  the  mean  fluctuations,  u'2,  there  were  no  significant  variations  in  the 
mean  statistical  quantities  when  compared  to  steady  flow. 

Most  recently,  Tardu,  Binder,  and  Blackwelder  (1987)  reported  their  preliminary 
results  of  a  study  to  understand  the  relationship  between  turbulent  bursting  and 
imposed  unsteadiness.  While  they  had  seen  some  trends,  their  observations  had  not 
yet  led  to  any  clear  quantitative  conclusions. 

Chehroudi  and  Simpson  (1985)  continued  the  study  of  Simpson  et  al.  in  sep¬ 
arating  turbulent  boundary  layers  using  a  rapidly  scanning  laser  anemometer.  In 
this  technique,  the  LDA  measuring  volume  was  rapidly  scanned  across  the  entire 
boundary  layer  by  a  steering  mirror  which  oscillated  at  frequencies  up  to  59  Hz. 
In  this  manner  they  were  able  to  obtain  space-time  correlation  coefficients,  from 
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which  they  could  obtain  an  integral  length  scale.  They  also  produced  scatter  plots 
of  the  velocity  within  a  spatial  interval  to  infer  the  percent  of  time  that  the  velocity 
exceeded  the  mean  value.  They  concluded  that  the  structure  of  the  flow  in  the 
backflow  region  was  strongly  related  to  the  local  outer  region  flow. 

Brereton,  Reynolds,  and  Carr  (1985)  studied  the  reponse  of  the  turbulent  bound¬ 
ary  layer  to  abrupt  changes  in  the  free-stream  velocity.  While  maintaining  a  steady 
upstream  velocity  in  their  water  facility,  they  abruptly  diverted  a  portion  of  the 
freestream  flow,  causing  a  streamwise  linear  decrease  of  the  freestream  velocity. 
Their  results  agreed  well  with  the  analytical  solution  to  Stokes’  first  problem,  where 
only  a  Stokes  layer,  of  thickness  t*  —  2 y/ut,  followed  the  “motion”  of  the  stationary 
plate,  and  the  remainder  of  the  flow  acted  as  a  slug  with  the  freestream.  This  was 
true  for  both  u/Z7<x>  and  u'2/!/^  when  scaled  with  y/t During  the  first-half  second 
of  the  transient,  the  displacement  thickness  of  the  entire  flow  responded  simulta¬ 
neously.  In  contrast,  during  relaxation  after  the  first  half  second,  the  displacement 
thickness  responded  on  a  time  scale  proportional  to  the  “local  freestream  time  of 
flight”. 

Cousteix  and  Houdeville  (1985)  continued  their  study  of  the  evolution  of  turbu¬ 
lence  and  skin  friction  in  an  oscillating  boundary  layer.  The  entire  flow  was  oscil¬ 
lated  such  that  the  amount  of  variation  of  the  freestream  velocity  (12  to  15  percent) 
was  maximum  near  the  leading  edge  and  then  decreased  linearly  with  streamwise 
position,  x.  While  the  time-averaged  properties  were  not  affected  much  by  the  un¬ 
steadiness,  they  observed  a  “quasi-periodic  spatial  oscillation”  of  the  momentum 
thickness  and  displacement  thickness.  They  related  this  effect  to  an  interaction 
between  the  unsteadiness  and  the  convection  of  turbulence  at  a  celerity  of  0.85f7oo. 
For  y+  values  on  the  order  of  10,  a^d  outside  the  boundary  layer,  u'2  was  in  phase 
with  tt,  but  in  the  middle  region  of  the  boundary  layer  the  phase  varied  up  to  180°. 
For  high  frequencies,  these  variations  were  confined  to  a  thin  Stokes  layer  on  the 
order  of  is  =  ( 2i//f2 W2.  While  the  time-averaged  profiles  exhibited  a  logarithmic 


region  over  the  entire  frequency  range,  only  at  low  frequencies  w as  the  friction  ve¬ 
locity  from  the  law-of-the-wall  related  to  the  actual  skin  friction  measured.  The 
amplitude  and  phase  of  the  skin  friction  did  not  correlate  well  with  flx/t/oo  hut  did 
correlate  well  with  Clu/(u £).  Finally,  they  saw  no  obvious  correlation  between  the 
unsteadiness  and  any  bursting  phenomena. 

Cook,  Murphy,  and  Owen  (1985)  and  then  Cook  (1986)  reported  on  turbulent 
boundary  layers  in  oscillating  flows.  The  experimental  work  was  performed  in  a  low- 
speed  air  tunnel  where  the  flow  rate  was  controlled  by  a  choked,  variable-throat- 
area  exit  nozzle.  They  compared  the  experimental  results  to  analytical  results 
using  two  different  turbulence  models,  one  developed  by  Cebeci-Smith  and  the 
other  by  Glushko.  While  the  results  from  neither  model  matched  the  experimental 
observations  very  well,  in  general,  the  results  from  the  Cebeci-Smith  model  yielded 
results  which  were  closer.  In  the  experiments,  the  amplitude  of  oscillation  and 
the  phase  difference  of  the  freestream  both  increased  with  increasing  frequency. 
While  the  amplitude  variation  was  greatest  for  low  values  of  x,  the  streamwise 
distance  from  the  entrance,  the  phase  variation  was  greatest  at  large  values  of  x. 
Inside  the  boundary  layer,  the  ratio  of  local  velocity  to  freestream  often  exceeded 
unity.  Meanwhile,  the  phase  of  the  inner  boundary  layer  region  tended  to  lead  the 
freestream  while  the  phase  of  the  outer  region  tended  to  lag  the  freestream.  Both  of 
these  trends  were  most  pronounced  at  larger  values  of  x  and  for  higher  frequencies. 
Finally,  the  turbulence  intensity,  \/(u/)2/I700,  deviated  from  steady  predictions  in 
the  outer  region  of  the  boundary  layer,  particularly  at  higher  frequencies. 

In  their  review  paper,  Reynolds  and  Carr  (1985)  proposed  a  simple  control 
volume  vorticity  balance  for  a  boundary  layer  subjected  to  an  unsteady  adverse 
pressure-gradient.  The  accumulation  of  spanwise  vorticity  is  seen  as  the  difference 
between  the  convective  inflow  of  vorticity  and  the  convective  outflow  of  vorticity, 
plus  the  vorticity  flux  at  the  surface  due  to  the  pressure  gradient.  While  the  span- 
wise  vorticity  embedded  in  the  boundary  layer  is  convected  downstream  by  the 
mean  flow,  there  is  a  tendency  for  the  vortices  to  travel  upstream  due  to  an  image 
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effect  from  the  surface.  When  the  boundary  layer  is  subjected  to  a  sudden  ad¬ 
verse  pressure  gradient,  three  things  occur:  (1)  the  image  effect  may  now  allow  the 
spanwise  vorticity  to  overcome  convection  and  move  upstream,  (2)  the  larger  flux 
of  vorticity  into  the  control  volume  than  out  causes  an  accumulation  of  vorticity, 
and  (3)  the  adverse  pressure  gradient  causes  a  flux  of  opposite  sign  vorticity  into 
the  control  volume.  Since  the  third  effect  tends  to  be  small,  the  vorticity  tends  to 
“grow  to  the  scale  of  the  boundary  layer  and  then  is  shed.” 

Cook  and  Murphy  (1987)  used  the  facility  described  previously  to  study  in  more 
detail  the  effect  of  amplitude  and  frequency  of  oscillation  on  the  phase-averaged, 
streamwise  normal  stress,  u/2.  The  amplitude  of  oscillation  ranged  from  6.5%  to 
16.9%  and  the  frequency  ranged  from  3.85 Hz  to  21Hz,  resulting  in  local  Strouhal 
numbers  ranging  from  0.54  to  4.05.  While  there  was  no  effect  of  amplitude  or 
frequency  on  the  mean  turbulence  levels,  they  observed  a  phase  shift  between  u 
and  Uqo  and  a  phase  shift  between  u/2  and  u.  Both  phase  shift  effects  varied  with 
frequency  and  amplitude  of  oscillation,  and  with  streamwise  position.  The  frequency 
where  effects  were  seen  was  more  than  an  order  of  magnitude  below  the  estimated 
turbulent  burst  frequencies.  Finally,  their  results  indicated  that  the  variation  in  the 
turbulence  level  was  more  than  could  be  accounted  for  by  the  variation  in  boundary 
layer  thickness  alone,  in  disagreement  with  the  suggestion  by  Cousteix  et  al.  (1985). 

Brereton  and  Reynolds  (1987)  expanded  the  work  of  Jayaraman  et  al.  with  two- 
component  LDA  measurements.  They  observed  that  most  time-averaged  quantities, 
including  turbulent  length  scales  and  the  dominant  production  term,  were  similar 
to  their  counterparts  in  steady  flow  and  were  invariant  with  frequency.  The  two 
time  averaged  terms  that  did  vary  with  frequency  were  the  d(ufcu,)/dxj.  tensor 
and  the  u'v^du/dy)  production  term.  In  contrast,  most  of  the  phase-conditioned 
quantities  (except  for  v)  were  strongly  dependent  on  frequency.  The  production  of 
kinetic  energy  was  very  anisotropic  in  that  u'u'  was  much  greater  them  v'v'  and  that 
u'u1  lead  v'v'  in  phase  angle.  Finally,  in  all  cases,  they  found  that  transfer  of  the 
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kinetic  energy  was  from  the  organized  unsteady  component  uti  to  the  unorganized 
component  u'u1. 

Gajdeczko  and  Ramaprian  (1987)  also  studied  a  periodic  turbulent  boundary 
layer  in  an  advert  pressure  gradient.  In  theix  wr.ter  tunnel,  the  wall  opposite  lIik 
test  surface  was  contoured  to  create  the  desired  adverse  pressure  gradient,  and  the 
periodic  flow  was  caused  by  a  rotating  sleeve  in  the  exit  section  of  the  channel.  Their 
2-component  LDA  provided  mean  and  phase-averaged  profiles.  The  time-averaged 
profiles  showed  the  effects  of  unsteadiness  in  the  strong  pressure  gradient  case,  in 
the  regions  where  flow  reversal  occurred.  In  their  phase-averaged  profiles,  the  shape 
of  the  profiles  changed  substantially  as  the  Strouhal  number,  based  on  streamwise 
position,  increased  from  1.89  to  28.43.  They  saw  poor  agreement  with  predictions 
from  a  computational  model  which  was  developed  for  unsteady  flow  with  a  zero 
pressure  gradient.  Thus  they  concluded  that  turbulent  production  occurs  much 
further  from  the  wall  in  flows  with  an  adverse  pressure  gradient. 

3.3  Vorticity  in  Unsteady,  Separating,  Turbulent  Flows 

The  work  presently  published  on  vorticity  in  unsteady  separating  turbulent  flows 
relates  mostly  to  spanwise  vortices  which  are  formed  during  dynamic  pitching  op¬ 
erations  and  those  that  are  shed  when  the  boundary  layer  separates.  Very  little, 
if  any,  work  has  been  done  with  skew-induced  streamwise  vorticity  in  an  unsteady, 
turbulent  boundary  layer. 

McCroskey,  Carr,  and  McAlister  (1976)  studied  dynamic  stall  on  oscillating  air¬ 
foils.  They  sought  to  understand  the  mechanism  that  led  to  the  shedding  of  the 
strong  “vortex-like”  structure  during  dynamic  stall.  In  all  their  study  cases,  they 
observed  that  the  difference  between  dynamic  and  static  stall  is  that  in  dynamic 
stall  a  “vortex-like  disurbance  develops  and  passes  over  the  airfoil  at  a  convection 
velocity  on  the  order  of  g  to  ^  of  U00.'n  They  concluded  that  dynamic  stall  orig¬ 
inates  with  a  breakdown  of  the  turbulent  boundary  layer,  not  with  bursting  of  a 


25 


laminar  separation  bubble.  As  the  frequency  was  varied,  the  phase  within  each 
cycle  changed  for  the  different  events  associated  with  the  formation  and  convection 
of  the  disturbance  vortex,  but  the  basic  progression  of  events  did  not  change.  In 
a  similar  way,  the  reduced  frequency  affected  both  the  strength  and  phase  of  the 
vortex  shedding.  In  contrast  to  this,  there  was  very  little  effect  of  Reynolds  number 
on  the  nature  of  the  dynamic  stall. 

Reisenthel,  Koga,  and  Nagib  (1985)  looked  at  control  of  separated  flows  with  a 
spanwise  flap.  The  flap  was  actuated  in  a  manner  that  generated  a  large,  spanwise 
vortical  structure  which  generally  caused  earlier  reattachment  and  better  pressure 
recovery.  Measurements  of  the  pressure  coefficient,  Cp,  indicated  a  sharp,  local  drop 
in  Cp  at  the  flap  with  a  much  quicker  recovery  to  a  Cp  of  0.  This  effect  increased  as 
the  frequency  increased  from  2.5 Hz  to  1 0Hz.  It  also  increased  as  the  flap  deflection 
angle  increased  to  40°,  after  which  no  further  increase  was  seen  for  larger  deflection 
angles.  These  various  effects  collapsed  into  a  single  curve  when  the  reattachment 
length  was  plotted  against  the  reduced  frequency,  fh/Uoo,  where  h  was  the  height 
of  the  flap  measured  normal  to  the  wall.  In  summary,  the  reattachment  length 
decreased  linearly  with  increasing  reduced  frequency  until  it  reached  a  value  of  0.4, 
after  which  no  further  improvement  was  seen. 

Sigurdson  and  Roshko  (1985)  looked  at  unsteady  excitation  of  the  reattaching 
flow  on  a  circular  cylinder.  A  closed  circular  cylinder  was  placed  in  a  wind  tunnel 
with  its  axis  aligned  parallel  with  the  freestream  flow.  Separation  occured  at  the 
circular  leading  edge  of  the  cylinder.  The  flow  was  excited  through  slits  at  the  lead¬ 
ing  edge  by  an  acoustic  speaker  placed  inside  the  cylinder.  The  pressure  coefficient, 
the  reattachment  length,  and  the  drag  were  used  as  global  evaluation  parameters. 
All  three  parameters  were  reduced  with  forcing.  The  maximum  effect  was  seen  for 
frequencies  that  corresponded  to  a  wavelength  comparable  to  the  separation  bub¬ 
ble  height.  The  optimum  excitation  frequency,  when  non-dimensionalized  with  the 
Cylinder  diameter,  fexD/Uoo,  was  in  the  range  of  2  to  3.  Alternately,  this  same 
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frequency  range,  when  scaled  with  the  local  velocity  at  separation  and  the  bubble 
height,  fexh/Us,  was  0.26  to  0.32. 

Luttges,  Robinson,  and  Kennedy  (1985)  looked  at  the  control  of  unsteady,  sepa¬ 
rated  flow  structures  on  airfoils  using  a  spanwise  vortex  generator  and  a  spanwhe 
air  jet.  They  first  observed  the  structures  that  formed  on  a  pitching  airfoil  with  no 
control  devices.  During  pitching,  a  vortex  formed  near  the  leading  edge.  It  con- 
vected  downstream  and  then  elicited  a  trailing  edge  vortex  of  opposite  sign  which 
led  to  a  “cataclysmic  stall” .  As  the  mean  angle  of  attack  was  increased,  the  leading 
edge  vortex  formed  earlier  and  convected  more  slowly.  In  contrast,  increasing  the 
pitching  amplitude  or  reduced  frequency  delayed  the  vortex  formation  and  led  to 
a  faster  convection  velocity.  Similar  effects  were  seen  on  a  stationary  airfoil  with  a 
vortex  initiated  by  either  a  vortex  generator  or  an  air  jet.  Since  the  air  jet  provided 
better  flexibility  and  dynamic  response,  it  was  used  for  most  of  the  experiments.  Fi¬ 
nally,  the  air  pulses  were  combined  with  the  pitching  of  the  airfoil.  The  effect  of  the 
air  pulses  varied  with  reduced  frequency  and  phasing,  but  for  reduced  frequencies 
greater  than  1.0  (meaning  there  was  always  at  least  one  vortex  over  the  airfoil)  flow 
adherence  to  the  airfoil  was  improved.  They  concluded  that  careful  programming 
of  air  pulses  could  eliminate  flow  conditions  that  lead  to  “cataclysmic  stall” . 

Lutges  (1987)  continued  his  work  with  the  formation  and  convection  of  leading 
edge  and  trailing  edge  vortices  during  unsteady  flow.  While  part  of  the  study  con¬ 
centrated  on  characterizing  the  vortices  on  a  flat  plate  and  a  NACA  0015  airfoil,  the 
other  portion  of  the  study  concentrated  on  understanding  the  unsteady  mechanisms 
that  allow  a  dragonfly  to  produce  lift  coefficients  on  the  order  of  5.  For  his  studies 
with  the  flat  plate  and  airfoil,  he  looked  at  mean  angles  of  attack  ranging  from 
10°  to  30°,  at  oscillating  angles  ranging  from  2°  to  9°,  and  at  reduced  frequencies 
ranging  from  0.4  to  2.0.  He  concluded  that  while  the  reduced  frequency  was  very 
significant  in  describing  the  results,  the  Reynolds  number  and  surface  geometry 
were  relatively  inconsequential. 
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Bhattacharjee,  Scheelke,  and  Troutt  (1985)  looked  at  modification  of  vortex  in¬ 
teractions  in  a  reattaching  separated  flow.  Forcing  was  provided  by  a  loud  speaker 
placed  on  the  wind  tunnel  wall  opposite  the  backward  facing  step  that  caused  the 
separation.  The  forcing,  as  observed  with  hot-w;re  probe  signals,  appeared  to  cause 
pairing  and  a  trend  toward  two-dimensional  flow  of  spanwise  vortical  structures. 
Forcing  modified  primarily  the  outer  portion  of  the  boundary  layer,  moving  higher 
velocity  fluid  closer  to  the  surface.  The  forcing  seemed  to  be  optimum  at  a  fre¬ 
quency  corresponding  to  the  passage  frequency  of  a  mixing  layer  of  the  width  of 
the  backward  facing  step.  This  corresponded  to  a  Strouha!  number  range  based  on 
the  step  height  of  0.2  to  0.4.  The  forcing  reduced  the  reattachmeni  length  by  about 
10%. 

In  a  similar  program,  Roos  and  Kegelman  (1985)  looked  at  control  of  coherent 
structures  in  reattaching  laminar  and  turbulent  shear  layers.  They  used  a  flap  which 
projected  along  the  flow  direction  from  the  edge  of  a  backward  facing  step.  The 
flap  was  oscillated  at  very  small  amplitudes  on  the  order  of  1%  of  the  unseparated 
boundary  layer  thickness,  and  at  several  frequencies  in  the  range  of  20 Hz.  This 
corresponded  to  a  Strouhal  number,  based  on  the  backstep  height,  of  0.22.  Gentle 
excitation  of  the  shear  layer  at  the  separation  point  provided  and  strong  “regular¬ 
izing”  process  and  reduced  the  reattachment  length  by  nearly  20%  in  laminar  flow 
and  nearly  30%  in  turbulent  flow.  While  the  changes  due  to  excitation  were  small 
in  the  mean  velocity  profile,  there  was  a  substantial  increase  in  the  Reynolds  shear 
stress  and  pressure  fluctuations,  indicating  intensified  turbulence  activity  within 
and  beneath  the  reattaching  shear  layer. 

Walker,  Helin,  and  Chou  (1985)  studied  unsteady  surface  pressure  measurements 
on  a  pitching  airfoil.  The  pitch  rates  on  the  NACA  0015  airfoil  ranged  from  230 ° /sec 
to  1380°/sec,  while  the  angle  of  attack  varied  from  0°  to  60°.  With  a  6m  chord  and 
a  range  of  freestream  velocities  from  20 ft/ sec  to  80 ft/ sec,  the  resulting  range  of 
dimensionless  pitch  rates  was  0.05  to  0.6.  Comparison  of  the  pressure  distribution 
and  the  integrated  lift  coefficients  revealed  parallel  and  consistent  trends.  For  a 
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constant  Reynolds  number,  increasing  the  pitch  rate  increased  ,  dC^da,  and 
the  angle  of  attack  at  .  Increasing  the  Reynolds  Number  at  constant  pitch 

rate  produced  the  opposite  effect.  The  characteristics  of  the  curves  were  similar 
when  plotted  for  constant  non-dimensional  pitch  rate,  a+  =  dec/U^. 

In  a  computational  study,  McCoy  and  Sarpkaya  (1986)  looked  at  the  stability 
of  vortex  motion  in  an  oscillating  flow.  Their  study  was  geared  towards  the  un¬ 
derstanding  of  cycle  to  cycle  variations  in  the  flow  observed  around  a  cylinder  in 
unsteady  flow.  They  traced  the  separate  histories  of  either  one  or  two  vortices 
placed  in  the  vicinity  of  a  circular  cylinder.  Their  conclusion  was  that  observed 
variations  in  the  vortex  motions  resulted  mostly  from  chaotic  motion  caused  by  the 
advection  of  vortices  in  a  time-dependent  flow. 

Gad-el-Hak  and  Ho  (1986)  studied  the  unsteady  vortical  flow  around  three- 
dimensional  lifting  surfaces.  They  used  flow  visualization  with  seven  different  lifting 
surfaces  in  a  tow-tank  to  characterize  the  formation,  convection,  interaction,  and 
shedding  of  spanwise  vortices  on  a  pitching  lifting  surface.  The  reduced  frequency, 
the  leading  edge  shape,  and  the  planform  of  the  lifting  surface  all  had  significant 
effects  on  the  nature  of  the  vortex  interactions. 

Gad-el-Hak  and  Blackwelder  (1987)  studied  control  of  the  discrete  vortices  from 
a  delta  wing,  using  oscillatory  injection  of  fluid  at  the  leading  edge.  They  used  both 
flow  visualization  and  hot  wire/film  in  both  a  water  tow  tank  and  a  wind  tunnel. 
The  spanwise  vortices  associated  with  dynamic  lift  were  significantly  altered  by  the 
perturbation  at  the  leading  edge.  The  effect  of  the  perturbation  was  the  largest 
when  the  speed  of  the  injected  fluid  was  near  the  ambient  velocity  and  when  the 
excitation  frequency  was  a  subharmonic  of  the  natural  shedding  frequency  of  the 
delta  wing. 
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3.4  Unanswered  Questions  in  the  Literature. 

One  of  the  most  obvious  questions,  which  appears  not  to  have  been  addressed  in 
the  literature,  is  addressed  by  this  study.  That  is,  how  do  vortex  generators  affect 
the  nature  of  a  boundary  layer  subjected  to  an  unsteady  adverse  pressure  gradient? 
To  answer  that  question,  one  must  address  many  more  specific  questions,  several  of 
which  are  posed  here. 

(1)  The  studies  of  unsteady  flows  have  shown  that  there  is  a  thin  Stokes  layer 
that  follows  the  pressure  gradient,  while  most  of  the  flow  acts  as  a  slug  flow  with 
the  freestream.  In  contrast,  vortex  generators  tend  to  act  in  the  wake  region  of  the 
boundary  layer.  How  much  interaction  will  there  be  between  these  two  processes, 
and  over  what  portion  of  the  boundary  layer  will  this  interaction  occur? 

(2)  Studies  of  streamwise  vortices  in  a  steady  adverse  pressure  gradient  indicate 
that  the  pressure  gradient  accelerates  vortex  growth,  which  decrease  the  core  vor- 
ticity.  What  will  happen  to  the  streamwise  vorticity  when  subjected  to  an  unsteady 
adverse  pressure  gradient?  If  the  streamwise  vorticity  is  rapidly  diffused,  can  it  still 
help  delay  boundary-layer  separation? 

(3)  Many  of  the  vortex  studies  emphasize  that  the  mechanisms  associated  with 
vortex  generators  are  dominated  by  convection  effects.  If  this  is  the  case,  what  role 
do  the  turbulent  terms  play  in  the  mechanisms  that  enable  a  vortex  generator  to 
delay  boundary-layer  separation? 

(4)  Finally,  what  are  the  characteristic  time  scales  associated  with  the  response 
of  the  boundary  layer,  with  and  without  vortex  generators,  to  an  unsteady  adverse 
pressure  gradient? 

These  are  the  questions  that  this  study  will  seek  to  answer. 


30 


4.  Experimental  Facility 


The  Unsteady  Flow  Facility  used  for  this  program  is  a  complex,  unsteady,  water 
channel  that  has  evolved  through  the  experimental  research  programs  of  Jayaraman, 
Parikh,  and  Reynolds  (1982)  and  Brereton  and  Reynolds  (1987).  During  the  period 
from  1985  to  1987,  major  upgrades  to  the  facility  were  designed  and  implemented  as 
a  cooperative  effort  among  three  graduate  students  and  Professor  Reynolds.  Further 
details  of  the  mechanical  design  of  the  test  section  and  computer-controlled  valves 
are  described  by  Henk  and  Reynolds  (1989).  In  a  similar  manner,  Carlson  and 
Reynolds  (1989)  document  in  more  detail  the  computerized  valve  control  system 
and  data  acquisition  system. 

As  shown  in  figure  4.1,  water  is  pumped  from  a  sump  tank  at  atmospheric  pres¬ 
sure  to  a  constant  head  tank  overhead,  where  the  water  level  is  maintained  by  a  weir 
at  2.7m  above  the  elevation  of  the  leading  edge  of  the  test  surface.  Water  exits  the 
constant  head  tank  and  enters  the  test  section  through  a  series  of  fiow-straightenmg 
devices  and  a  20:1  contraction.  A  new  boundary  layer  forms  on  the  upper  surface 
as  the  flow  passes  over  an  elliptical  leading  edge.  For  turbulent  experiments,  the 
boundary  layer  is  then  tripped  at  a  streamwise  postion  of  12.8cm.  As  the  boundary 
layer  develops  along  the  upper  surface,  suction  slots  on  the  side-walls  minimize  the 
growth  of  the  side  boundary  layers  while  manual  valves  on  the  lower  surface  draw 
off  sufficient  flow  to  maintain  a  steady,  mean,  zero  pressure  gradient.  With  this 
steady  flow  as  a  mean,  unsteady  conditions  are  imposed  by  computer-controlled 
valves  which  are  also  located  on  the  lower  surface.  Since  the  valves  are  computer- 
controlled  and  can  be  placed  at  any  location  in  the  channel,  nearly  any  arbitrary, 
oscillatory,  pressure-gradient  can  be  established. 

Near  the  leading  edge  a  stationary  helium-neon  LDA  system  monitors  the  up¬ 
stream  velocity,  while  in  the  test  section  a  traversing  multiple-component  argon-ion 
LDA  system  measures  two  or  three  components  of  the  velocity  in  the  boundary 
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layer.  An  IBM  PC  performs  valve  control  while  an  AST  AT  accomplishes  the  data 
acquisition. 

The  major  subsystems  of  the  facility  in  its  present  configuration  are:  (1)  the  test 
section,  (2)  the  flow  conditioning  systems,  (3)  the  velocity  measurement  systems, 
and  (4)  the  computerized  control  and  data  acquisition  systems.  In  the  following 
sections,  each  subsystem  will  be  described  in  greater  detail. 

4.1  The  Test  Section. 

Basic  Tunnel  Construction.  The  test  section  of  the  water  channel  is  a 
stainless-steel  external  structural  frame  with  inset  transparent  panels  on  the  sides 
and  upper  surface,  and  drain  panels  on  the  lower  surface  for  either  manual  or 
computer-controlled  valves.  The  channel  was  designed  in  this  manner  to  allow  ease 
of  access  to  the  interior  of  the  channel,  maximum  flexibility  for  adapting  to  differ¬ 
ent  experimental  configurations,  maximum  optical  access  for  flow  visualization  and 
laser  measurement  systems,  and  easy  relocation  of  the  position  of  drain  valves  on 
the  lower  surface.  The  channel  is  35.6cm  wide,  12.9cm  high,  and  3.65m  long. 

Sidewall  Boundary  Layer  Control.  The  eight  side  windows  on  each  side  are 
constructed  from  lucite,  are  39.4cm  wide  by  16.5cm  high,  and  can  be  easily  removed 
and  replaced  to  allow  access  to  all  interior  parts  of  the  channel.  The  stainless  steel 
frames  surrounding  each  window  interlock  to  form  a  side  suction  slot  0.25mm  wide 
every  46cm  along  each  side  of  the  channel.  A  needle  valve  located  at  the  base  of 
each  slot  regulates  the  suction  flow  rate.  The  adjustment  of  each  valve  is  based  on 
a  combination  of  LDA  measurements,  flow  visualization,  and  the  measurment  of  its 
discharge  rate  relative  to  the  other  valves. 

The  upper  inset  panels  are  functionally  divided  into  three  sections  as  (l)  a  1.37m 
laminar  test  section  with  elliptical  leading  edge,  (2)  a  1.37m  turbulent  test  section 
with  two  smaller  drop-in  test  surfaces,  and  (3)  a  0.91m  exit  section  with  a  deflector 
plate  to  guide  the  exiting  flow. 
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Laminar  Test  Section.  The  leading  edge  of  the  laminar  test  section  is  a 
stainless-steel  semi-ellipse  with  a  16:1  ratio  of  major  to  minor  axes.  It  was  ma¬ 
chined  on  a  Matsuura  N-C  mill  and  then  hand  polished.  Flow  approaching  above 
the  stagnation  streamline  is  bypassed  to  the  sump  tank  through  a  honeycomb  struc¬ 
ture  and  four  manual  valves.  These  valves  are  carefully  adjusted  using  a  hydrogen 
bubble-wire  in  front  of  the  leading  edge  and  by  observing  the  smooth  release  of  dye 
at  a  slot  which  is  located  14.9cm  behind  the  leading  edge.  The  slot  is  formed  at 
the  junction  between  the  stainless-steel  leading-edge  structure  and  a  sheet  of  1.9cm 
thick  float  glass.  Float  glass  was  chosen  for  the  laminar  inset  plate  for  optimum  flow 
visualization  since  mechanical  access  through  the  test  surface  was  not  required.  For 
turbulent  experiments,  this  first  1.37m  section  serves  as  the  development  section 
for  the  turbulent  boundary-layer,  by  tripping  it  at  a  streamwise  location  of  12.8cm. 
The  trip  is  a  crimped  stainless-steel  wire  .13 mm  in  diameter,  which  when  rotated 
forward  serves  as  a  hydrogen  bubble  wire  at  the  leading  edge. 

Turbulent  Test  Section.  The  turbulent  test  surface  and  the  exit  section  were 
each  fabricated  from  1.9cm  lucite  with  spanwise  stiffening  ribs  every  23cm.  The 
lucite  construction  allows  for  two  smaller  and  easier  to  remove  test  surfaces  in  the 
turbulent  test  section,  each  26.6cm  by  11.4cm.  The  vortex  generators  used  in  this 
program  were  mounted  on  one  of  these  smaller  test  surfaces  located  at  a  streamwise 
position  of  1.49m.  The  configuration  of  the  vortex  generator  test  surface  will  be 
described  in  detail  in  Section  6.1.  The  exit  section  has  a  1.3cm  float-glass  plate 
situated  at  a  15°  angle  to  the  axis  of  the  water  channel.  This  provides  for  a  smooth 
exit  of  the  flow  at  the  end  of  the  channel  while  also  allowing  optical  axis  directly 
upstream  through  the  end  of  the  channel. 

Passive  Pressure  Gradient  Control.  The  lower  surface  of  the  water  chan¬ 
nel  (the  surface  opposite  the  test  surface)  is  divided  into  24  separate  drain  panel 
locations.  At  each  location,  either  a  manual  or  a  computer-controlled  valve  is  in¬ 
stalled.  The  1.9cm  diameter  manual  valves  maintain  a  mean,  steady,  zero  pressure 
gradient  while  the  10cm  diameter  computer-controlled  valves  alter  the  streamwise 
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gradient  of  the  freestream  velocity.  All  the  drain  plates  on  the  lower  surface  are 
constructed  from  0.64cm  PVC.  The  drain  hole  pattern  over  the  manual  drain  valves 
is  two  rows  of  holes,  3mm  in  diameter  staggered  on  6mm  centers  (113  holes)  with 
a  total  area  of  8.95cm2.  In  contrast,  the  plates  over  the  computer-controlled  valves 
have  9  rows  of  6mm  diameter  holes  staggered  on  13mm  centers  (167  holes)  for  a 
total  area  of  52.9cm2.  Each  drain  plate  configuration  was  designed  for  optimum 
spanwise  uniformity  while  keeping  the  throat  area  of  the  valve  as  the  controlling 
flow  resistance. 

Active  Pressure  Gradient  Control.  The  computer-controlled  valves  regulate 
the  flow  by  means  of  two  adjacent  disks  with  identical  hole  patterns  as  one  disk 
rotates  against  the  second  stationary  disk.  As  the  rotating  disk  moves  from  closed 
to  open  position,  the  two  hole  patterns  gradually  align  with  each  other  so  the  flow 
area  varies  in  a  linear  fashion  with  the  angular  position  of  the  valve.  At  the  full  open 
postion,  the  flow  area  is  8.3cm2  for  the  low  flow  valve  and  29.9 cm2  for  the  high  flow 
valve.  Each  high  flow  valve  can  pass  about  0.8m3 /min  of  water  at  the  design  head 
of  2.7 m.  The  linear  variation  of  flow  area  yields  a  linear  variation  in  the  freestream 
velocity.  Water  exiting  the  channel  through  each  drain  plate  is  funneled  through  a 
lucite  adapter  for  a  smooth  transition  from  the  14cm  by  34cm  rectangular  exit  to 
the  10cm  diameter  valve.  After  passing  through  the  valve,  water  is  directed  back 
to  the  sump  tank  through  a  10cm  diameter  PVC  pipe  in  a  manner  that  minimizes 
the  air  entrainment,  while  providing  a  free  atmospheric  pressure  interface  at  the 
disk  assembly  of  the  valve.  The  valve  is  linked  by  sprocket  and  chain  to  a  Motion 
Technologies  Model  00-01203-002  Printed  Circuit  Motor  which  is  capable  of  driving 
the  valve  at  frequencies  up  to  10 Hz.  Henk  and  Reynolds  (1989)  describe  the  valve 
construction  in  more  detail,  and  the  valve  control  will  be  described  more  in  section 
4.4. 
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4.2  Flow  Conditioning  Systems 


The  water  entering  the  test  section  from  the  nozzle  is  carefully  conditioned  to 
minimize  freestream  turbulence,  to  maintain  temperature  within  0.1°C,  to  eliminate 
microscopic  dissolved  air  bubbles,  and  to  control  the  size  of  particles  in  the  flow. 

Flow  Straightening.  Straightening  of  the  flow  and  minimization  of  free-stream 
turbulence  begins  in  the  constant  head  tank.  A  constant  area  stainless-steel  adapter 
smoothly  adjusts  the  flow  geometry  from  round  15cm  diameter  pipe  to  rectangular. 
The  flow  then  turns  and  decelerates  through  three  honeycomb  blocks  as  the  flow 
area  increases  to  92cm  by  92cm  in  the  stilling  chamber,  as  shown  in  figure  4.1.  Each 
honeycomb  section  has  a  nominal  cell  size  of  1cm  and  is  13cm  thick.  In  the  stilling 
chamber,  the  flow  passes  through  three  22mesh  screens  spaced  10cm  apart.  The 
screen  material  is  0.1905mm  stainless  steel  wire  which  results  in  69%  open  area. 
The  screens  are  stretched  tightly  over  stainless  steel  frames,  and  can  be  removed 
from  the  stilling  chamber  periodically  for  cleaning.  Finally,  the  flow  passes  through 
a  20:1  contraction  into  the  entrance  of  the  water  channel.  The  resulting  freestream 
turbulence  level  is  approximately  0.5%. 

Temperature  Control.  During  operation  of  the  water  channel,  the  water  tem¬ 
perature  is  stabilized  and  maintained  at  16. 7°  C  ±0.1°C.  A  2.24 kW  pump  operates 
in  parallel  with  the  main  5.6 kW  pump,  directing  flow  from  the  sump  tank  to  a 
water/ethylene-glycol  heat-exchanger  which  is  chilled  by  a  35 kW  Trane  refrigera¬ 
tion  system.  While  the  water-flow  through  the  heat-exchanger  runs  continuously, 
the  temperature  control  system  enables  or  disables  the  refrigeration  system  in  re¬ 
sponse  to  a  thermistor  located  in  a  region  of  the  sump  tank  where  very  high  mixing 
occurs. 

Deaeration.  Air  is  removea  during  system  start-up  using  a  back-flow  purge 
system,  and  then  during  operation  dissolved  air  is  continuously  removed  with  a 
deaeration  system.  The  purge  system  is  necessary  during  start-up  to  ensure  that 
no  trapped  air  bubbles  in  any  of  the  suction  holes  or  slots  alter  the  flow  rate  through 
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each  valve.  Once  the  system  is  running,  microscopic  air  bubbles  must  be  removed 
These  bubbles  come  from  air  which  was  originally  dissolved  in  the  tap  water  and  air 
that  is  entrained  at  the  air/water  interface  in  the  sump  tank.  For  this  purpose,  a 
deaeration  system  runs  continuously.  A  portion  of  the  water  leaving  the  main  pump 
is  bypassed  through  a  venturi  which  drops  the  pressure  sufficiently  for  boiling  to 
occur  at  room  temperature.  The  water  vapor  tends  to  nucleate  on  the  microscopic 
air  bubbles  which  then  combine  and  grow  in  size  as  they  convect  along  a  horizontal 
pipe  2m  in  length.  A  mechanical  vacuum  pump  then  draws  off  the  air  and  water 
vapor  in  a  settling  chamber,  while  the  water  is  returned  to  the  main  pump  inlet 
through  a  smaller  l.lArW  pump. 

Particulate  Control.  The  particle  content  in  the  flow  is  controlled  both  while 
the  system  is  filling  and  during  operation.  Tap  water  entering  the  system  passes 
through  a  125mtcron  and  a  bmicron  filter  in  series  so  that  the  only  particles  entering 
the  water  channel  are  5 microns  or  smaller.  This  is  the  optimum  particle  size  for 
light  scattering  with  LDA  systems  in  water.  To  prevent  the  build-up  of  any  larger 
particles  in  the  system  while  it  is  operating,  all  the  tanks  in  the  system  are  kept 
carefully  covered,  and  a  set  of  three  parallel  filters  in  the  chiller  loop  removes  any 
impurities  that  do  get  into  the  system. 

4.3  Velocity  Measurement  Systems 

Two  LDA  systems  are  used  for  velocity  measurement.  This  enables  a  single¬ 
component  system  to  monitor  the  steady  upstream  conditions  while  a  multiple- 
component  system  measures  velocities  in  the  boundary-layer  downstream.  In  ad¬ 
dition,  the  two  systems  can  be  located  side  by  side  in  a  steady,  freestream  flow  to 
verify  the  calibration  of  each  system. 

Single  Component  Upstream  Monitoring.  The  single-component  system 
uses  a  5 mW  helium-neon  laser  with  Dantec  transmitting  and  receiving  optics.  A 
custom  fixture  secures  the  laser  and  Dantec  optics  on  a  platform  that  rests  on  the 
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water  channel.  Two  turning  mirrors  mounted  on  the  fixture  direct  the  transmitted 
beams  down  and  across  the  channel,  and  two  additional  turning  mirrors  on  the 
fixture  direct  the  scattered  light  back  up  to  the  receiving  optics.  This  fixture  allows 
easy  placement  and  removal  of  the  system  anywhere  on  the  channel,  as  necessary. 
The  LDA  signal  is  processed  by  a  Dantec  model  55K24  tracker,  which  gives  a 
direct  digital  reading  of  the  mean  velocity  and  sends  an  analog  signal  to  the  data 
acquisition  system. 

Multicomponent  Measurements.  The  multiple-component  LDA  system  uses 
a  4 W  Lexel  model  95  argon-ion  laser  with  TSI  transmitting  and  receiving  optics. 
While  many  different  optical  arrangements  were  tried  in  the  preliminary  phases 
of  this  program,  two  arrangements  evolved  which  were  used  for  final  data  acquisi¬ 
tion,  depending  on  whether  two  or  three  velocity  components  were  to  be  measured. 
Schematics  of  the  two-component  system  are  shown  in  figures  4.2  through  4.4,  and 
schematics  of  the  three-component  system  are  shown  in  figures  4.5  and  4.6.  In  ad¬ 
dition,  tables  4.1  and  4.2  summarize  the  key  optical  components  in  the  two  LDA 
systems.  In  both  configurations,  the  optics  are  located  on  a  breadboard  which  al¬ 
lows  automated  simultaneous  traverses  of  the  transmitting  and  receiving  optics  in 
the  direction  normal  to  the  test  surface.  In  addition,  linear  bearings  allow  manual 
motion  in  the  spanwise  direction,  and  trolley  wheels  along  two  I-beams  allow  man¬ 
ual  motion  in  the  streamwise  direction.  Also,  in  both  the  three-component  and  the 
four-component  systems,  the  standard  TSI  dichroic  color  separator  was  modified 
by  removing  two  beam  steering  prisms.  This  enabled  a  single  beam  to  enter  the 
separator  module  on-axis  and  allowed  the  green  beam  to  exit  on-axis,  while  the  blue 
beam  was  displaced  25 mm. 

Two-Component  LDA.  The  two-component  LDA  system  uses  a  modified  TSI 
dichroic  color  separator  (see  figure  4.4)  to  create  a  four-beam,  forward  scatter  sys¬ 
tem.  Two  green  beams  (514. 5nm)  create  a  measuring  volume  along  the  z-axis  for 
measuring  the  u-velocity  component,  while  two  blue  beams  (488nm)  create  a  mea¬ 
suring  volume  for  measuring  the  u-velocity  component.  One  blue  beam  is  offset 
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back  to  the  center  to  allow  measurements  to  be  made  very  close  to  the  test  surface. 
Figure  4.7a  shows  the  configuration  of  the  four  beams.  This  particular  configuration 
was  chosen  because  the  green  beams  create  a  measuring  volume  that  contains  no 
contamination  from  v  or  w,  which  in  some  of  the  separating  flows  under  investigation 
could  be  a  significant  percentage  of  u.  The  orientation  of  the  v  measuring  volume 
accepts  about  3°  of  w  contamination  in  order  to  make  measurements  very  near  the 
test  surface.  One  of  the  green  beams  is  Bragg-shifted  at  38 MHz  and  one  of  the  blue 
beams  is  shifted  at  40 MHz.  This  allows  measurement  of  reversing  flows  in  both 
velocity  components,  and  gives  frequency  separation  in  addition  to  color  separation 
for  processing  the  scattered  light.  The  scattered  light  passes  through  a  50 micron 
pinhole  in  the  field  stop  system  then  is  dichroically  separated  and  processed  by 
separate  photo-multiplier  tubes  (PMT’s).  Each  PMT  signal  is  then  downmixed 
to  200 kHz  by  TSI  model  9186 A  mixers,  and  routed  to  TSI  model  1090  trackers 
through  a  Krone-Hite  model  3202  low-pass  filter  set  at  800 kHz.  The  filter  has  an 
attenuation  of  3 dB  at  the  indicated  frequency  with  a  roll-off  of  80 dB/ decade. 

Three-Component  LDA.  A  new  three-component  LDA  system  collected  all 
the  primary  data  in  this  project.  It  also  uses  a  modified  TSI  dichroic  color  separator 
(see  figure  4.4)  to  yield  three  green  beams  in  forward-scatter  and  two  blue  beams 
in  side  scatter.  Figure  4.7b  summarizes  the  configuration  of  the  five  beams. 

Forward  Scatter  Measurement  of  u  and  v.  The  three-beam  forward-scatter  system 
meaures  u  and  v,  each  with  about  3°  of  w  contamination.  The  forward  shooting  u 
beam  is  shifted  at  40MHz  and  the  upward  shooting  v  beam  is  shifted  at  38 MHz 
so  that  reversing  flows  in  both  u  and  v  can  be  measured  and  so  that  frequency 
separation  can  be  used  as  the  primary  means  to  separate  the  u  and  v  signals  during 
processing  of  the  scattered  light.  The  doppler  bursts  from  both  green  measuring 
volumes  are  detected  by  a  single  PMT  after  passing  through  the  dichroic  color 
separator  and  a  narrow-band  color  filter.  The  PMT  signal  is  split  to  enter  two  TSI 
9186A  mixers.  Since  one  signal  is  down-mixed  at  200 kHz  relative  to  40 MHz  and 
one  relative  to  38 MHz,  frequency  separation  occurs  when  each  mixer  output  passes 
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through  the  effective  low-pass  filter  of  each  tracker.  The  u  and  v  signals  coming 
out  of  the  low-pass  filters  are  then  processed  by  TSI  trackers  in  the  same  manner 
as  the  two-component  system  above. 

Side  Scatter  Measurement  of  w.  The  two  transmitted  blue  beams  of  the  side- 
scatter  system  enter  through  the  transparent  upper  test-surface  of  the  water  channel 
and  are  oriented  so  as  to  measure  very  nearly  pure  w.  During  the  alignment  proce¬ 
dure,  the  measurement  volume  is  oriented  at  a  45°  angle  to  the  x  and  z  axes  rather 
than  measuring  pure  w,  because  it  is  impossible  to  observe  the  beam  crossing  in  the 
receiving  optics  when  viewed  “straight-on”.  After  the  optics  are  sufficiently  aligned 
such  that  a  doppler  signal  can  be  detected  by  a  tracker  or  counter,  the  optics  are 
carefully  rotated  back  to  a  position  where  pure  w  is  measured  for  final  optimization 
of  the  signal.  To  ensure  that  the  path  length  in  air  and  water  does  not  change  as 
the  system  traverses  in  the  direction  normal  to  the  surface,  a  water-tight  air-tube 
traverses  up  and  down  in  a  miniature  tank  of  water,  which  is  also  maintained  at 
constant  temperature.  The  side-scattered  blue  light  passes  through  the  same  field- 
stop  system  as  the  green  forward-scattered  light.  It  is  then  separated  from  the  green 
light  by  the  dichroic  color  separator,  and  then  passes  through  a  narrow  band-pass 
color  filter  into  a  PMT.  Selection  of  the  appropriate  pinhole  in  the  field-stop  system 
is  very  critical  to  good  turbulence  measurements,  and  this  selection  becomes  much 
more  complex  for  a  system  using  a  combination  of  forward  scatter  and  side  scatter. 
This  critical  pinhole  issue  will  be  discussed  at  length  in  section  5.3 

Side  Scatter  Signal  Processing.  Since  the  count  rate  is  much  lower  in  side  scatter, 
the  PMT  output  is  processed  by  a  TSI  model  1980  Counter  after  being  downmixed 
to  200 kHz.  The  low-pass  filter  in  the  counter  has  an  attenuation  of  lOOdB/decade. 
In  contrast  to  trackers  which  need  nearly  continuous  doppler  bursts  to  track  prop¬ 
erly,  a  counter  works  well  with  lower  count  rates.  The  analog  output  from  the  two 
trackers  and  the  counter  then  feed  into  the  analog-to-digital  input  of  the  data  ac¬ 
quisition  system,  where  the  appropriate  mathematical  operations  convert  voltages 
into  velocities. 
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To  achieve  reasonable  count  rates  in  side  scatter,  the  flow  had  to  be  seeded.  While 
several  different  seeding  materials  were  tried,  the  seed  that  was  least  expensive  yet 
gave  good  results  was  Titanium  Oxide  (Ti02).  Initially,  about  4cm3  of  seed  was 
added  to  the  5.7 m3  of  water  in  the  system.  While  the  data  validation  rate  for  the 
forward  scattered  signals  varied  very  little  with  the  addition  of  seed,  the  count  rate 
on  the  side  scattered  signal  increased  from  about  50cps  to  600cps.  Because  the 
filters  in  the  system  tended  to  slowly  remove  seed,  occasionally  1cm3  of  seed  was 
added  to  maintain  the  count  rate  in  the  freestream  at  about  600cps. 

One  of  the  beams  was  shifted  at  40MHz  to  allow  measurement  of  reversing 
flows.  Since  both  the  green  u  signal  and  the  blue  w  signal  were  shifted  at  40 MHz, 
and  since  the  dichroic  color  separation  in  the  receiving  optics  is  not  100  percent 
efficient,  each  signal  was  checked  for  contamination  by  the  other  signal.  When 
the  counter  was  adjusted  to  count  bursts  from  the  w  signal  at  a  rate  of  400  to 
800cps,  the  w  beams  were  blocked  while  the  u  beams  were  not.  The  counter  was 
not  able  to  validate  any  stray  bursts  from  the  u  signal.  A  similar  test  confirmed  that 
the  trackers  were  unable  to  track  stray  w  signals  when  the  u  beams  were  blocked. 
Therefore,  cross-contamination  was  not  a  problem. 

4.4  Computerized  Control  and  Data  Acquisition 

Both  the  computer  valve-control  system  and  the  data  acquisition  system  operate 
on  either  IBM  or  IBM-compatible  personal  computers.  Both  of  these  systems  run 
with  menu-driven  software  developed  at  Stanford,  and  the  valve-control  system  uses 
control  boards  which  were  also  developed  at  Stanford.  The  details  of  these  systems 
are  described  by  Carlson  and  Reynolds  (1989). 

Valve  Control.  The  valve-control  system  operates  on  an  IBM  PC,  with  custom- 
designed  control  boards,  and  a  Sorensen  DCR  40- 125 A  Power  Supply.  The  system 
allows  the  user,  through  a  menu-driven  system,  to  select  the  waveform  and  fre¬ 
quency  for  a  set  of  up  to  eight  valves.  Waveforms  presently  programmed  are  a  sine 
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wave,  square  wave,  and  triangular  wave,  but  since  the  program  generates  its  own 
look-up  table,  virtually  any  waveform  could  be  generated.  Once  the  user  selects  the 
waveform  and  frequency  for  the  set  of  valves,  the  amplitude  and  phase  of  each  of  the 
eight  valves  can  be  set  independently.  Typically,  upstream  and  downstream  valves 
are  operated  out  of  phase  but  with  balanced  amplitudes  so  the  flow  downstream  of 
the  upstream  valve  oscillates  while  the  flow  upstream  remains  steady.  The  valve 
controller  sends  two  pulses  to  the  data-acquisition  system.  The  first  pulse  is  contin¬ 
uous  for  triggering  phase-sampling  at  specified  phases  of  the  cycle,  and  the  second 
pulse  identifies  the  zero-phase  point  of  each  cycle. 

Data  Acquisition.  The  data  acquisition  system  operates  on  a  10 MHz  AST 
AT  with  an  80 MByte  hard  disk.  Analog-to-digital  processing  of  the  LDA  signals 
is  accomplished  by  a  Tecmar  model  TM-40-PGH  board  and  the  traverse  is  run 
through  the  computer  digital-out.  The  menu-driven  software  allows  the  operator 
to  specify  all  the  parameters  of  sampling  through  a  set  of  input  control  files.  The 
system  is  presently  capable  of  simultaneously  collecting  velocity  data  from  three 
LDA  systems.  From  the  three  input  signals,  the  three  velocity  components  are  cal¬ 
culated,  then  individual  phase  ensembles  are  collected  for  the  appropriate  statistical 
quantities.  When  the  LDA  is  operated  with  two  components,  two  mean  velocities 
and  three  Reynolds  stress  terms  are  collected  at  each  phase  point.  When  the  LDA 
is  operated  with  three  velocity  components,  three  mean  velocities,  six  Reynolds 
stress  terms,  and  ten  triple-product  terms  are  collected  for  each  phase  point.  All 
three  channels  are  sampled  on  command  within  a  period  of  70  microseconds.  The 
equivalent  correlation  length  associated  with  this  time  delay  is  y+  <  0.5  (based 
on  Uqq  =  0.5 m/sec  and  ur  =  0.021m/sec).  This  is  much  smaller  than  the  effec¬ 
tive  measurement  volume  (y+  ~  10),  and  thus  does  not  degrade  the  data.  The 
data-acquisition  system  can  control  the  traverse  system  so  that  an  entire  set  of 
profiles  is  taken  of  each  of  these  quantities,  at  each  phase,  at  each  y- location  in  the 
boundary-layer. 


41 


Unsteady  Flow  Facility 


Figure  4.1  Schematic  of  unsteady  flow  facility 


Figure  4.2  Two-component  LDA  system 
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Figure  4.3  Two-color  transmitting  and  receiving  optics 


Standard  Dichroic  Color  Separator  (TSI  9112) 


- blue-green 

- green 

. blue 


Modified  Dichroic  Color  Separator  (prisms  removed) 


Figure  4.4  Dichroic  color  separator 
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Figure  4.6  Three-component  transmitting  optics 
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- green 

- 2  blue 


Figure  4.7b  Beam  configuration  for  three-component  system 
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Table  4.1  Summary  of  Two-Component  LDA  Components 


Component 

Description 

Part  Number 

TM1 

lm  dia.  mirror,  45° 

TSI  9109 

TM2 

4 in  dia  mirror 

NRC  630A-4 

TM3 

6m  dia.  mirror 

NRC  630A-6 

TM4 

4m  dia.  mirror 

NRC  630A-4 

LI 

450mm  lens/178mm  dia. 

TSI  9169-450 

L2 

500mm  lens/83mm  dia. 

TSI  9167-500 

Transmitting  Optics 

see  figure  4.3 

Beam  Expander 

3.75X 

TSI  9189 

Receiving  Optics 

see  figure  4.3 

Table  4.2  Summary  of  Additional  Three- Component  LDA  Components 


Component 

TM5 

TM6 

TM7 

L3 

Color  Separator 
U-V  Transmitting  Optics 
W  Transmitting  Optics 
Receiving  Optics 


Description 
lin  dia.  mirror,  45° 

4m  dia.  mirror 
4m  dia.  mirror 
200mm  lens /76mm  dia. 
see  figure  4.4 
see  figure  4.0 
see  figure  4.6 
see  figure  4.3 


Part  Number 
TSI  9109 
NRC  600A-4R 
NRC  600A-4R 
NRC  KPX229AR.14 
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5.  Facility  Validation 


Facility  validation  involves  verification  that  the  quality  of  the  flow  meets  the 
basic  requirements  of  the  experiment  and  that  the  measurement  techniques  and  in¬ 
strumentation  are  yielding  correct,  consistent  information.  The  key  issues  regarding 
quality  of  the  flow  are  repeatability  of  the  steady,  mean  flow;  streamwise  unifor¬ 
mity  of  the  velocity  for  zero  pressure  gradient  work;  low  freestream  turbulence; 
spanwise  uniformity  in  the  flow;  and  the  absence  of  any  vortex-type  flow  structures 
that  would  detract  from  the  two-dimensional  nature  of  the  boundary  layer.  Issues 
regarding  the  velocity  measurements  are  the  fundamental  accuracy  of  the  absolute 
magnitude  of  the  mean  velocity,  the  correctness  of  the  RMS  fluctuating  quantities, 
and  the  proper  spatial  and  temporal  resolution  when  two  fluctuating  quantities  are 
correlated.  To  verify  the  flow  quality  one  would  like  to  have  a  well-proven  veloc¬ 
ity  measurement  system,  and  conversely  to  verify  the  LDA  system  one  would  like 
to  have  a  well  known,  well  behaved  flow.  In  this  program,  where  both  the  water 
channel  and  the  LDA  optical  configuration  were  new,  a  systematic  validation  pro¬ 
gram  was  required  to  verify  both  the  quality  of  the  flow  and  the  credibility  of  LDA 
measurements. 

5.1  Characterization  of  the  Freestream  Velocity 

For  initial  validation  of  LDA  measurements,  the  water  channel  simulated  a  nom¬ 
inally  steady,  two-dimensional  flow  over  a  flat  plate.  First,  the  flow  bypass  around 
the  leading  edge  was  adjusted  by  flow  visualization  for  optimum  stagnation  con¬ 
ditions  at  the  leading  edge.  Then  the  opposite  wall  suction  was  adjusted  for  a 
nominally  uniform  freestream  velocity,  using  the  LDA  as  an  indicator  of  uniform 
relative  velocity.  With  a  relatively  good  flow  along  a  single  streamline,  then,  the 
measurement  volumes  of  the  TSI  system  and  the  Dantec  system  were  moved  within 
several  centimeters  of  each  other  for  comparison  of  absolute  velocity  measurements. 
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Although  velocity  measurement  with  an  LDA  system  is  fundamental  in  nature, 
(that  is,  the  time  to  traverse  between  fringes  in  an  interference  pattern  is  measured 
rather  than  inferred  from  a  secondary  measurement  such  as  heat  transfer)  careless 
determination  of  crossing  angles  can  lead  to  errors  as  large  as  10%  in  the  abso¬ 
lute  measurement  of  velocities.  By  projecting  the  transmitting  beams  of  the  LDA 
through  the  channel  wall  and  across  the  width  of  the  lab,  lengths  from  20cm  to  1.7m 
could  be  used  to  determine  the  crossing  angles,  which  typically  are  between  6°  and 
12°,  to  a  precision  of  about  1%.  After  this  careful  calibration,  the  measurements 
from  the  TSI  and  Dantec  systems  agreed  with  each  other  to  within  about  0.5%. 

RMS  Measurements.  Since  two  very  different  LDA  systems  were  able  to  inde¬ 
pendently  verify  the  absolute  magnitude  of  velocity  measured  by  the  other  system, 
confidence  in  the  mean  value  of  LDA  measurements  is  very  high.  The  remaining 
question  about  fundamental  LDA  measurements  was  the  RMS  value  of  fluctuations 
in  the  u-component  of  velocity,  u'2,  and  what  it  implied  about  the  freestream  turbu¬ 
lence  intensity,  \fvP-  jU .  While  some  of  the  data  presented  in  subsequent  chapters 
indicate  that  the  apparent  turbulence  intensity  is  about  1%  to  1.5%,  studies  of 
various  LDA  configurations  indicated  that  the  actual  turbulence  intensity  of  the 
water  channel  is  lower,  about  0.5%  to  0.7%,  The  difference  appears  to  come  from 
electronic  noise  introduced  with  the  use  of  Bragg  cells  and  the  down-mixing  cir¬ 
cuitry.  This  was  demonstrated  by  setting  up  the  TSI  optics  in  a  simple  two-beam 
system  with  no  color  separation  and  no  Bragg  cells.  With  careful  adjustment  of  the 
collimator,  the  freestream  turbulence  level  was  measured  to  be  0.5%  to  0.7%  at  a 
freestream  velocity  of  0.7 m/ sec.  Without  changing  any  flow  conditions,  the  optics 
were  then  modified  to  a  simple  two-beam  system  with  a  Bragg  cell.  The  appar¬ 
ent  freestream  turbulence  intensity  then  changed  to  between  1%  to  1.5%,  and  no 
optimizing  of  the  system  seemed  to  improve  that.  Since  the  turbulence  intensities 
inside  the  boundary-layer  are  much  greater  than  this,  this  limitation  in  the  LDA 
system  was  deemed  to  be  acceptable. 
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Spanwise  Uniformity  of  the  Flow.  Early  in  the  testing  of  the  new  water 
channel,  studies  of  the  spanwise  uniformity  revealed  several  regions  in  the  —z  por¬ 
tion  of  the  channel  where  the  boundary  layer  was  thicker  than  over  the  rest  of  the 
span.  The  problem  was  traced  to  leakage  from  the  leading  edge  bypass  back  into 
the  dye  slot  near  the  leading  edge.  When  this  was  corrected,  the  mean  velocity 
showed  no  major  spanwise  variations.  These  results  are  shown  in  figure  5.1.  The 
spanwise  uniformity  and  symmetry  of  the  flow  behind  the  vortex  generators  under 
unsteady  flow  conditions  will  be  presented  in  chapter  6,  both  for  mean  velocities 
and  RMS  measurements. 

5.2  Mean  Velocity  Measurements  in  a  Steady,  Zero  Pressure-Gradient 
Turbulent  Flow 

With  the  fundamental  LDA  issues  resolved,  the  next  logical  step  to  validate 
was  the  behavior  of  the  turbulent  boundary-layer  in  a  zero  pressure-gradient  flow, 
and  the  ability  of  the  LDA  system  to  measure  the  mean  two-dimensional  velocity 
components.  Figure  5.2  shows  the  streamwise  velocity  history  of  the  flow  in  the 
configuration  used  for  the  following  validation  data.  The  opposite  wall  suction  was 
adjusted  for  a  constant  freestream  velocity  at  about  0.5 mj sec.  Figure  5.2  shows 
that  the  variation  from  this  mean  value  is  less  than  0.5%  throughout  the  entire  test 
section. 

During  data  acquisition,  a  steady  zero  pressure-gradient  profile  is  taken  at  each 
spanwise  position  in  order  to  establish  the  apparent  position  of  the  LDA  measuring 
volume  and  to  determine  the  integral  boundary-layer  parameters.  In  figure  5.3,  one 
such  profile  is  shown  in  wall  coordinates.  The  outer  region  is  fit  to  Coles’  Law, 

“+  =  0^1^+>  +  5-0  +  0iIM,(f)'  <51> 

where 

*(S) =“■’(!!)  <*•*> 
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while  the  linear  sublayer  is  fit  to 


u+  =  j/+.  (5.3) 

In  each  profile  that  is  analyzed,  the  apparent  position  of  the  measuring  volume  is 
first  inferred  from  the  best  fit  of  a  line  to  the  inner  points.  Then,  the  outer  region 
is  fit  by  linear  regression  to  equation  (5.1),  which  determines  the  values  of  6  and 
uT.  It  is  apparent  from  the  figure  that  the  data  fits  the  equations  very  well.  In 
addition,  the  value  of  n  from  the  linear  regression  for  the  case  shown  was  0.489. 
This  value  is  relatively  close  to  the  value  of  0.55  suggested  by  Coles  (1956)  for  a  zero 
pressure  gradient  turbulent  boundary  layer.  Finally,  figure  5.3  also  shows  very  good 
agreement  with  predictions  by  Spalart  (1986).  Therefore,  the  mean  values  of  the 
streamwise  velocity  component,  u,  appear  to  be  very  consistent  with  all  predictions. 

The  u-component  of  the  mean  velocity  is  much  more  difficult  to  measure,  since 
a  very  small  u-contamination  (less  than  0.5°)  or  a  slight  error  in  the  value  of  the 
downshift  frequency  as  read  by  the  data  acquisition  system  can  completely  obliterate 
the  correct  mean  value  of  v.  Figure  5.4  shows  v  profiles  taken  at  several  spanwise 
locations.  These  profiles  are  in  general  agreement  with  the  profiles  reported  by 
Brereton  and  Reynolds  (1987). 

The  tu-component  of  the  mean  velocity  is  also  difficult  to  validate,  since  it  is 
not  present  in  a  two-dimensional  flow.  However,  the  LDA  side-scatter  measuring 
volume  can  be  oriented  to  measure  u  cos  0  +  w  sin  0,  where  0  ~  45°.  By  orienting 
the  w  side-scatter  optics  at  45°  in  a  two-dimensional  flow,  the  u  velocity  component 
was  measured  simultaneously  with  the  normal  forward  scatter  u  optics  and  with 
the  45°  side-scatter  w  optics.  The  two  values  agreed  within  the  1%  uncertainty  of 
the  velocity  calibrations. 

5.3  Measurement  of  Second  Order  Statistics 

The  most  difficult  measurements  to  make,  and  hence  the  most  difficult  to  val¬ 
idate,  are  the  fluctuating  quantities  and  correlations  in  a  turbulent  flow.  After 
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setting  up  an  LDA  system,  one  can  easily  obtain  measurements  that  have  a  mea- 
sureable  fluctuation,  but  getting  meaningful  turbulent  quantities  takes  some  care. 
The  key  indicators  of  good  turbulent  measurements  are  the  frequency  distributions 
or  power  spectra,  known  values  of  correlations  between  u'  and  v'  in  shear  layers, 
and  published  results  for  standard  experiments  such  as  steady  flow  over  a  flat  plate. 

Measurement  of  Normal  Reynolds  Stresses.  The  best  “first  test”  of  the 
goodness  of  fluctuating  quantities  is  to  look  at  the  spectral  content  of  the  apparent 
turbulent  kinetic  energy.  Figure  5.5a  shows  the  energy  spectrum  of  u'2  and  v'2 
under  steady  flow  at  a  position,  y/S  ^  0.4,  and  figure  5.5b  shows  the  same  for  w ,2. 
The  u12  spectrum  shows  a  roll-off  at  about  1  Hz.  Consistent  with  Hinze’s  (1981) 
discussion  of  homogeneous  shear-flow  turbulence,  there  is  a  region  where  the  slope  is 
—  1  and  a  region  where  it  is  —5/3.  The  v12  spectrum  shows  a  roll-off  at  a  somewhat 
higher  frequency  and  thus,  within  the  frequency  limitations  of  the  FFT  taken,  only 
the  region  with  slope  of  —1  is  seen. 

Both  the  u12  and  the  v'2  spectra  and  the  relationship  between  the  two  are  virtually 
identical  to  those  reported  by  Brereton  and  Reynolds  (1987).  In  addition,  the 
spectra  for  v'2  and  w'2  are  very  similar  in  magnitude  and  shape,  and  this  is  consistent 
with  the  calculations  of  Kim,  Moin,  and  Moser  (1987)  in  the  outer  portion  of  the 
boundary  layer.  Finally,  for  an  unsteady  driven  flow,  it  is  important  to  compare 
the  spectra  of  the  driven  turbulent  flow  with  that  of  the  steady  flow.  Figure  5.6 
shows  the  spectra  of  u'2  while  the  flow  is  driven  at  0.5 Hz.  While  the  flow  shows  a 
clear,  sharp  peak  at  0.5 Hz,  the  roll-off  location  and  characteristic  slopes  have  not 
changed  from  the  steady  flow  spectra. 

One  final  evaluation  of  the  normal  stresses  is  comparison  .vith  published  results 
for  turbulent  steady  flow  over  a  flat  plate.  These  profiles  are  shown  in  figure  5.7a  and 
5.7b  along  with  the  flat  plate  simulation  results  of  Spalart  (1986)  for  a  momentum 
thickness  Reynolds  number  of  1410.  Except  for  some  wall  interference  with  values 
of  w' /ur,  the  two  are  in  very  good  agreement.  Note  that  figure  5.7a  shows  good 
near  wall  agreement.  For  comparison  in  the  outer  region  of  the  boundary  layer, 
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the  results  are  replotted  in  figure  5.7b  against  y/S,  to  accomodate  the  difference 
in  Reynolds  numbers  (Re#  =  1840  for  this  program).  The  values  shown  for  u' /ur 
and  v1  juT  are  also  in  very  good  agreement  with  the  experimental  results  by  Laufer 
(1950)  and  Brereton  and  Reynolds  (1987). 

Shear  Stress  Measurements.  Measurement  of  the  shear  stresses  adds  another 
level  of  complexity.  The  one-point  correlations  will  only  correlate  well  if  the  corre¬ 
sponding  measuring  volumes  are  indeed  located  at  nearly  the  same  point  in  space. 
Blackwelder  (1983)  noted  that  the  active  length  of  a  hot-wire  probe  needed  to  be  on 
the  order  of  y+  <  20  to  obtain  good  near-wall  correlations.  For  a  properly  aligned 
LDA  system,  the  diameter  of  the  measuring  volume  is  normally  much  less  than 
this,  but  the  length  of  the  measuring  volume  is  much  greater,  so  a  pinhole  is  used 
to  limit  this  length.  While  an  off-axis  forward-scatter  system  places  some  rather  se¬ 
vere  constraints  on  the  pinhole  selection,  a  combination  of  side-scatter  and  forward 
scatter  places  additional  conflicting  constraints  on  the  pinhole  requirements. 

Figure  5.8  demonstrates  why  there  are  conflicting  requirements  on  selection  of 
the  optimum  pinhole.  The  critical  issue  is  the  size  of  the  image  of  the  limiting 
aperature  as  it  is  projected  onto  the  measuring  volume.  Figure  5.8  shows  that  the 
pinhole  image  is  extended  by  the  inverse  of  the  cosine  of  15°  or  a  factor  of  1.04  as  it  is 
projected  on  the  side-scatter  measuring  volume,  but  it  is  extended  by  the  inverse  of 
the  sine  of  15°  or  a  factor  of  3.85  as  it  is  projected  on  the  forward-scatter  measuring 
volume.  Therefore,  to  obtain  the  same  effective  size  of  measuring  volumes,  different 
pinhole  sizes  are  required.  The  following  are  the  key  factors  and  decisions  that  went 
into  the  selection  of  pinholes: 

(1)  There  are  three  pinhole  locations,  the  blue  PMT,  the  green  PMT,  and  the 
field-stop  system.  Also,  the  imaging  optics  of  the  PMTs  magnify  the  pinhole  image 
by  a  factor  of  2.5,  while  the  field  stop  system  magnifies  the  pinhole  image  by  a 
factor  of  5. 
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(2)  Correlations  between  u  and  v  are  the  most  difficult  because  the  maximum 
uV  values  occur  in  small  eddies  very  close  to  the  surface.  This  task  is  simplified 
somewhat  in  the  three-component  system  because  the  same  PMT  is  used  for  both 
u  and  v.  Therefore,  a  50  micron  pinhole  placed  in  the  green  PMT,  with  a  magni¬ 
fication  of  2.5  and  a  projection  factor  of  3.85,  has  an  image  length  of  481mtcrons. 
For  the  run  conditions  of  this  experiment,  the  measurement  volume  length  in  wall 
coordinates  is  therefore  £+  =  9,  which  should  be  satisfactory  for  good  ( u'v ')  cor¬ 
relations.  Observations  of  the  doppler  bursts  on  a  dual-trace  storage  oscilloscope 
confirmed  that  the  u  and  v  trackers  were  both  tracking  the  same  bursts. 

(3)  In  a  similar  fashion,  a  200micron  pinhole  placed  in  the  blue  PMT,  with  a 
magnification  of  2.5  and  a  projection  factor  of  1.04,  has  an  image  of  520 microns 
which  in  wall  coordinates  is  £+  =  10. 

(4)  The  field-stop  system  will  control  the  spatial  proximity  of  the  tv  measuring 
volume  to  the  u  and  v  measuring  volumes.  Due  to  the  difference  in  magnification 
factors,  a  lOOmtcron  pinhole  in  the  field-stop  system  also  results  in  a  520 micron 
image  on  the  side-scatter  optics.  The  best  compromise  possible  was  to  let  the 
smaller  image  of  the  green  PMT  50micron  pinhole  look  through  the  center  of  this 
larger  pinhole.  In  practice,  since  traversing  requires  near  perfect  alignment  of  two 
sets  of  transmitting  optics  and  one  set  of  receiving  optics,  the  100 micron  pinhole 
in  the  field-stop  system  would  have  required  correction  of  the  alignment  during 
the  traverse.  Since  the  unsteady  profiles  each  took  26  hours  of  continuous  data 
acquisition,  this  was  not  acceptable,  and  a  200micron  pinhole  had  to  be  used  in  the 
field-stop  system. 

(5)  The  net  result  of  the  pinhole  tradeoffs  is  that  the  u  and  v  measuring  volumes 
are  co-located  with  an  effective  length  of  £+  =  9.  The  tv  measuring  volume  has 
an  effective  length  of  £+  =  10,  but  at  worse  case  the  u  and  v  measuring  volumes 
could  be  offset  from  the  tv  measuring  volume  by  as  much  as  £+  =  ±20.  This  trade¬ 
off  was  deemed  to  be  acceptable  since  the  correlations  of  u  and  v  with  w  generally 
would  result  from  larger  scale  structures  in  the  outer  boundary  layer  induced  by  the 
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vortex  generators.  Also,  this  measuring  volume  is  still  small  compared  to  the  X-wire 
probes  used  by  Pauley  and  Eaton  (1988)  who  still  obtained  very  good  correlations. 

A  key  indicator  of  quality  shear  stress  measurements  as  reported  by  both  Bar- 
low  and  Johnston  (1985)  and  Brereton  and  Reynolds  (1987),  is  the  behavior  of 
correlations  between  u'  and  v'  very  near  the  test  surface.  In  figure  5.9,  the  corre¬ 
lation  coefficient,  —u'v'/y/vfiyfr/i,  is  plotted  along  with  the  results  from  Spalart 
(1986).  In  figure  5.10,  v' /u'  is  plotted  in  a  similar  manner,  for  values  of  y+  near  the 
wall.  Barlow  observed  that,  with  properly  aligned  optics,  the  conclation  coefficient 
reached  a  value  of  0.4  by  y+  values  of  10  to  20.  In  addition,  he  observed  that  the 
ratio  of  v' /u'  remained  linear  as  y  decreases  until  y+  values  of  about  7  to  10.  In 
figure  5.9  the  correlation  coefficient  reaches  a  value  of  0.4  at  a  y+  value  of  40  and  in 
figure  5.10  the  ratio  of  v'  fu'  remains  linear  as  y  decreases  until  a  y+  value  of  about 
20.  While  these  values  are  about  three  times  larger  than  the  y+  values  reported 
by  Barlow,  the  friction  velocity  for  this  flow  is  also  about  three  times  larger  than 
Barlow’s  flow.  Therefore,  the  two  trends  appear  to  occur  at  the  same  location  in 
actual  physical  space.  Since  this  distance  is  on  the  order  of  several  beam  diameters, 
the  limitation  may  be  physical  interference  between  the  test  surface  and  the  trans¬ 
mitted  beams.  Given  this  qualification  or  stipulation,  the  correlations  between  u' 
and  v1  shown  in  figures  5.9  and  5.10  indicate  that  the  cross  terms  and  correlations 
are  valid. 

5.4  Measurement  of  Turbulent  Triple  Products 

The  new  data  acquisition  system  provided  the  capabilty  of  collecting  the  third- 
order  statistics,  even  for  phase-averaged  unsteady  conditions,  without  any  penalties 
to  the  other  data  collection  requirements.  Since  these  measurements  were  available 
basically  “free  of  charge”,  they  were  collected  for  whatever  information  they  could 
provide.  Klewicki  and  Falco  have  recently  published  results  on  the  times  required  for 
statistical  convergence  of  various  turbulent  quantities.  In  general,  they  indicate  that 
the  third  order  statistics  take  four  times  longer  to  converge  than  the  second  order 


statistics.  The  steady  triple  product  measurements  of  the  terms  normally  found  in 
a  two-dimensional  boundary  layer  appear  to  be  nearly  converged,  and  are  ir  general 
agreement  with  results  published  by  Murlis,  Tsai,  and  Bradshaw  (1982)  and  Pauley 
and  Eaton  (1988).  These  results,  shown  in  figures  5.11  and  5.12,  differ  only  in 
that  they  show  peak  values  near  the  surface  not  shown  in  the  other  experiments. 
The  location  of  these  peak  values,  however,  seems  consistent  with  the  fact  that 
maximum  production  occurs  very  near  the  surface.  Finally,  for  unsteady  flows  this 
program  accepted  a  compromise  on  convergence  of  the  second  order  statistics  and 
had  to  do  without  the  triple  products  because  of  the  unacceptably  long  sampling 
times  required  for  convergence. 

5.5  Uncertainty  Analysis 

Estimates  of  the  uncertainly  in  each  ot  the  measured  and  calculated  quantities 
were  made,  and  are  summarized  in  Table  5.1.  These  estimates  are  given  at  the  95 
percent  confidence  level  according  to  the  recommendations  of  Kline  and  McClintock 
(1953)  and  follow  the  recent  procedural  recommendations  of  Moffat  (1988).  The 
details  of  this  process  are  shown  in  Appendix  A.  Because  the  uncertainty  analysis 
requires  insightful  judgement,  guidance  was  also  sought  from  similar  recent  pro¬ 
grams  here  at  Stanford,  in  particular  Brereton  and  Reynolds  (1987),  Pauley  and 
Eaton  (1988),  and  Eaton  and  Johnston  (1980). 
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Table  5.1  Uncertainty  Estimates  for  Measured  Variables 


Variable 

Typical  Value 

(U) 

1.0 -10"1 

(U) 

1.0-10"° 

(V) 

3.0  •  10"2 

(V) 

1.5  •  10"1 

(W) 

3.0  •  10“2 

(W) 

5.0  •  10"2 

(u'u') 

1.5  •  10"2 

( v'v ') 

1.5  •  10'2 

(w'w1) 

1.5  •  10“2 

{ u'v ') 

2.5  •  10"3 

(u'xv1) 

5.0  •  10"3 

(t l'w') 

2.5  •  10“3 

Uncertainty  Interval 


Absolute 

Relative 

6.3  •  10~3 

6.3  •  10~2 

1.5  •  10"2 

1.5  -10"2 

to 

o 

1 

N3 

3.9-  10"1 

1.5  •  10"2 

9.8  •  10"2 

<N 

1 

o 

w 

3.9  •  10"1 

1.5  •  10"2 

2.9-  10"1 

3.7  •  10"3 

2.5  •  10"1 

3.7  •  10"3 

2.5  ■  10"1 

3.7  •  10“3 

2.5  •  10"1 

p 

V-* 

o 

1 

3.6  •  10"1 

1.8  •  10"3 

3.6-  10"1 

9.1  • 10"4 

3.6 -10"1 
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Figure  5.1  Streamwise  velocity  history 
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Figure  5.9  Correlation  coefficient  for  uV  at  zero  pr  ssure  gradient 


Figure  5.10  Ratio  of  v'/u1  at  zero  pressure  gradient 
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Figure  5.11  u'u'u1  /U^  at  zero  pressure-gradient 
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Figure  5.12  u'u'v1  /U^,  u'v'v1  /U^,  and  v'v'v' /U^.  at  zero  pressure-gradient 


6.  Experimental  Plan 


This  chapter  will  describe  the  vortex  generator  test  apparatus,  the  manipulation 
of  the  freestream  velocity  to  create  the  desired  unsteady  pressure  gradient,  selec¬ 
tion  of  a  realistic  parameter  space  for  the  experimental  problem,  and  some  of  the 
procedures  necessary  to  obtain  a  quality  data  set. 

6.1  Vortex  Generator  Test  Apparatus 

The  vortex  generator  test  apparatus  inserts  into  the  forward  portion  of  the  tur¬ 
bulent  test  surface  and  supports  an  array  of  six  half-delta-wing  counter-rotating 
vortex  generators.  While  the  design  scheme  allowed  for  computer-controlled  actu¬ 
ation  of  the  vortex  generators,  in  this  program  the  surfaces  were  manually  rotated 
as  a  set  to  the  desired  angle  of  attack.  In  addition,  the  test  plate  has  dye  slots  for 
introducing  dye  into  the  boundary  layer  just  ahead  of  the  leading  edge  of  the  vortex 
generators. 

Each  vortex  generator  surface  is  a  half-delta-wing  with  a  chord  length  of  3cm 
and  a  height  of  1.2cm.  Details  of  the  vortex  generator  geometry  are  shown  in  figure 
6.1  The  half-delta-wing  waj  chosen  primarily  because  that  shape  has  produced 
satisfactory  results  in  other  programs  at  Stanford.  The  dimensions  were  based 
on  recommendations  by  Pearcey  (1961)  for  an  array  of  generators  with  th  '  best 
combination  of  spacing  and  range: 

D/d  =  4 
D/h  -  10 

t/h  =  2.5.  (6.1) 

While  Pearcey’s  recommendations  tend  to  favor  a  co-rotating  array,  the  counter¬ 
rotating  array  was  chosen  because  the  vortices  tend  to  convect  straight  along  the 


69 


channel.  In  contrast,  co-rotating  arrays  tend  to  “walk”  sideways  because  of  the 
effect  of  the  image  vortices. 

Once  a  counter-rotating  array  was  chosen,  the  number  of  pairs  of  generators  was 
basically  fixed  by  the  size  constraints  of  the  water  channel.  The  pairs  were  to  be 
spaced  so  that  the  effect  of  image  vortices  would  simulate  an  infinite  array  (see 
figure  6.2).  In  addition,  an  odd  number  of  pairs  was  chosen  so  that  one  pair  would 
be  located  on  the  channel  center  line  for  ease  in  LDA  measurements.  Since  with 
five  vortex  generator  pairs  the  size  of  each  generator  surface  would  have  been  too 
small  relative  to  the  boundary-layer  thickness,  the  logical  choice  was  three  vortex 
generator  pairs.  With  three  pairs  scaled  according  to  equations  6.1,  the  1.2cm  height 
of  each  vortex  generator  compares  favorably  to  a  nominal  momentum  thickness  of 
0.4cm  and  a  boundary-layer  thickness  of  4cm. 

The  basic  mounting  plate  for  the  vortex  generators  is  a  26.6cm  by  11.4cm  lucite 
insert  which  fits  into  the  turbulent  test  surface.  Figure  6.3  shows  a  photograph  of 
the  flow  side  of  the  insert  plate  with  the  vortex  generators  set  at  an  angle  of  attack 
of  about  18°,  and  figure  6.4  shows  the  opposite  side  of  the  test  plate,  revealing  the 
control  mechanism.  Each  generator  blade  fits  into  a  slot  in  a  12.7mm  stainless-steel 
shaft.  A  sprocket  on  the  opposite  end  of  the  shaft  rotates  the  assembly.  Figure  6.5 
shows  a  cross-section  schematic  of  the  shaft  assembly.  Pump  seals  are  used  to  allow 
rotation  while  preventing  leakage.  Each  shaft  is  supported  by  two  nylon  bearings, 
one  in  the  lucite  plate  and  one  in  the  aluminum  sprocket  mount  plate.  This  plate 
also  is  used  to  compress  the  pump  seals  by  means  of  adjustable  rings  which  allow 
each  shaft  to  be  individually  positioned  flush  with  the  flow  surface. 

Figure  6.3  also  reveals  that  there  are  four  separate  dye  slots  adjacent  to  each 
other,  each  with  its  own  plenum.  These  allow  introduction  of  dye  onto  all  or  just 
part  of  the  surface,  or  the  opportunity  to  use  different  color  dyes  to  study  the  mixing 
between  vortices. 
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Figure  6.4  shows  the  drive  chain  assembly  used  to  rotate  the  generators  to  the 
desired  angle  of  attack.  A  steel-reinforced  polyurethane  belt  engages  the  sprocket  on 
each  shaft  in  a  serpentine  manner,  resulting  in  the  counter-rotation  of  each  vortex 
pair.  This  particular  approach  was  chosen  over  a  variety  of  alternate  designs  because 
the  chain  and  sprocket  assembly  ideally  has  no  backlash,  and  because  the  same 
drive  approach  works  very  successfully  on  the  computer-controlled  valves  described 
in  chapter  4.  In  practice,  the  eccentric  loads  on  each  shaft  result  in  too  much  friction 
for  this  design  to  work  well  for  rapidly  actuated  vortex  generators.  Nevertheless, 
this  design  was  quite  satisfactory  for  manually  setting  the  set  of  vortex  generators 
to  a  desired  angle  of  attack.  Checks  for  repeatability  indicated  that  the  angle  of 
attack  of  each  generator  can  be  set  to  within  about  one  degree  of  the  angle  set  for 
the  array. 

6.2  Tailoring  the  Freestream  Flow  and  Pressure  Gradient 

Since  the  stated  goal  of  this  program  was  to  study  the  effect  of  vortex  gener¬ 
ators  placed  in  an  unsteady,  separating,  turbulent  boundary  layer,  the  computer- 
controlled  valves  modify  the  freestream  velocity  such  that  the  tr  oulent  boundary 
layer  separates  downstream  of  the  vortex  generators  when  they  ,.re  at  an  angle  of 
attack  of  0°.  A  family  of  steady  streamwise  velocity  gradients  capable  of  creat¬ 
ing  this  separation  are  shown  in  figure  6.6.  This  figure  shows  four  different  steady 
velocity  gradients.  The  first  shows  constant  velocity  or  a  zero  pressure  gradient 
condition.  In  the  subsequent  three  cases,  a  pair  of  valves  are  manipulated  together 
to  create  the  adverse  pressure  gradients.  The  upstream  valve  (located  between 
1.95m  and  2.05m)  was  opened  while  a  downstream  valve  was  closed  a  proportional 
amount  such  that  the  velocity  at  the  leading  edge  remained  constant  at  0.5 m/sec. 
As  this  suction  through  the  upstream  valve  is  increased,  the  velocity  and  pressure 
gradients  become  steeper  until,  in  the  fourth  case,  the  boundary  layer  separates  at 
about  2.05m. 
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For  an  unsteady,  adverse  pressure  gradient,  the  two  valves  oscillate  with  the 
same  waveform  and  frequency  but  180°  out  of  phase.  By  monitoring  the  upstream 
and  downstream  velocities  on  a  dual-trace  oscilloscope,  the  relative  amplitudes  and 
phases  of  the  two  valves  can  be  adjusted  and  balanced  for  the  exact  flow  conditions 
desired.  While  the  facility  has  the  capability  of  imposing  the  pressure-gradient  in 
virtually  any  time  waveform  desired,  a  square  wave  was  chosen  for  this  experiment 
so  that  frequency  was  not  a  parameter.  The  corners  of  the  square  wave  were  slightly 
rounded  over  0.2  seconds  of  the  10-second  cycle  (71%  change  in  0.05  seconds,  91% 
change  in  0.1  seconds,  98%  change  in  0.15  seconds,  and  lOOin  0.2  seconds).  Figure 
6.7  shows  the  resulting  unsteady  adverse  pressure  gradient.  Note  that  at  tUo/l  =  0, 
there  is  no  pressure  gradient.  At  tUo/l  =  1  the  gradient  is  the  strongest,  and  then 
as  the  tunnel  adjusts  to  the  blockage  caused  by  the  separation,  the  gradient  relaxes 
a  small  amount. 

The  location  of  the  adverse  pressure  gradient  with  respect  to  the  location  of 
the  vortex  generators  was  chosen  for  optimum  interaction  between  the  streamwise 
vortices  and  the  separating  boundary-layer.  It  was  also  important  that  the  vortex 
generators  were  located  in  a  region  of  steady  flow  such  that  the  streamwise  vortices 
convected  into  the  region  of  unsteady  flow.  Laser  sheet  flow  visualization  in  steady 
flow  gave  a  good  indication  of  the  size  and  strength  of  the  vortex  structures  that 
are  formed  by  the  vortex  generators.  Figure  6.8  shows  spanwise  cross-sections  of 
the  structures  behind  the  vortex  generators  at  various  streamwise  positions  and  for 
three  blade  positions:  18°  (common  flow  up),  0°,  and  18°  (common  flow  down).  The 
common  flow  down  case  at  a  location  60cm  downstream  of  the  vortex  generators 
appears  to  be  optimum  in  terms  of  a  vortex  that  is  not  too  diffuse  but  has  grown 
to  where  the  structures  fill  most  of  the  boundary  layer.  Therefore,  the  location  of 
the  adverse  pressure  gradient  was  chosen  so  that  separation  tends  to  occur  about 
60cm  downstream  of  the  vortex  generators.  Figure  6.9  summarizes  the  relationships 
between  the  locations  on  the  test  surface  of  the  vortex  generators  and  the  imposed 
pressure  gradient. 
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6.3  Upstream  Conditions. 


This  experiment  was  designed  such  that  the  velocities  varied  in  a  near  linear 
fashion  along  the  pressure  gradient  ramp,  and  such  that  they  were  steady  in  time 
upstream  of  the  ramp.  To  verify  this,  and  to  document  the  upstream  conditions  for 
computational  purposes,  unsteady  profiles  were  taken  at  the  streamwise  positions 
x  =  1.88m  and  x  =  2.00m.  The  resulting  profiles  for  the  separation  portion  of  the 
cycle  are  shown  in  figures  6.10  and  6.11. 

6.4  Spanwise  Uniformity 

Chapter  5  demonstrated  that  the  spanwise  uniformity  of  the  flow  was  very  good 
under  steady  conditions.  Unsteady  flow  acts  as  an  amplifier  and  greatly  exagger¬ 
ates  any  inconsistencies  in  the  flow.  The  same  can  be  said  of  an  adverse  pressure 
gradient.  Therefore,  one  would  expect  that  it  is  far  more  difficult  to  maintain  good 
spanwise  uniformity  under  unsteady  adverse  pressure  gradient  conditions.  Figure 
6.12  shews  the  spanwise  uniformity  of  the  streamwise  velocity  component,  (u),  with 
vortex  generators  set  at  a  =  0°,  and  figure  6.13  shows  similar  results  for  a  =  18°. 
In  both  cases,  at  tU$/l  =  0  the  velocity  is  quite  uniform  across  the  span.  Also, 
at  tUo/l  =  1,  the  initial  response  to  the  suction  is  quite  uniform  in  both  cases. 
However,  as  the  flow  adjusts  to  the  separation  and  to  the  resulting  olockage  in  the 
tunnel,  some  non-uniformity  does  develop.  Figure  6.12  indicates  the  pressure  gra¬ 
dient  relaxes  more  on  the  positive  z  axis,  resulting  in  a  spanwise  gradient  in  the  —  z 
direction.  This  effect  is  even  stronger  with  the  vortex  generators  set  to  a.  —  18°. 
While  the  reason  for  this  effect  was  not  determined,  the  spanwise  flow  caused  by 
this  spanwise  gradient  is  seen  very  clearly  in  the  secondary  flow  vectors  presented 
in  chapter  7.  Finally,  figures  6.14  through  6.16  show  the  unsteady  spanwise  uni¬ 
formity  of  (v),  (uV),  and  (t/V)  respectively.  Note  that  while  ( u'u ')  and  (v'v1)  do 
show  larger  fluctuations  on  the  —  z  side,  which  is  consistent  with  a  stronger  adverse 
pressure  gradient,  that  {v)  does  not  show  much  non-uniformity.  This  suggests  that 
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the  suction  created  by  the  automated  valve  is  very  uniform,  and  that  the  valves  are 
not  the  cause  of  the  spanwise  variation. 

6.5  Experimental  Approach 


This  chapter  has  already  identified  most  of  the  experimental  parameters  of  the 
experiment.  The  dimensions  of  the  vortex  generators  have  been  defined  and  the 
characteristics  of  the  unsteady  adverse  pressure-gradient  have  been  explained.  The 
question  that  remains  is  where  to  take  data.  The  basic  experimental  objective  was 
to  completely  map  out  a  symmetrical  spanwise  section  of  the  unsteady  flow  field  at 
a  single  streamwise  position,  so  that  mean  and  turbulent  transport  terms  could  be 
estimated. 

Streamwise  Location.  The  streamwise  location  selected  for  the  data  plane 
was  simply  determined  by  how  close  the  three-component  optics  could  be  located 
on  the  downstream  side  of  the  point  of  separation.  Because  of  a  support  rib  on 
the  top  of  the  test  section,  and  because  of  the  water  tank  required  for  optical  path 
length  compensation  on  the  top  of  the  tunnel,  the  closest  location  was  x  —  2.08 m. 
The  most  significant  length  for  the  experiment  is  the  distance  downstream  of  the 
pressure  ramp,  which  commences  at  x  =  1.95 m.  Therefore  t  =  13cm,  as  indicated 
in  figure  6.9.  Also,  this  defines  a  useful  time  scaling  parameter,  the  convection  time 
based  on  the  initial  freestream  velocity.  For  convenience  in  scaling  time  only,  since 
Uq  —  0.5m/ sec,  the  length  will  be  taken  as  12.5cm  so  that 

i  _  0.125m 

U0  0.50 m/sec 


—  0.25sec 


(6.2) 


Time  Period  of  Unsteady  Cycle.  The  factor  that  most  directly  affected  how 
many  spatial  data  points  to  take  was  the  time  required  for  ^ach  unsteady  phase- 
average.  This,  of  course,  depended  on  the  time  period  of  each  cycle  and  the  number 
of  cycles  required  for  reasonable  convergence  of  the  turbulence  quantities.  Pre¬ 
liminary  tests  at  various  locations  in  the  flow  indicated  that  most  of  the  “mean” 
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velocity  change  occurred  in  about  0.5see,  and  most  of  the  change  in  the  turbulence 
quantities  occured  in  about  2sec,  but  some  adjustment  occurred  over  the  next  sev¬ 
eral  seconds.  Therefore,  in  order  for  the  flow  to  not  have  the  characteristics  of 
the  driving  frequency,  a  ten  second  cycle  was  chosen.  Since  the  valve  control  soft¬ 
ware  had  not  been  set-up  to  accommodate  diilerent  times  for  the  separation  and 
recovery  phases,  this  meant  there  was  a  five-second  separation  phase  followed  by  a 
five-second  recovery  phase. 

Convergence  Requirements.  Convergence  requirements  were  discussed  briefly 
in  chapter  5,  but  the  study  referenced  there  by  Xlewicki  and  Falco  (1988)  did  not 
address  the  added  complication  of  phase-conditioned  sampling.  Using  the  ten- 
second  cycle  discussed  above,  preliminary  data  sets  were  taken  at  selected  locations 
in  the  boundary  layer,  for  phase-conditioned  averages  based  on  200,  500,  1000,  and 
3000  cycles.  The  200-sample  averages  showed  reasonably  good  convergence  of  the 
“mean”  velocities,  but  very  poor  convergence  of  the  Reynolds  stress  terms.  The 
500-sample  averages  showed  very  good  convergence  of  the  “mean”  velocities  and 
reasonable  convergence  of  the  turbulence  quantities.  As  would  be  expected,  the 
1000-sample  and  3000-sample  averages  showed  proportionately  better  convergence. 
The  minimum  acceptable  number  of  cycles,  therefore,  was  500.  While  more  cycles 
would  have  been  desired,  a  compromise  was  made  because  of  the  time  constraints 
discussed  in  the  following  section.  This  compromise  was  possible  because  time 
history  information  was  used  to  make  up  for  the  incomplete  convergence  of  the 
turbulence  quantities  using  phase  smoothing  of  the  data,  a  concept  which  will  be 
presented  in  section  7.2. 

Data  Grid  Selection.  The  sampling  requirements  defined  in  the  previous  two 
paragraphs  result  in  a  sampling  time  requirement  of  5000  seconds,  or  one  hour  and 
23  minutes  at  each  individual  point  in  space,  and  two  vortex  generator  postions  were 
to  be  run  at  each  point.  Faced  with  these  unpleasant  realities  and  the  awareness 
that  equipment  malfunctions  do  occur  in  any  complex  facility,  a  less  dense  data 
grid  was  selected  than  would  have  been  selected  for  a  steady  experiment.  The 
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minimum  spanwise  requirements  were  to  map  between  the  center  line  and  the  line 
of  symmetry  between  vortex  generator  pairs.  (In  terms  of  the  vortex  generator 
parameters  defined  in  figure  6.1,  0  <  z/D  <  0.5.)  Across  the  span,  five  spanwise 
locations  were  selected  as  the  minimum  realistic  number  that  could  capture  the  key 
features  of  the  flow.  In  the  y  direction,  6  logarithmically  spaced  points  were  selected 
to  map  the  inner  region  of  the  boundary  layer,  and  12  evenly  spaced  points  were 
selected  to  map  out  the  outer  region.  Once  again,  this  seemed  to  be  a  minimum 
reasonable  number  that  could  capture  most  of  the  structure  of  the  flow  field.  Table 
6.1  summarizes  the  basic  data  grid  for  this  experiment. 

The  data  grid  of  Table  6.1  was  followed  in  a  manner  that  provided  the  greatest 
breadth  of  information  as  early  as  possible  in  the  data  acquisition  process.  Each 
spanwise  position  required  a  major  realignment  of  the  w  optics,  so  a  complete 
data  set  was  taken  at  each  spanwise  position.  At  each  position,  a  sequence  of  one 
or  more  quick  steady  profiles  was  taken  to  verify  that  the  u  and  v  optics  were 
properly  aligned,  paying  special  attention  to  the  ttV  correlation  coefficient  as  a 
sensitive  indicator  of  good  alignment.  Next  the  w  optics  were  aligned  to  the  u  and 
w  optics.  Then,  a  documentation  set  of  steady  profiles  were  taken  u  —  0°  and 
a  =  18°.  Finally,  if  all  the  cteady  data  looked  good,  the  unsteady  profiles  were 
taken.  Allowing  for  automated  traverse  time,  each  unsteady  profile  took  about  26 
hours.  The  data  acquisition  process  was  checked  every  several  hours,  and  could  be 
stopped  and  restarted  as  necessary. 
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Table  6.1  Primary  Data  Grid 


Spanwisc  position,  z/D 


alpha 

-1/2 

-1/4 

0 

1/8 

1/4 

3/8 

1/2 

Steady 

0° 

21 

21 

21 

21 

21 

21 

21 

Steady 

18° 

21 

21 

21 

21 

21 

21 

21 

Unsteady 

0° 

6 

6 

18 

18 

18 

18 

18 

Unsteady 

18° 

6 

6 

18 

18 

18 

18 

18 

Nofe:  The  number  in  each  column  indicates  the  number  of  points  in  that  particular 
profile 
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J 

D  = 1 2  cm 

d  =  3  cm 

J,  □ 

h  =  1 .2  cm 

vg 

'vg  =  3  cm 

Figure  6.1  Dimensions  of  vortex  generator  assembly 


Vortex  Generator  Configuration 
in  Water  Channel 


Figure  6.2  Effect  of  image  vortices 
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A  -  vortex  generator  blade 
B  -  drop-in  test  surface 
C  -  set  screw 
D  -  sprocket  shaft 
E  -  nylon  bearing 
F  -  pump  seal  seat 
G  -  pump  seal  head 
H  -  stainless-steel  washer 
I  -  retaining  ring 


K  -  sprocket  mount  plate 
L  -  nylon  bearing 
M  -  lock  ring 
N  -  adjustment  ring 
0  -  cable  chain 
P  -  chain  sprocket 
Q  -  stainless-steel  nut 
R  -  set  screw 
S  -  shaft  coupling 


Figure  6.5  Cross-section  of  vortex  generator  assembly 
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Figure  6.6  Steady  streamwise  veloctiy  gradients 


Streamwise  position  (cm) 

Figure  6.7  Unsteady  streamwise  velocity  gradients 
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Streamwise  Common  Flow 
Position  (cm)  Down,  20° 


Common  Flow 
Up,  20° 
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Figure  6.8  Streamwise  development  of  embedded  streamwise  vorticity 
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Figure  6.9  Location  of  vortex  generators  and  adverse  pressure  gradient 
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Figure  6.10  Unsteady  upstream  profile  at  x 
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Spanwise  position,  z/D 

Figure  6.12  Spanwise  uniformity  of  unsteady  (u)  at  a  = 
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Figure  6.13  Spanwise  uniformity  of  unsteady  (u)  at  a 


(u'lO/Uo2  at  y/h=5 


Spanwise  position,  z/D 


7.  Results  I:  Unsteady  Measurements 


This  chapter  will  present  the  basic  phase-averaged  measurements  of  velocity  and 
Reynolds  stresses.  Ttie  major  analysis  and  evaluation  of  the  results,  however,  will 
be  reserved  for  chapter  8.  While  steady,  zero  pressure-gradient  profiles  were  taken 
as  part  of  the  set-up  validation  process  at  each  spanwise  location,  these  results  are 
usually  redundant  with  the  unsteady  results  obtained  at  t  =  0-.  Therefore,  steady- 
flow  data  will  not  be  shown  in  this  chapter.  However,  some  steady  results  will  be 
discussed  in  chapter  8  as  a  point  of  departure  from  other  steady  experiments  with 
streamwise  vorticity. 

The  key  features  of  the  experimental  plan  described  in  chapter  six  were  that  the 
primary  unsteady  data  for  this  experiment  consists  of  two  5  by  18  grids  at  a  single 
streamwise  location,  one  grid  for  the  vortex  generators  at  an  angle  of  attack  of  0° 
and  one  at  an  angle  of  attack  of  18°.  The  grid  consisted  of  5  spanwise  positions 
and  18  positions  normal  to  the  surface.  Each  unsteady  cycle  consisted  of  a  (l)  step 
change  from  zero  pressure  gradient  to  adverse  pressure  gradient,  (2)  5  seconds  of 
constant  adverse  pressure  gradient,  (3)  a  step  change  back  to  zero  pressure  gradient, 
and,  finally,  (4)  a  5-second  recovery  period  at  zero  pressure  gradient.  During  each 
cycle,  the  data  acquisition  system  maintained  phase  averages  at  400  separate  phases, 
that  is,  every  0.025sec.  While  the  second  half  of  the  cycle  produced  information  on 
the  boundary  layer  recovery,  all  the  results  discussed  in  this  report  pertain  only  to 
the  separation  half  of  the  cycle. 

It  is  a  challenge  to  know  just  how  to  display  adequately  both  time  development 
and  spatial  relationships  in  the  data.  Since  the  data  was  measured  at  one  point  in 
space  at  a  time,  it  seems  most  natural  to  begin  with  the  time  histories  of  the  various 
quantities  at  selected  spatial  locations.  This  then  leads  logically  to  an  explanation 
of  how  data  at  a  given  point  in  time  was  assembled  over  an  entire  spatial  plane. 
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Finally,  the  chapter  will  conclude  with  limited  flow  visualization  results  to  reinforce 
the  information  conveyed  by  the  data. 

7.1  Time  History  of  Single  Point  Measurements 

Figures  7.1  through  7.5  present  the  time  history  of  the  phase-averaged  quantities 
(u),  (v),  { u'u '),  ( v'v '),  and  —(u'v1)  for  the  first  two  seconds  of  the  separation  cycle. 
Consistent  with  the  non-dimensionalization  described  in  chapter  six,  one  second  is 
equivalent  to  4  convection  times,  IUq /l.  The  two-second  time  length  was  chosen  for 
display  since  most  of  the  significant  changes  in  the  flow  have  occurred  by  this  time. 
In  each  figure,  three  locations  of  y/h  are  shown,  one  which  is  fairly  representative 
of  the  inner  region  of  the  boundary  layer,  one  representative  of  the  inner  portion  of 
the  wake,  and  one  representative  of  the  outer  portion  of  the  wake. 

The  time  history  curves  were  created  by  joining  all  the  phase-averaged  data 
points  collected  between  t  =  —0.25 sec  to  t  =  2.0 sec.  Since  the  sampling  program 
maintained  a  separate  phase-averaged  bin  every  0.025  seconds,  there  are  actually 
90  separate  phase-averaged  data  points  represented  in  each  trace.  Data  symbols 
were  eliminated  from  these  plots  for  clarity.  These  time  histories  are  included  to 
demonstrate  the  relative  magnitude  of  change  in  each  quantity  as  well  as  the  relative 
time  over  which  the  changes  occurred. 

In  general,  the  mean  velocities  respond  much  more  quickly  than  the  turbulent 
stresses.  With  respect  to  (u),  there  is  a  fast  response  which  occurs  in  the  same 
time  frame  as  the  opening  of  the  control  valve,  followed  by  a  slower  adjustment 
period.  In  general,  the  velocity  continues  to  decrease  during  the  adjustment  period, 
but  in  certain  regions  of  the  downwash,  it  actually  increases  during  this  adjustment 
period.  The  normal  velocity  component,  (v),  responds  very  little  in  the  inner  re¬ 
gion  of  the  boundary  layer.  In  addition  to  responding  much  more  slowly  than  the 
mean  velocities,  the  turbulent  stresses  decrease  in  the  inner  boundary  layer  as  the 
separation  develops,  but  increase  substantially  in  the  outer  region  of  the  boundary 


layer.  The  time  response  of  {v) ,  (") ,  and  {a'u')  will  be  characterized  in  considerable 
detail  in  chapter  8. 

7.2  Phase-Smoothing  of  the  Data 

The  time  histories  shown  in  figures  7.1  through  7.5  also  convey  information  about 
the  degree  of  convergence  obtained  for  each  quantity  during  the  500  cycles  over 
which  phase  averages  were  taken.  The  quantity  (u)  seems  tc  be  fully  converged 
and  ( v )  seems  to  be  almost  converged.  In  contrast,  the  three  turbulence  quantities 
shown  are  not  fully  converged.  As  discussed  in  chapter  6,  the  number  of  spatial 
points  and  number  of  averaging  cycles  chosen  required  some  compromises  in  order 
to  obtain  a  reasonable  data  set  in  a  reasonable  time  period.  However,  to  overcome 
this  lack  of  convergence,  one  can  take  advantage  of  the  fact  that  the  phase-averaged 
quantities  should  evolve  “smoothly”  in  time. 

For  this  experiment,  a  piece-wise  parabolic  least-squares  fit  was  used  to  approxi¬ 
mate  the  expected  “smooth”  behavior  of  each  variable.  In  general,  9  or  11  adjacent 
phases  were  fit  using  a  NAPL  Singular  Value  Decomposition  (SVD)  subroutine, 
which  produced  a  least-squares  parabolic  fit.  (NAPL  is  the  Numerical  Analysis 
Package  Library  developed  by  Argonne  National  Laboratory.)  The  central  time 
value  of  the  phases  being  smoothed  was  then  substituted  into  the  least-squares  fit. 
The  value  obtained  from  that  fit  was  then  used  as  a  “phase-smoothed”  value  for 
most  of  the  remaining  data  presentation  and  analysis.  While  this  process  is  a  data 
smoothing  process,  in  many  respects  it  is  simply  an  intelligent  extension  of  the 
averaging  process  used  in  all  turbulence  studies. 

To  demonstrate  the  validity  of  this  technique,  figure  7.6  shows  time  histories  of  (u) 
and  —(u'v1)  with  the  values  obtained  by  phase-smoothing  shown  as  symbols  plotted 
over  the  time  traces.  Figure  7.6a  demonstrates  that  the  phase-smoothing  process 
yields  values  that  are  consistent  with  the  well-converged  values  of  (u).  In  a  similar 
fashion,  figure  7.6b  shows  that  the  phase-smoothed  values  are  also  representative  of 
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what  one  would  expect  for  the  converged  -\-V)  observations.  Since  this  process 
yields  credible  values  for  both  converged  and  non-converged  data,  it  was  applied  to 
all  the  remaining  data  in  this  chapter,  and  will  be  referred  to  as  “phase-smoothed” 
data. 

7.3  Spanwise  presentations  of  the  data 

While  the  time  history  data  just  presented  were  useful  in  seeing  how  the  velocities 
and  Reynolds  stresses  at  a  point  evolve  in  time,  one  is  usually  also  interested  in  how 
these  quantities  vary  in  profile  format,  both  in  the  direction  normal  to  the  surface 
and  in  the  spanwise  direction.  Therefore,  using  phase-smoothing  to  improve  the 
convergence  of  the  data,  the  velocities  and  Reynolds  stresses  were  reassembled  into 
complete  spanwise  planes  of  data  for  a  single  instant  in  time  and  for  a  specific 
position  of  the  vortex  generator. 

Frequently,  spanwise  sets  of  data  are  presented  as  contour  plots.  Contour  plots 
work  best  for  dense  data  grids,  but  several  disadvantages  of  contour  plots  are:  (1) 
the  actual  data  points  are  not  shown,  (2)  splining  through  data  before  contour 
plotting  can  produce  some  bizarre  and  anomalous  features  on  the  contours,  and  (3) 
data  smoothing  effects  are  concealed  from  the  reader  (and  from  the  data  analyst  if 
he  is  not  careful).  Since  time  requirements  dictated  a  rather  coarse  data  grid,  the 
results  from  this  experiment  seemed  particularly  vulnerable  to  the  weaknesses  of 
contour  plotting.  Therefore,  in  this  report,  contour  plotting  will  be  minimized.  An 
alternative  display,  which  will  be  used  in  all  the  remaining  figures  in  this  chapter, 
is  a  three-dimensional  perspective  view  of  the  profiles  from  each  of  the  five  primary 
spanwise  locations. 

In  each  of  the  figures,  from  figures  7.7  through  figure  7.32,  four  complete  spanwise 
data  sets  are  presented.  The  upper  two  data  sets  represent  cases  with  the  vortex 
generators  at  an  angle  of  attack  of  0°  and  the  lower  two  data  sets  represent  cases 
with  the  vortex  generators  at  an  angle  of  attack  of  18°.  Also,  the  two  data  sets  on 
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the  left  represent  a  single  phase  in  time,  and  the  two  sets  on  the  right  represent 
a  later  single  phase  in  time.  In  this  way,  each  figure  provides  the  reader  with  a 
comparison  to  a  different  phase  of  the  cycle  and  a  comparitive  view  of  the  effect 
of  the  vortex  generators.  In  most  cases,  the  curve  drawn  is  a  5-point  piecewise 
parabolic  least-squares  fit  to  the  data.  For  simplicity  in  reading,  all  the  figures 
from  7.7  to  7.32  follow  this  identical  pattern. 

Streamwise  velocity  component.  Figures  7.7  through  7.10  present  the  de¬ 
velopment  of  the  streamwise  velocity  component,  (u).  At  OIq /l  =  0,  the  profiles 
are  quite  full  both  with  and  without  vortex  generators.  In  contrast  to  the  stronger 
vortices  studied  by  Pauley  and  Eaton  (1988),  these  profiles  do  not  show  a  large 
velocity  deficit  in  the  vortex  core.  A  closer  evaluation  of  the  boundary  layer  in¬ 
tegral  parameters  in  chapter  8  will  highlight  other  differences  in  these  profiles.  In 
succeeding  figures,  the  profiles  are  seen  to  react  erratically  to  the  intial  change  in 
the  pressure  gradient.  Without  the  vortex  generators,  the  boundary  layer  begins  to 
separate  at  tU0/l  —  6.  Since  this  experiment  was  intentionally  designed  to  create 
a  very  mild  separation,  the  boundary  layer  remains  slightly  attached  at  z/D  =  0 
and  at  z/D  =  0.5,  but  reverse  flow  occurs  at  the  other  three  locations.  With  vor¬ 
tex  generators  on,  the  flow  remains  attached  in  the  downwash  regions,  zjD  =  0, 
z/D  =  0.125,  and  z/D  =  0.25,  but  the  boundary  layer  separates  very  quickly  in 
the  upwash  regions  of  the  vortex,  z/D  —  0.375  and  z/D  —  0.5.  Reverse  flow  first 
begins  at  z/D  =  0.5  at  tUo/l  =  3,  and  by  tUo/l  =  8  the  separated  region  extends 
three  generator  heights  away  from  the  surface.  Additional  turbulent  profiles  in  wall 
coordinates  (u+  vs.  y+),  will  be  presented  in  chapter  8,  along  with  discussions 
regarding  the  skin  friction  and  friction  velocity. 

Secondary  velocities.  Figures  7.11  through  7.14  show  the  analogous  behavior 
of  the  secondary  flow,  which  is  best  seen  as  the  combined  vector  sum  of  (v)  and 
(u/).  Because  the  zero-velocity  offset  voltage  is  very  hard  to  get  precisely  correct, 
small  errors  in  (v)  and  (w)  show  up  as  large  variations  in  the  secondary  flow.  For 
this  reason,  at  each  location  the  values  of  (u)  and  ( w )  at  tUo/l  —  0  and  a  —  0° 
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were  subtracted  as  the  offset  error.  The  small  uniform  vectors  shown  in  figure  7.11 
at  this  reference  state  are  intended  to  communicate  that  the  secondary  flow  there 
is  assumed  to  be  small  and  uniform.  At  tUo/l  =  0  and  a  =  18°,  the  secondary 
flow  associated  with  the  vortex  pair  is  clearly  seen.  The  downwash  and  upwash 
regions  are  both  evident,  and  the  vortex  core  appears  to  be  at  about  zjD  =  0.32 
and  y/h  =  1.2.  Also,  the  vortex  is  quite  weak  in  comparison  to  the  vortex  strength 
used  in  most  studies  in  the  literature.  The  maximum  secondary  velocity  is  about 
only  five  percent  of  the  streamwise  velocity.  At  a  =  0°,  as  the  separation  develops, 
(v)  increases  dramatically  as  the  boundary  layer  separates.  A  small  component  of 
—  ( w )  is  also  seen  at  a  =  0°  as  the  separation  develops.  This  is  consistent  with  the 
observations  in  chapter  6  that  the  tunnel  tends  to  pull  in  the  —  z  direction  during  the 
adverse  pressure  gradient  portion  of  the  cycle.  At  a  =  18°,  the  circular  secondary 
flow  associated  with  the  vortex  can  be  seen  until  tUo/l  —  3.  After  this,  the  vortex 
seems  to  lose  its  identity  and  some  large  spanwise  flows  are  seen  in  conjunction  with 
the  separation.  The  significance  of  these  latter  observations  will  be  dealt  with  at 
the  end  of  chapter  8. 

Reynolds  stresses.  Figures  7.15  through  7.32  show  the  development  of  each  of 
the  components  of  the  Reynolds  stress  tensor.  The  normal  stresses,  (u'u,))  (v'v1), 
( w'w '),  and  the  primary  shear  stress  —  (uV),  all  behave  in  the  the  same  general 
manner.  The  trends  are  qualitatively  similar  to  those  seen  by  Pauley  and  Eaton 
(1988).  At  a  =  0°,  the  peak  concentration  of  each  stress  tends  to  spread  outward 
away  from  the  surface  as  the  boundary  layer  separates.  Also,  the  magnitude  of  each 
stress  increases  during  the  adverse  pressure  gradient.  At  a  =  18°,  the  concentration 
of  each  stress  is  pushed  inwards  in  the  downwash  region  and  redistributed  outwards 
in  the  upwash  region.  This  was  described  well  by  Mehta  and  Bradshaw  (1988) 
as  “convection  of  pre-existing  boundary-layer  turbulence  by  the  vortices.”  As  the 
separation  develops,  the  magnitude  of  each  stress  increases  even  more  at  a  =  18° 
than  the  increase  seen  at  a  =  0°.  This  is  particularly  true  of  ( u'u '}  and  (v'v1). 
Note  that  during  the  separation  at  a  =  18°  the  shear  stress,  —  (uV),  even  develops 
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regions  of  negative  stress  in  the  region  of  zjD  =  0.5.  With  regard  to  the  secondary 
shear  stresses,  ( u'w ')  and  ( v'w '),  very  little  if  any  should  be  seen  at  a  =  0°,  and  very 
little  is  seen.  There  is  a  significant  concentration  of  (u'w1)  in  the  cross-flow  region 
of  the  vortex,  as  would  be  expected.  This  stress  is  greatly  amplified  during  the 
adverse  pressure  gradient.  In  addition,  the  small  concentration  of  (v'w1)  at  a  —  18° 
is  also  amplified  significantly  during  the  separation. 

7.4  Flow  Visualization 

To  add  some  additional  insight  into  the  behavior  of  the  streamwise  vortices, 
movies  were  taken  of  the  vortices  during  an  unsteady  cycle,  using  laser  sheet  fluo¬ 
rescence.  Water  containing  fluoroscein  sodium  salt  (C2oHioNa2C>5)  was  introduced 
through  the  dye  slots  just  upstream  of  the  vortex  generators.  The  fluid  containing 
this  dye  tended  to  wrap  around  the  vortex  generator  and  mark  the  region  of  con¬ 
centrated  vorticity.  A  sequence  of  four  photos  from  this  movie  is  shown  in  figures 
7.33  through  7.36.  In  all  the  photos  the  bright  line  is  the  laser  sheet  illuminating 
the  surface  and  the  “bar”  seen  is  a  support  rib  on  the  outer  side  of  the  channel. 
Figure  7.33  shows  the  fluid  marked  by  the  dye  during  the  zero  pressure  gradient 
phase.  These  structures  should  resemble  those  shown  in  figure  6.8.  Figure  7.34 
shows  the  initial  response  to  the  applied  suction.  While  for  the  moment  the  struc¬ 
tures  generally  have  not  changed,  one  streak  of  fluid  marked  by  the  dye  is  being 
pulled  out  of  the  boundary  layer.  Figure  7.35  shows  the  fluid  in  the  structures  now 
being  pulled  away  from  the  boundary  layer  until  finally  in  figure  7.36,  the  marked 
fluid  now  extends  about  four  times  its  original  distance  from  the  surface. 
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Figure  7.12  Secondary  velocity 
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Figure  7.22  Profiles  of  ( w'w ')  at  tUo/l  =  1.2  and  HJq/1 


tlL/l  =  8.0 


vo  =  i/tu  oo  =  i/ni 


Profiles  of  (u'w1)  at  Wq/1  =  0  and  tUQ/l  =  0.4 


tU„/l  =  1.2  tl)n/l  =  2.0 


Profiles  of  { u'w ')  at  Wq/I  =  1.2  and  tU^/l 


FO 


tu„/l  =  4.0  tiyi  =  8.0 


Profiles  of  { v'w ')  at  Wq/I  —  4  and  Wq/1 


Figure  7.33  Flow  visualization  at  Wq/I  =  0 


Figure  7.34  Flow  visualization  at  Wq/1  =  O1- 
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Figure  7.35  Flow  visualization  at  Wq/1  —  1 


Figure  7.36  Flow  visualization  at  Wq/I  =  3 
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Remember,  when  discoursing  about  water, 
to  induce  first  experience,  then  reason. 
—  Leonardo  da  Vinci 
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8.  Results  II:  Analysis  and  Discussion 


Chapter  7  fulfilled  one  program  objective,  which  was  to  document  profiles  of  the 
velocities  and  Reynolds  stresses.  This  chapter  will  address  the  remaining  objectives, 
which  are  to  characterize  the  time  response  of  the  boundary  layer  and  to  understand 
the  mechanisms  by  which  vortex  generators  modify  its  response.  The  time  response 
characteristics  that  will  be  presented  offer  useful  information  for  the  design  of  a 
flow  control  system  and  provide  insight  into  the  mechanisms  of  flow  modification. 

As  a  starting  point  for  the  analysis,  some  results  at  steady,  zero  pressure  gradi¬ 
ent  conditions  will  be  compared  to  similar  results  from  other  steady  experiments. 
Then  the  time  response  characteristics  of  the  integral  parameters  and  of  (u),  (v), 
and  ( u'u ')  will  be  discussed.  The  variation  in  skin  friction  coefficient  and  behavior 
of  the  mean  velocity  in  wall  coordinates  are  then  presented  for  documentation  pur¬ 
poses.  Next,  evaluation  of  the  transport  of  streamwise  vorticity,  spanwise  vorticity, 
and  momentum  provide  additional  insight  into  the  mechanisms  that  make  vortex 
generators  work.  Finally,  a  look  at  turbulent  kinetic  energy  production  and  the 
anisotropy  of  the  Reynolds  stress  tensor  provides  some  insight  into  the  response  of 
the  turbulence. 

8.1  Steady  Flow:  Comparison  to  previous  experiments 

Skin  Friction  and  Mean  Turbulent  Profiles.  Use  of  the  universal  law  of  the 
wall  lends  insight  into  the  two-dimensionality  of  the  boundary  layer  and  provides  a 
means  to  estimate  the  spanwise  skin  friction  variation  with  and  without  vortex  gen¬ 
erators.  Eibeck  and  Eaton  (1985)  concluded  that  their  weak  vortices  only  affected 
the  wake  region  of  the  boundary  layer,  whereas  Pauley  and  Eaton  (1988)  in  addition 
saw  significant  deviation  in  the  logarithmic  region.  Steady  zero  pressure  gradient 
profiles  were  taken  as  part  of  the  initial  validation  at  each  spanwise  position.  These 
profiles  are  presented  in  wall  coordinates  in  figures  8.1a  and  8.1b. 
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For  all  the  analysis  in  this  chapter,  steady  and  unsteady,  the  skin  friction  at  the 
surface  was  estimated  using  a  best  fit  of  the  near-wall  data  to  the  universal  law  of  the 
wail.  To  accomplish  this,  several  different  approaches  were  evaluated:  (l)  a  fit  to  the 
“inner  law”,  u+  =  y+,  (2)  a  fit  to  the  outer  region  following  the  recommendations 
of  Coles  (1968),  and  (3)  a  fit  to  the  entire  region  covered  by  the  universal  law  of 
the  wall  using  interactive  graphics  and  an  approximate  expression  for  the  law  of  the 
wall.  Of  the  three  approaches,  the  third  gave  the  most  consistent  and  satisfactory 
results.  Use  of  the  “inner  law”  was  unsuccessful  because  not  enough  data  points 
were  taken  in  the  linear  region.  The  Coles  method,  using  computer  codes  developed 
by  Jayaraman  and  Brereton,  gave  good  results  for  the  attached  boundary  layer 
without  vortex  generators.  However,  the  vortex  generators  sufficiently  distorted 
the  boundary  layer  that  the  least-squares  fit  totally  misinterpreted  the  location  of 
the  logarithmic  region  and  wake.  Fortunately,  interactive  graphics  enable  one  to 
make  a  very  good  visual  determination  of  a  best  fit. 

For  this  experiment,  the  interactive  graphics  package  Quattro  was  used.  ( Quattro 
is  a  registered  trademark  of  Borland  International,  Inc.)  The  universal  law  of  the 
wall  was  expressed  in  analytical  form  using  an  approximate  implicit  expression 
developed  by  Spalding  (1961), 

y+  =  «+  +  e~kcleku+  -  1  -  ku+  -  l{ku+)2  -  i(lfcu+)3  -  i-(fcu+)4},  (8.1) 

l  2  6  24  J 

where  k  =  0.41  and  c  =  5.0.  Then,  allowing  for  a  small  correction  in  the  location 
of  the  wall  (on  the  order  of  1  y+  unit)  a  best  fit  of  the  inner  wall  data  is  obtained 
visually  by  varying  the  skin  friction  (or  friction  velocity)  and  the  offset  from  the 
wall.  For  boundary  layers  with  no  embedded  streamwise  vorticity,  this  method 
was  compared  to  the  Coles  method  using  the  computer  codes  by  Jayaraman  and 
Brereton.  The  comparison  showed  that  the  visual  method  is  very  sensitive  and 
repeatable,  yielding  friction  velocities  consistent  with  the  Coles  method  to  about 
2  percent.  The  visual  “universal  law”  method  also  works,  however,  for  boundary 
layers  with  highly  distorted  wake  regions,  where  the  Coles  method  does  not.  (Note: 
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A  more  useful  expression  than  equation  8.1  is  given  by  Liakopoulos  (1984).  This 
was  discovered  after  completing  the  analysis  described  above.) 

The  resulting  values  for  steady  flow  of  the  skin  friction  coefficient  and  the  friction 
velocity  are  shown  in  figures  8.2a  and  8.2b.  The  skin  friction  coefficient,  Cf,  and  the 
friction  velocity,  uT,  are  each  related  to  the  wall  shear  stress,  rw,  by  their  definitions, 


Relative  to  a  value  of  Cf  of  0.004  without  vortex  generators,  Cf  varies  from  about 
0.0047  at  the  centerline  to  about  0.0052  in  the  strongest  downwash,  and  to  about 
0.003  in  the  upwash  region.  These  values  are  similar  in  magnitude  and  in  the  degree 
of  variation  to  the  “common  flow  down”  study  of  Pauley  and  Eaton  (1988)  and  the 
“common  flow  up”  study  of  Mehta  and  Bradshaw  (1988).  In  general,  the  variation 
in  Cf  was  less  in  this  experiment  than  in  the  other  experiments  referenced.  This  is 
consistent  with  all  observations  that  indicate  the  induced  vortex  in  this  experiment 
was  quite  weak  relative  to  other  experiments.  Of  the  cases  shown  by  Pauley,  the 
variation  in  Cf  is  more  similar  to  his  downstream  case  where  the  variation  from 
minimum  to  maximum  values  is  a  little  less  than  a  factor  of  2.  It  is  also  important 
to  note  that  in  the  other  experiments  which  were  all  in  air,  the  skin  friction  was 
actually  measured,  whereas  in  this  experiment  it  was  only  inferred  from  the  velocity 
profiles. 

The  velocity  profiles  in  figure  8.1a  indicate  a  nearly  uniform  wake  region  with 
the  vortex  generators  at  an  angle  of  attack  of  0°,  whereas  major  variations  are  seen 
in  the  wakes  of  figure  8.1b  with  the  vortex  generators  at  an  angle  of  attack  of  18°. 
The  spanwise  variations  seen  in  figure  8.1a  are  greater  than  spanwise  variations  seen 
before  the  vortex  generator  assembly  was  installed  in  the  tunnel,  and  are  attributed 
to  the  streamwise  vorticity  and  total  pressure  wake  that  is  formed  around  the  vortex 
generators  even  at  an  angle  of  attack  of  0°.  With  the  vortex  generators  at  18°, 
however,  the  entire  wake  region  is  severely  distorted.  As  would  be  expected,  the 
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upwash  region  has  a  large  amount  of  sluggish  fluid  and  manifests  a  large  wake. 
Conversely,  the  downwash  and  cross-flow  regions  have  more  energetic  fluid  near 
the  surface  and  show  very  little  evidence  of  any  wake.  Since  the  linear  spacing  of 
data  points  resulted  in  very  few  points  in  the  logarithmic  region,  it  is  hard  to  be 
dogmatic  about  how  much  this  region  is  distorted  by  the  vortex  generators.  The 
evidence,  though,  seems  to  indicate  that  even  with  these  weak  vortex  pairs,  the 
inner  portion  of  the  boundary  layer  is  not  altogether  two-dimensional,  consistent 
with  the  observations  of  Pauley  and  Eaton. 

The  Integral  Parameters.  The  integral  parameters  displacement  thickness, 
6*,  momentum  thickness,  6,  and  shape  factor,  H ,  are  useful  since  they  summarize 
in  a  single  number  the  condition  of  the  entire  boundary  layer.  Standard  definitions 
of  these  quantities  were  used,  but  are  restated  below  for  clarity. 

6 

0 

<8-4> 

J  Uoo  v  Uoo  / 

0 

H  =  j.  (8.5) 

Th».  ; y*i.«.vise  var??.*»cn  or  the  displacement  thickness  and  momentum  thickness  is 
shown  in  figure  8.3a  and  the  spanwise  variation  of  the  shape  factor  is  shown  in 
figure  8.3b.  In  all  cases  where  integration  was  performed  in  this  program,  the 
integrands  were  first  plotted  to  ensure  that  they  were  smooth,  continuous  functions 
with  no  anomalies  at  the  origin  or  boundary- layer  edge.  Then,  a  parabolic  numerical 
integration  scheme  was  used,  which  allowed  for  the  uneven  spacing  of  data  points 
near  the  wall. 
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The  momentum  thickness,  with  vortex  generators  at  an  angle  of  attack  of  0°, 
varies  across  the  span  from  3.8mm  to  4.2mm.  As  with  the  variation  in  Cj,  this 
variation  is  presumed  to  be  caused  by  vorticity  generated  by  the  vortex  generators  at 
a  =  0°.  Using  an  average  value  for  6  of  4mm,  along  with  the  freestream  velocity  of 
0.5 m/sec,  and  the  kinematic  viscosity  of  water  at  16.7°C,  the  momentum  thickness 
Reynolds  number,  Reg,  is  determined  to  be  1840  ±  90.  This  number  will  be  used 
to  characterize  the  state  of  the  boundary  layer  at  the  start  of  the  adverse  pressure 
gradient  cycle.  The  variation  in  the  integral  parameters  with  the  vortex  generators 
at  an  angle  of  attack  of  18°  is  very  much  as  would  be  expected.  In  the  downwash  and 
crossflow  regions  the  shape  factor  is  about  1.2,  which  is  well  below  the  typical  values 
of  1.4  to  1.5  for  this  range  of  Reynolds  numbers.  This  is  indicative  of  an  energetic 
boundary  layer  with  very  little  sluggish  flow.  In  contrast,  in  the  upwash  region 
the  momentum  thickness,  displacement  thickness,  and  shape  factor  all  increase, 
indicating  that  sluggish  fluid  from  near  the  surface  has  been  moved  into  the  outer 
portion  of  the  boundary  layer. 

Transport  of  Turbulent  Kinetic  Energy.  As  part  of  the  validation  of  the 
data  acquisition  systems,  chapter  5  presented  the  triple  product  measurements  for  a 
two-dimensional  boundary  layer.  Shabaka,  Mehta,  and  Bradshaw  (1985)  presented 
their  triple  product  data  in  terms  of  transport  velocities  of  normal  and  shear  stresses. 
Pauley  and  Eaton  (1988)  presented  their  results  in  a  similar  fashion.  In  general,  data 
collection  times  during  this  experiment  were  not  long  enough  for  good  convergence 
of  the  triple  products.  This  was  particularly  true  for  t u'w1  statistics.  Also,  the 
coarser  data  grid  selected  for  this  study  does  not  provide  good  resolution  on  a  vector 
map.  Nevertheless,  secondary  transport  velocities  of  the  quantitity  u'u'  +  v'v1  were 
computed  and  are  presented  in  figure  8.4,  to  demonstrate  that  a  similar  transport 
phenonenon  was  observed  with  this  experiment.  The  secondary  transport  velocities 
were  computed  from  the  following: 

_  (u'uV)  +  (v'v'v1) 

~  (u'u')  +  (v'v1) 
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Figure  8.4  shows  the  same  general  qualitative  trend  as  reported  by  Pauley  and  Eaton 
(1988)  for  their  “common  flow  down”  vortex  pair.  It  also  reveals  the  limitations 
of  the  triple  product  data  taken  as  part  of  this  experiment.  Therefore,  no  further 
data  analysis  or  interpretations  will  be  based  on  triple  product  results. 


8.2  Time  Response  of  the  Boundary  Layer 


Chapter  7  presented  the  characteristic  time  history  of  the  mean  velocities  and  the 
Reynolds  stresses.  This  section  will  now  characterize  the  nature  of  the  boundary- 
layer  response  to  the  step  change  in  pressure  gradient.  The  analysis  will  show  that 
there  is  an  immediate  response  to  the  pressure  gradient  throughout  the  boundary 
layer  which  is  characteristic  of  the  elliptical  nature  of  the  governing  equations, 
followed  by  a  slow  or  convective  response  characteristic  of  the  hyperbolic  nature  of 
the  governing  equations.  This  section  will  begin  with  contour  plots  of  the  streamwise 
velocity  followed  by  the  time  development  of  the  integral  parameters,  then  will 
present  the  parameters  which  describe  the  step  response  in  (u),  (v),  and  (u'u1). 

Contours  of  Streamwise  Velocity.  Contour  plots  of  the  streamwise  velocity 
component  provide  a  good  visualization  of  the  changes  in  the  boundary  layer.  Fig¬ 
ures  8.5  through  8.9  show  contours  of  (u),  with  vortex  generators  at  0°  and  18°,  as 
the  boundary  layer  responds  to  the  impulsively  started  adverse  pressure  gradient. 
The  “jagged  corners”  on  the  contour  plots  result  from  the  coarse  data  grid,  and  no 
attempt  was  made  to  conceal  this.  Figure  8.5  shows  the  nature  of  the  boundary 
layer  with  zero  pressure  gradient.  While  it  is  thinned  in  the  downwash  region  and 
thickened  in  the  upwash  region,  there  is  not  a  large  velocity  deficit  apparent  in  the 
core  of  the  embedded  vortex,  as  was  seen  in  the  experiments  of  Pauley  and  Eaton 
(1988).  In  this  respect,  this  vortex  more  resembles  the  experiments  of  Bradshaw  et 
al.,  where  the  vortex  generators  were  placed  before  the  converging  nozzle  of  their 
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wind  tunnel.  Succeeding  figures  show  the  growth  of  the  boundary  layer  as  the  ad¬ 
verse  pressure  gradient  is  introduced.  In  figure  8.6,  a  small  region  of  reverse  flow 
is  already  apparent  in  the  upwash  region.  In  figure  8.7,  reverse  flow  occurs  across 
most  of  the  span  without  vortex  generators  in  addition  to  the  upwash  region  with 
vortex  generators.  This  trend  then  continues  in  figures  8.8  and  8.9.  This  sequence 
of  figures  follows  the  development  of  the  boundary  layer  through  Wq/1  =  8.  Since 
all  of  the  salient  features  of  the  flow  adjustment  process  occur  in  that  time  frame, 
most  of  the  analysis  and  discussion  that  follows  will  cover  this  same  time  period. 

Response  of  the  Integral  Parameters.  The  integral  parameters  are  one 
means  of  summarizing  the  condition  of  the  entire  boundary  layer.  It  was  impossible 
to  compare  the  time  responses  of  the  integral  parameters  across  the  entire  span  be¬ 
cause  the  boundary  layer  grew  in  the  upwash  region  beyond  the  traverse  capabilities 
of  the  three-component  LDA  system.  The  integrals  in  equations  8.3  and  8.4  only 
have  meaning  if  the  integrand  goes  to  zero  on  the  outer  edge  of  the  boundary  layer. 
Since  this  was  not  the  case  across  most  of  the  span,  the  time  response  of  the  inte¬ 
gral  parameters  is  shown  only  for  the  center  line  position.  Figure  8.10a  shows  the 
response  of  the  displacement  thickness  and  the  momentum  thickness,  while  figure 
8.10b  shows  the  response  of  the  shape  factor. 

If  a  simple  time  constant  for  these  integral  parameters  is  defined  as  the  length 
of  time  required  for  63%  of  the  change  to  occur,  then  the  time  constant  of  the 
momentum  thickness,  in  both  cases,  is  on  the  order  of  1.51  /Uq.  The  same  is  true  for 
the  displacement  thickness  when  the  vortex  generators  are  at  an  angle  of  attack  of 
18°,  except  there  is  a  secondary  adjustment  which  occurs  over  a  longer  time  period 
but  is  small  in  magnitude.  Since  the  freestream  velocity  is  decelerating  from  Uq 
to  0.6C/qj  this  value  of  the  time  constant  of  1.51 /Uq  is  approximately  equal  to  the 
time  of  flight  of  fluid  in  the  freestream.  In  contrast  to  this,  the  time  constant  for 
the  displacement  thickness  when  the  vortex  generators  are  at  an  angle  of  attack 
of  0°  is  approximately  31/Uq.  The  behavior  of  the  shape  factor  in  figure  8.10b 
represents  the  combined  effect  of  the  changes  in  the  displacement  thickness  and 
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momentum  thickness,  which  are  occurring  on  different  time  scales.  The  analysis 
of  the  streamwise  velocity  time  response,  will  provide  more  information  into  the 
nature  of  the  boundary-layer  response. 

Time  Response  of  the  Streamwise  Velocity.  Figures  7.1a  through  7.1f 
showed  the  time  history  of  the  streamwise  velocity.  In  each  case,  there  is  an  initial 
response  during  which  the  velocity  decreases  quite  rapidly,  followed  by  a  slow  ad¬ 
justment  period  during  which  the  velocity  continues  to  change.  While  in  most  cases 
the  velocity  continues  to  decrease,  it  actually  increases  in  portions  of  the  downwash 
region.  The  response  of  (u)  will  be  modelled,  therefore,  in  a  two-step  fashion  as 
depicted  by  figure  8.11. 

The  Mathematical  Model.  There  are  various  ways  that  the  velocity  change  can 
be  normalized.  After  experimenting  with  several  different  schemes,  the  following 
was  chosen. 


$  =  Aoo  =  U0  -  Uf.  (8.8) 

.Aoo 

In  equation  8.8,  Ago  represents  the  amplitude  of  the  freestream  disturbance, 
which  for  this  experiment  was  0.4.  This  suggests  that  the  linearized  system  response 
could  be  represented  as 


$  =  1  - 

or  equivalently, 


(8.9a) 


uo  -  u  =  Aoo(l  -  r°^-  (8.96) 

The  justifications  for  this  normalization  scheme  are:  (1)  this  follows  the  general 
form  of  the  solution  to  Stokes’  first  problem  (except  the  actual  solution  is  a  com¬ 
plimentary  error  function),  (2)  the  freestream  disturbance,  A<x>,  can  be  viewed  as 
a  system  input,  (3)  this  scheme  avoids  ambiguities  that  arise  when  the  change  is 
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normalized  on  the  total  local  change  because  the  change  is  not  always  monotonic 
in  time,  and  finally  (4)  this  scheme  in  general  provided  a  very  good  collapse  of  the 
data. 

To  demonstrate  the  actual  behavior  of  $  relative  to  the  mathematical  model 
shown  in  figure  8.11,  representative  values  of  $  have  been  plotted  in  figure  8.12a 
and  8.12b  for  the  centerline  position,  zj  D  =  0.  Note  that  these  figure  show  the  same 
key  features  as  the  mathematical  model.  These  features  are:  (l)  a  region  where  the 
response  is  very  quick,  characterized  by  an  exponential  rate,  /?,  and  a  time  delay, 
to;  (2)  an  initial  roll-off  where  the  slope  decreases  drastically,  characterized  by  the 
time  of  the  roll-off,  r2,  and  the  corresponding  value,  $r;  and  (3)  a  slow  adjustment 
period  where  the  value  of  $  may  increase  or  decrease,  characterized  by  the  final 
value,  and  a  time  constant,  73,  which  is  the  time  to  reach  0.63(4>^  —  $r). 

Determination  of  the  Time  Response  Parameters.  Due  to  the  smooth,  well- 
converged  nature  of  the  streamwise  velocity  data,  a  straight-forward  computerized 
scheme  was  used  to  identify  the  time  response  parameters  at  each  y  and  2  location. 
A  segment  of  the  curves  from  figures  8.12a  and  8.12b  are  shown  again  in  figures  8.13a 
and  8.13b  with  data  symbols  to  demonstrate  that  there  are  sufficient  data  points 
along  the  quick  rise  portion  of  each  curve  to  properly  determine  /?.  Basically,  the 
computer  scheme  marched  along  the  data  until  $  changed  by  5%.  The  program 
then  computed  the  slope  between  each  pair  of  points  and  identified  the  steepest 
slope.  This  steepest  slope  was  then  fit  to  an  exponential  expression  from  which  /? 
and  tq  were  determined.  An  alternate  routine  was  also  tested,  which  used  a  3-point 
least-squares  fit.  It  gave  similar  results,  but  actually  had  more  scatter  than  the  first 
scheme  which  simply  used  the  two  points  with  the  steepest  slope.  To  determine  the 
location  of  the  roll-off  point,  72,  the  computer  scheme  computed  predicted  values 
based  on  a  linear  fit  through  the  two  points  with  the  steepest  slope.  The  value 
of  T2  was  defined,  after  some  experimentation,  as  the  point  where  this  predicted 
value  exceeded  the  actual  value  of  the  data  by  0.2.  Finally,  the  computer  program 
marched  through  the  remainder  of  the  data  and  determined  the  maximum  value  of 
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the  quantity  f  —  $r),  then  went  back  and  located  the  time,  r^,  where  the  value 
was  0.63($ f  —  $r). 

The  resulting  values  of  these  parameters  are  plotted  in  figures  8.14a  through 
8.14f.  Figure  8.14a  presents  the  roll-off  value,  $r.  Note  that  all  the  data  sets,  for 
all  spanwise  positions  and  both  angles  of  attack,  collapse  onto  the  same  curve.  The 
vortex  generators  do  not  appear  to  affect  this  feature  of  the  response.  Due  to  the 
normalization  scheme  chosen,  the  outer  portion  of  the  boundary  layer  rolls  off  at 
a  value  of  $r  of  1,  while  close  to  the  surface  the  roll-off  value  is  proportionately 
smaller.  In  contrast  to  $r,  which  was  not  affected  by  the  vortex  generators,  the 
values  of  the  quantity  —  $r)  as  shown  in  xigure  8.14b,  are  very  clearly  affected 
by  the  vortex  generators.  While  the  value  of  this  quantity  is  quite  consistent  for  the 
three  spanwise  positions  at  an  angle  of  attack  of  0°,  it  is  substantially  reduced  in  the 
downwash  region  and  substantially  increased  in  the  upwash  region.  These  trends, 
in  conjunction  with  the  time  scales  that  follow,  provide  insight  into  the  mechanisms 
that  make  vortex  generators  work. 

Figures  8.14c  through  8.14f  present  the  time  scales  for  the  step  response  of  (u). 
Figure  8.14c  shows  that  the  value  of  the  exponential  rate  for  the  initial  response, 
/?,  is  about  the  same  for  all  data  sets,  and  varies  from  a  value  of  about  -4.5  near 
the  surface  to  about  -7  in  the  outer  portions  of  the  boundary  layer.  Due  to  the 
finite  rise-time  of  the  control  valve,  the  freestream  response  should  have  a  value 
of  /?  of  -8.81.  Within  experimental  uncertainty,  it  appears  that  the  values  of  /? 
computed  from  the  data  are  approaching  this  value  at  the  edge  of  the  boundary 
layer.  Once  again,  this  feature  of  the  response  does  not  apear  to  be  affected  by 
vortex  generators.  On  the  other  hand,  the  value  of  the  initial  phase  delay,  to,  does 
vary  significantly  between  data  sets,  but  the  trends  are  not  well  defined.  Figure 
8.14d  shows  that  the  phase  delay  is  greater  in  the  outer  portion  of  the  boundary 
layer.  The  variation  in  the  phase  delay  at  an  angle  of  attack  of  0°  is  probably  due  to 
the  opening  characteristics  of  the  flow  control  valve.  There  may  be  a  subtle  trend 
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indicating  that  the  vortex  generators  reduce  the  magnitude  of  the  phase  delay  but 
more  data  would  be  required  to  substantiate  this. 

Figure  8.14e  shows  that  the  value  of  fa  —  tq)  is  constant  at  a  value  of  about 
0.351 /Uq  through  the  boundary  layer  for  all  data  sets.  This  indicates  that  the 
duration  of  the  initial  response  is  the  same  throughout  the  boundary  layer  and  is 
unaffected  by  the  presence  of  vortex  generators.  Since  90%  of  the  change  in  the 
control  valve  position  occurs  in  0.41 /Uq,  it  appears  that  the  waveform  introduced  by 
the  control  valve  propagates  quite  uniformly  through  the  boundary  layer.  This  trend 
is  consistent  with  the  observations  from  figure  8.14a  and  8.14c.  That  is,  because 
the  rate  of  change  is  less  in  the  inner  region  while  the  duration  of  the  change  is 
the  same  throughout  the  boundary  layer,  a  lower  roll-off  value  is  seen  in  the  inner 
region.  Finally,  figure  8.14f  shows  that  the  the  time  constant  of  the  convective 
response,  73,  generally  has  a  value  of  about  31 /Uq,  except  in  the  outer  portion  of 
the  boundary  layer  where  the  vortex  generators  do  cause  some  variations.  Note  that 
the  large  variations  in  the  value  of  tq  in  the  inner  regions  of  the  boundary  layer, 
shown  in  figure  8.14d,  do  not  appear  to  affect  T3.  This  is  additional  evidence  to 
suggest  that  the  time  scale  of  the  convective  or  hyperbolic  response  is  independent 
of  the  time  scale  of  the  initial  or  elliptical  response. 

Conclusions  regarding  ( u }.  (1)  Vortex  generators  do  not  affect  the  basic  response 
rates  of  the  boundary  layer.  They  only  affect  the  magnitude  of  the  convective 
portion  of  the  response.  (2)  There  is  an  initial  response  throughout  the  boundary 
layer  typical  of  the  elliptical  nature  of  the  governing  equations,  which  is  unaffected 
by  the  presence  of  vortex  generators.  The  duration  of  this  initial  response  reflects 
the  waveform  of  the  control  valve  and  is  uniform  throughout  the  boundary  layer, 
but  because  the  rate  of  change  is  somewhat  less  in  the  inner  region,  the  magnitude  of 
the  initial  change  is  also  proportionately  less  in  the  inner  region.  (3)  There  is  also  a 
slow  or  convective  response  throughout  the  boundary  layer,  typical  of  the  hyperbolic 
nature  of  the  governing  equations,  during  which  the  magnitude  is  substantially 
modified  by  the  presence  of  vortex  generators.  (4)  Because  the  magnitude  of  the 
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convective  response  is  substantially  reduced  in  the  downwash  region,  a  far  greater 
percentage  of  its  response  occurs  during  the  initial  response.  The  apparent  effect 
of  vortex  generators  in  the  downwash  region,  therefore,  is  that  they  substantially 
decrease  the  time  required  for  the  boundary  layer  to  adjust  to  the  step  change  to 
an  adverse  pressure  gradient. 

Time  Response  of  ( v ).  The  mathematical  model  described  in  the  previous 
parag  aphs  and  summarized  in  figure  8.11  was  also  applied  to  the  step  response  of 
{ v ) ,  with  several  small  changes.  These  changes  were:  (1)  the  convective  response  was 
more  complex  for  (v)  than  for  (u),  so  only  the  initial  response  has  been  modelled, 
(2)  due  to  less  convergence  in  the  data,  a  3-point  least-squares  fit  was  necessary, 
and  (3)  based  on  figure  7.2,  a  value  of  0.06  was  used  for  the  disturbance  almplitude, 
yloo-  Only  three  points  could  be  used  in  the  least-squares  fit  because  the  initial 
response  occurs  over  a  period  which  is  spanned  by  only  three  or  four  phase  points. 

Figures  8.15a  and  8.15b  show  the  initial  response  rate  and  phase  delay  for  the 
step  response  of  ( v ).  Once  again,  all  the  data  sets  show  the  same  trend,  so  vortex 
generators  do  not  appear  to  affect  the  initial  response  of  ( v }.  Due  to  the  effect  of  the 
wall,  the  magnitude  of  the  response  in  (v)  is  very  small,  thus  the  response  rate,  (3, 
approaches  0  near  the  surface.  In  the  outer  region,  however,  (3  approaches  the  rise 
time  of  the  control  valve,  as  was  observed  for  (u).  The  phase  delay,  ro,  is  i  datively 
constant  across  the  boundary  layer.  The  larger  values  near  the  surface  are  probably 
not  significant  because  of  the  near-zero  slope  of  the  response  curve.  Therefore,  it 
appears  that  the  initial  response  in  { v )  is  governed  by  the  same  mechanisms  that 
govern  the  initial  response  of  { u ). 

Time  Response  of  ( u'u' ).  From  the  time  histories  shown  in  chapter  7,  the  step 
response  characteristics  of  ( u'u ')  appear  to  be  representative  of  all  the  turbulent 
quantities.  The  turbulent  quantities  show  only  a  single  response  which  occurs  at  a 
slow  rate  relative  the  initial  response  of  the  streamwise  velocity.  For  comparitive 
purposes,  this  response  was  modelled  in  the  same  manner  as  the  initial  response  of 
(u)  and  (u),  except  for  the  following  modifications:  (1)  the  initial  response  actually 
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describes  the  total  response,  (2)  due  to  less  convergence  and  a  slower  response,  a 
9-point  least-squares  fit  was  used,  (3)  based  on  the  maximum  change  that  occurs, 
a  value  of  0.015  was  used  for  the  disturbance  amplitude,  A0 0,  and  (4)  the  absolute 
value  of  $  was  used  because  (u'u1)  decreases  in  the  inner  region  of  the  boundary 
layer  but  increases  in  the  outer  region  of  the  boundary  layer. 

Figures  8.16a  and  8.16b  show  the  response  rate  and  phase  delay  of  (u'u1).  In 
general,  the  response  rate  is  an  order  of  magnitude  smaller  than  the  initial  response 
rate  of  the  streamwise  velocity.  In  the  outer  region  of  the  boundary  layer,  the 
inverse  of  0  would  correspond  to  a  time  constant  in  the  range  2.51 /Uq  to  41/ Uq, 
which  is  about  the  order  of  magnitude  of  the  convective  response  of  (u).  The 
turbulence  does  respond  more  quickly,  however,  in  the  innermost  portion  of  the 
boundary  layer.  Also,  this  rate  may  be  slightly  increased  in  the  downwash  portions 
of  the  embedded  vorticity.  Finally,  figure  8.16b  shows  the  same  typical  values  for 
the  phase  delay  in  the  inner  region  as  were  seen  for  the  initial  phase  delay  of  the 
streamwise  velocity,  but  somewhat  larger  values  in  the  outer  region.  It  appears, 
therefore,  that  the  turbulence  in  the  inner  region  of  the  boundary  layer  begins  to 
respond  at  the  same  time  as  the  initial  response  of  (u),  except  at  a  much  slower  rate. 
In  contrast,  the  turbulence  in  the  outer  portion  of  the  boundary  layer  responds  to 
the  growth  of  the  boundary  layer,  which  occurs  on  the  slower  convective  time  scales 
of  the  freestream  velocity. 

8.3  Time  Development  of  Cf  and  Turbulent  Mean  Profiles 

This  section  presents  the  variation  in  the  skin  friction  coefficient  and  the  turbu¬ 
lent  mean  profiles  in  wall  coordinates,  analogous  to  the  results  presented  in  section 
8.1  for  steady  flow.  In  response  to  the  adverse  pressure  gradient,  the  shape  of  the 
each  profile  changes  drastically.  Precisely  the  same  methodology  was  used  to  de¬ 
termine  the  friction  velocity  for  unsteady  flow,  up  to  the  time  where  reverse  flow 
begins.  After  that  time  no  value  of  Cj  could  be  inferred  from  the  universal  law  of 
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the  wall  so  no  turbulent  profiles  are  presented.  The  values  obtained  for  Cj,  with¬ 
out  and  with  vortex  generators,  are  shown  in  figures  8.17a  and  8.17b.  The  same 
information  is  presented  in  figures  8.18a  and  8.18b  in  terms  of  the  friction  velocity, 
ur. 

Variations  in  the  skin  friction.  In  figure  8.17a  for  a  =  0°,  the  skin  friction 
is  initially  quite  uniform  at  tUo/l  =  0.  At  Wq/I  —  2,  reverse  flow  appears  to  be 
imminent  and  has  already  begun  near  the  spanwise  location,  zjD  =  0.125.  The 
initial  “pull-away”  at  this  location  was  seen  very  clearly  in  the  velocity  profiles  of 
figure  7.8.  After  Wq/I  =  2,  reverse  flow  occurs  across  most  of  the  span.  While  there 
is  no  reverse  flow  at  z/ D  =  0  and  z/ D  =  0.5,  C j  drops  from  a  value  of  0.004  to  about 
0.0001.  Figure  8.17b  reveals  that  the  boundary  layer  with  embedded  streamwise 
vorticity  responds  much  differently.  In  the  upwash  region,  reverse  flow  occurs  very 
rapidly  at  zjD  =  0.5,  and  eventually  occurs  after  tUo/l  =  4  at  zjD  —  0.375.  At 
the  remaining  three  spanwise  positions,  there  is  no  tendency  towards  reverse  flow. 
A  spanwise  average  of  these  three  positions  would  show  that  Cj  drops  from  an 
initial  value  of  0.0048  to  about  0.001  at  tUo/l  =  2,  and  then  holds  quite  steady  at 
that  value.  The  spanwise  variations  in  time  at  these  three  locations  are  quite  likely 
secondary  effects  due  to  the  blockage  effect  in  the  tunnel  from  the  upwash  regions. 

Variations  in  the  Mean  Turbulent  Profiles.  Using  the  friction  velocities 
determined  above,  the  phase-averaged  “mean”  turbulent  profiles  are  presented  for 
four  dimensionless  times  in  figures  8.19  through  8.22.  The  observations  for  figures 
8.19a  and  8.19b,  for  tUo/l  =  0,  would  be  virtually  the  same  as  the  observations  made 
for  the  steady  profiles  at  zero  pressure  gradient.  As  the  adverse  pressure  gradient 
affects  the  velocity  profiles  at  subsequent  time  intervals,  the  corresponding  velocity 
profiles  show  several  general  trends.  (1)  At  a  —  0°,  the  profiles  deviate  markedly 
from  the  universal  law  of  the  wall,  and  have  no  discernible  overlap  region.  This  is 
consistent  with  the  observations  of  Jayaraman,  Parikh,  and  Reynolds  (1982)  and 
with  White’s  (1974)  characterization  of  a  reversing  flow.  Even  near  reverse  flow,  the 
profiles  at  different  spanwise  positions  are  quite  consistent  with  each  other,  except  in 
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the  outermost  portion  of  the  boundary  layer.  (2)  At  a  =  18°,  the  profiles  are  similar 
on  the  inner  portion,  but  once  again  do  not  have  an  overlap  region  consistent  with 
the  universal  law  of  the  wall.  There  is  substantial  spanwise  variation  in  the  profile 
shape  past  y+  =  60,  until  at  Wq/I  —  8,  where  the  three  profiles  in  the  downwash 
region  are  surprisingly  similar.  Therefore,  while  the  effect  of  the  vortex  generators 
may  leave  the  inner  portion  of  the  boundary  layer  two-dimensional  at  zero  pressure 
gradient,  as  observed  by  Eibeck  and  Eaton  (1985),  this  is  most  certainly  not  the 
case  in  the  presence  of  a  strong  adverse  pressure  gradient. 

8.4  Evolution  of  Streamwise  Vorticity 

The  analysis  will  now  turn  to  the  transport  equations  for  additional  insight  into 
the  mechanisms  by  which  vortex  generators  modify  the  response  of  the  boundary 
layer.  Previous  studies  with  vortex  generators  have  explored  in  depth  the  stream- 
wise  vorticity  equation,  as  originally  formulated  by  Perkins  (1970),  so  this  seems  to 
be  a  logical  place  to  start.  The  results  will  show  that  after  the  step  change  occurs 
in  the  pressure  gradient,  the  core  of  the  embedded  vortex  diffuses  very  rapidly. 

The  Governing  Equation.  The  streamwise  turbulent  vorticity  equation  is 
obtained  by  expanding  the  terms  of  equation  2.11  or  2.12,  for  the  i  =  1  component 
of  vorticity.  This  yields 

(W*)  +  (tt)(«x),x  +(v)(«x),y  +  ( U/)(W*)„  =  (Wx)(ti),x  +(Wy)(u),y 
+  (Uz){u),z+v((ux),xx+{ux),yy+(ux),zz)  +  (uV),** 


(u  tV  ),xy  +(u  U  )  ,yz  (v  W  )  ,yy  -\-{ti  W  )  >Zz  (w  W  )  >yz  •  (8.10) 

The  three  components  of  the  vorticity  are 

(wx)  =  (w),y-(w)„  (8.11) 

(wy)  =  (u),z-{w),x  (8.12) 

(wz)  =  (y),i  — (u),y  .  (8.13) 

147 


In  equation  8.10,  the  left-hand  side  represents  the  rate  of  change  of  streamwise  vor- 
ticity  plus  the  convection  of  vorticity  out  of  the  fluid  element  by  the  three  phase- 
averaged  velocity  components.  On  the  right-hand  side,  the  first  three  terms  rep¬ 
resent  vortex  stretching  and  bending  and  the  next  three  terms  represent  viscous 
diffusion  into  the  fluid  element.  The  last  six  terms  represent  turbulent  transport 
of  vorticity.  They  are  the  turbulent  counterparts  to  the  convection  and  stretching 
terms,  and  can  be  grouped  in  several  different  ways.  In  equation  8.10  they  are 
grouped  in  pairs  as  the  turbulent  flux  in  the  x,  y,  and  z  directions.  Alternately, 
Perkins  (1970)  grouped  them  as 

P2  =  (( uV '),z  ~{ulw,)ty)ix  (8-14) 

Pi  =  ((*>'*''}  -  (u>V»,j,z  (8.15) 

P\  =  (v'w'),z2  -{t >'w'),yy  .  (8.16) 

In  a  zero  pressure  gradient  flow,  Pauley  and  Eaton  (1988)  found  that  the  terms  P3 
and  P4  were  dominant  in  changing  the  structure  of  embedded  vorticity.  For  this 
flow,  where  there  are  significant  gradients  in  the  streamwise  direction,  any  of  the 
terms  can  act  to  create  or  alter  the  streamwise  vorticity.  Since  this  experiment  was 
not  configured  so  as  to  provide  estimates  of  changes  in  the  streamwise  direction, 
the  emphasis  in  the  analysis  that  follows  will  be  on  what  happens  to  the  streamwise 
vorticity  rather  than  a  detailed  analysis  of  the  individual  transport  terms. 

Differentiating  Data.  A  major  hazard  associated  with  vorticity  estimates  in 
an  experimental  program  is  that  differentiation  of  data  is  required.  To  estimate 
the  vorticity,  a  single  derivative  is  required,  and  to  estimate  terms  in  the  equations 
above,  up  to  three  derivatives  or  the  product  of  two  derivative  terms  are  required. 
The  following  considerations  apply  to  each  of  the  derivative  terms  that  will  be 
discussed  in  this  section  and  the  following  sections.  (1)  A  method  was  selected 
that  found  the  slope  of  the  expected  behavior  (that  is  a  fit  to  the  data)  rather  than 
finding  slopes  of  the  inconsistencies  in  the  data.  For  the  data  in  this  program,  the 
scheme  that  appeared  to  be  most  satisfactory,  was  to  first  make  a  5-point  piecewise 
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parabolic  least-squares  fit  to  the  data,  and  then  use  a  cubic  spline  through  the 
resulting  curve  to  estimate  the  derivative.  (2)  The  fit  of  the  splines  through  the 
smoothed  data  points  was  checked  at  each  stage  of  the  differentiation  to  ensure  that 
the  expected,  behavior  was  indeed  represented.  (3)  Wherever  possible,  the  expected 
behavior  was  compared  to  analytically  predicted  curves  and  baseline  experimental 
curves.  This  was  particularly  important  for  second  derivatives,  where  the  the  shapes 
of  the  curves,  degree  of  curvature,  and  location  of  inflection  points  are  all  key 
features  of  the  final  results. 

Determination  of  ux.  Figures  8.23  and  8.24  show  splines  in  the  y  direction 
through  the  smoothed  values  of  { tv ),  and  figures  8.25  and  8.26  show  splines  in  the  z 
direction  through  the  smoothed  values  of  { v ).  Since  the  splines  are  reasonably  well- 
behaved,  the  derivatives  from  these  splines  were  used  to  calculate  the  streamwise 
vorticity,  uz.  These  figures  also  demonstrate  several  important  features  about  the 
embedded  streamwise  vorticity.  (l)  First,  figure  8.25a  is  plotted  on  the  same  scale 
as  figures  8.23  and  8.24.  Notice  that  in  general  the  y  gradients  of  ( w )  and  the 
z  gradients  of  ( v )  are  of  the  same  general  magnitude.  This  is  characteristic  of 
embedded  vorticity  generated  by  vortex  generators  and  in  contrast  to  the  spanwise 
vorticity  embedded  in  a  boundary  layer,  where  one  of  the  two  gradients  is  much 
larger  than  the  other.  (2)  The  exception  to  this  trend  is  very  near  the  surface  where 
the  { w )  component  is  larger  and  then  must  decrease  rapidly  to  meet  the  “no-slip” 
condition  at  the  surface.  This  large  gradient  of  opposite  sign  in  { w )  will  be  seen  as  a 
strong  vorticity  of  opposite  sign  near  the  surface.  (3)  The  scale  on  subsequent  plots 
of  { v )  was  necessarily  changed  to  accomodate  the  growing  value  of  this  component 
as  the  boundary  layer  grows  rapidly  due  to  the  adverse  pressure  gradient.  (4) 
Figure  8.25a  clearly  identifies  the  location  of  the  core  of  the  embedded  vortex,  as 
the  magnitude  of  the  velocity  component  passes  through  zero  at  a  spanwise  position 
of  zj D  ~  0.32.  (Note,  since  all  derivatives  were  non-dimensionalized  with  respect 
to  the  vortex  generator  height,  z  is  scaled  here  a s  z/h.  Since  the  spacing  between 
vortex  pairs,  D,  is  equal  to  10/*,  the  conversion  back  to  z/D  is  quite  simple.) 
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The  streamwise  vorticity  calculated  from  the  splined  data  shown  above  is  pre¬ 
sented  in  figures  8.27  throught  8.31.  In  this  case,  contour  plotting  was  the  most 
practical  graphical  means  to  present  the  results.  The  contour  plotting  interpolation 
was  done  from  the  original  data  grid,  since  splining  to  create  a  denser  grid  also 
seemed  to  create  fictitious  regions  of  vorticity.  The  contour  information  was  also 
checked  against  tabulations  of  the  vorticity  to  confirm  that  the  contours  shown  were 
realistic  and  consistent.  Figures  8.27  through  8.31  show  the  progression  in  time  of 
the  streamwise  vorticity,  for  the  case  a  =  18°  only.  At  dimensionless  times  of 
tUo/l  =  0,  0.4,  1.2,  and  2,  the  embedded  vortex  (with  negative  vorticity)  is  clearly 
seen.  The  induced  vorticity  of  opposite  sign,  due  to  the  “no-slip”  condition  on  ( w ), 
is  also  seen  very  clearly.  The  concentration  of  positive  vorticity  centered  at  zjh  =  5 
is  generally  accepted  to  result  from  the  convection  of  the  wall-generated  vorticity 
by  the  primary  vortex. 

After  several  convection  times,  the  primary  embedded  vortex  suddenly  loses  its 
strength.  At  tUo/l  =  3  and  4,  the  primary  vortex  appears  to  be  diffused  and  weak¬ 
ened  in  strength.  Past  this  time,  the  negative  vorticity  appears  to  be  “wrapped¬ 
around”  a  vortex  of  opposite  sign  that  appears  from  the  outer  portion  of  the  bound¬ 
ary  layer.  The  weakening  of  the  vorticity  is  the  expected  consequence  of  the  adverse 
pressure-gradient,  shown  by  the  first  term  on  the  right-hand  side  of  equation  8.10, 
{ uix){u),x .  The  adverse  pressure-gradient,  caused  by  a  negative  value  of  (u),z,  “com¬ 
presses”  the  vortex  tube,  making  it  “fatter”  and  weaker.  While  this  same  effect  was 
seen  by  Westphal  et  al.  (1987),  they  concluded  that  the  core  diffused  more  rapidly 
than  could  be  accounted  for  by  the  streamwise  velocity  gradient  term.  Their  anal¬ 
ysis  indicated  that  the  anisotropy  between  the  secondary  normal  stresses,  v'v'  and 
w'w1,  was  responsible  for  the  increased  diffusion  of  vorticity.  Since,  under  unsteady 
conditions,  the  turbulence  responds  quite  slowly  (as  shown  in  section  8.2),  its  role 
in  this  diffusion  process  is  less  clear  than  in  the  case  of  steady  flow. 

To  characterize  the  change  in  vortex  strength,  the  vorticity  of  the  primary  vortex 
was  integrated  to  estimate  the  total  circulation.  First,  the  data  was  splined  to 
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create  a  denser  grid,  then  following  the  practice  of  Westphal  et  al.  the  limits  of 
numerical  integration  were  defined  by  where  the  vorticity  dropped  to  less  than  10 
percent  of  the  peak  value.  The  results  shown  in  figure  8.32  indicate  that,  while  the 
vorticity  spreads  out,  the  circulation  stays  relatively  constant  until  the  vortex  loses 
its  identity  and  integration  becomes  impractical. 

The  source  of  the  positive  vorticity  that  becomes  more  prominent  at  Wq/I  =  6 
is  not  totally  clear.  Since  there  are  gradients  in  all  three  directions,  nearly  all 
the  terms  in  equation  8.10  could  participate.  Since  figures  8.7  through  8.9  show 
that  there  is  a  strong  spanwise  gradient  of  (u),  a  likely  source  of  positive  vorticity 
appears  to  be  the  bending  of  the  spanwise  vorticity  into  positive  streamwise  vorticity 
by  this  strain-rate.  This  is  described  by  the  third  term  on  the  right-hand  side  of 
equation  8.10,  (u2)(u),z.  Howevei,  after  closer  inspection  one  observes  that  part 
of  this  strain-rate  appears  with  opposite  sign  in  the  second  term  of  equation  8.10. 
Expanding  the  two  terms  together  yields: 

(w„)  <u)  ,y  +  (w*)  (u)  „  =  (« „  - u>,* )  <u>  ,y  +  (u - u,y )  («) 

=  {u,z  u,y  —W,x  u,y  -j-V,x  U,z  —U,y  U,z  ) 

=  (~t V,XU,y+V,XU,z). 

Therefore,  it  is  very  difficult  to  identify  which  term  is  the  source  of  the  apparent 
vorticity  of  opposite  sign. 

8.5  Evolution  of  Spanwise  Vorticity 

In  their  review  paper,  Reynolds  and  Carr  (1985)  stressed  the  importance  of 
spanwise  vorticity  in  understanding  the  mechanisms  of  a  boundary  layer  separating 
due  to  an  adverse  pressure  gradient.  The  accumulation  of  spanwise  vorticity  was 
seen  as  the  difference  between  the  convective  inflow  and  outflow  of  vorticity  plus  the 
vorticity  flux  at  the  surface  due  to  the  pressure  gradient.  This  section  will  attempt 
to  identify  the  key  transport  terms  and  how  they  are  modified  by  the  embedded 
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streamwise  vorticity.  The  results  will  show  that  diffusion  and  turbulent  transport  in 
the  y  direction  are  the  dominant  terms,  and  that  the  majority  of  vorticity  transport 
occurs  in  the  viscous  inner  region  of  the  boundary  layer. 

Following  the  procedures  for  differentiating  data  that  were  outlined  in  section 
8.4,  the  spanwise  vorticity  was  estimated  from  the  term  —  (u),y.  That  is 

(u>z)  =  (v)iX  ~{u)iy 
-  ~  (u)>y  • 

The  time  development  of  the  spanwise  vorticity,  thus  determined,  is  shown  in 
figures  8.33  and  8.34.  These  figures  vividly  show  that  a  healthy  attached  boundary 
layer  is  characterized  by  a  very  high  concentration  of  spanwise  vorticity  at  the 
surface.  The  scale  on  these  figures  was  chosen  to  magnify  the  vorticity  in  the 
outer  portion  of  the  boundary  layer  at  the  expense  of  “chopping”  the  value  near 
the  wall.  At  zero  pressure  gradient,  the  vorticity  at  the  surface  peaks  at  about 
lOk/Uo.  This  is  consistent,  as  it  should  be,  with  the  slope  inferred  by  the  fit  to  the 
law  of  the  wall,  u^/iA  In  contrast,  in  the  case  where  reverse  flow  is  occurring,  the 
vorticity  is  very  nearly  uniformly  distributed  throughout  the  entire  boundary  layer. 
Of  course,  part  of  the  negative  vorticity  has  inevitably  diffused  into  the  surface 
due  to  the  adverse  pressure  gradient.  The  task  for  the  vortex  generators  if  they 
are  to  preclude  separation,  therefore,  is  to  force  a  concentration  of  vorticity  back 
into  the  surface  against  the  natural  processes  that  are  occurring.  Evaluation  of  the 
terms  in  the  equation  for  transport  of  spanwise  vorticity  will  lend  insight  into  the 
mechanisms  involved  in  this  process. 

The  Governing  Equation.  The  spanwise  vorticity  equation  is  obtained  by 
expanding  equation  2.12  for  »  =  3.  For  convenience  in  the  analysis  that  follows, 
each  term  will  be  given  a  name  that  also  serves  to  emphasize  the  nature  of  that 
term.  With  that  preface,  the  spanwise  vorticity  equation  is 

(wz)  +  C\  +  C2  +  C3  =  S\  -f-  S2  +  S$  4-  D\  +  D2  +  +  T\  +  T2  +  X3.  (8.17) 
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In  equation  8.17,  the  convection  terms  are 


^1  —  {u){wz),x 
C2  =  {v)(u>z),  y 

C3  =  (t v)(u)g)  tZ  • 

The  vortex  stretching  and  straining  terms  are 

51  -  (w*)(to),* 

52  =  (Uy)(w)ty 

53  =  (Uz)(w),z  ■ 

The  viscous  diffusion  terms  are 

D\  =  i/(ujz),xx 

D2  =  ^M,yy 

D3  —  u(Uz),zz  • 

And  the  turbulent  transport  terms  are 


T\  —  (u  u  )  ,Xy  —  (u  v  )  ,XX 
T2  =  (u  V  ),yy—(v  v'),xy 

Tz  -  (uw  ),yz  -( Vfw'),XZ  • 

Following  the  procedures  outlined  in  section  8.4  for  taking  derivatives  of  the  data, 
all  the  terms  except  those  requiring  x  derivatives  were  estimated.  The  three  terms 
that  appeared  dominant  and  of  the  same  order  of  magnitude  were  C2,  D2,  and 
T2.  Note  that  in  T2,  (t>V),zy  could  not  be  estimated,  so  only  ( u'v'),yy  was  used 
to  estimate  T2.  Since,  from  the  arguments  made  by  Reynolds  and  Carr,  C\  should 
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also  be  significant,  this  suggests  that  a  simplified  model  for  transport  of  spanwise 
vorticity  might  be 

(u>2)  +  C\  +  C2  —  D2  +  T2.  (8.18) 

While  in  general,  data  was  taken  at  only  one  streamwise  position,  an  unsteady 
validation  profile  was  taken  at  a  location  upstream  of  the  primary  data  location,  but 
half-way  down  the  adverse  pressure  gradient  ramp  (see  figure  6.11).  By  differencing 
data  between  these  two  locations,  crude  estimates  of  C\  for  a  =  0°  could  be  made. 
These  estimates  indicate  that  for  tUo/l  <  4,  the  unsteady  term  tends  to  dominate 
the  left-hand  side  of  the  equation  8.18,  while  for  tU^/l  >  4,  C\  tends  to  dominate 
the  left-hand  side  of  the  equation. 

Estimates  of  Transport  Terms.  For  comparison,  figures  8.35  through  8.38 
show  the  time  evolution  of  Z?2,  figures  8.39  through  8.42  show  the  evolution  of 
T2,  and  figures  8.43  through  8.46  show  the  evolution  of  C2.  In  reviewing  these 
graphs,  one  should  remember  that  the  uncertainty  in  these  terms  requiring  several 
derivatives  of  the  data  is  on  the  order  of  100  percent.  However,  for  trends  between 
families  of  curves,  the  relative  uncertainty  between  them  is  probably  less.  Also,  in 
figure  8.43a,  C2  is  zero  by  the  arbitrary  decision  explained  in  chapter  7,  to  correct 
for  offset  errors  by  using  the  values  of  ( v )  and  (w)  at  Wq/1  =  0  and  a  =  0°  as 
zero  reference  values.  Finally,  some  of  the  trends  shown  by  the  spanwise  data  set  at 
z/D  =  0.125  seem  to  be  inconsistent  with  the  trends  of  the  other  four  spanwise  data 
sets.  This  was  the  last  data  set  taken,  and  it  was  taken  with  a  replacement  laser  after 
the  primary  laser  was  returned  to  the  manufacturer  for  warranty  repairs.  While  all 
procedures  were  followed  to  preclude  potential  inconsistencies,  some  trends  shown 
at  z/D  =  0.125  may  not  be  valid. 

In  making  observations  about  the  transport  terms,  it  will  be  useful  to  think  of 
terms  being  beneficial  when  they  make  the  left-hand  side  of  equation  8.18  more 
negative  near  the  surface,  since  this  indicates  an  increase  in  negative  vorticity  near 
the  surface.  The  major  observations  regarding  the  transport  of  spanwise  vorticity 
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are:  (l)  All  three  of  these  terms  are  large  only  in  a  region  that  is  very  near  the 
surface.  Note  that  for  this  reason  y/h  is  on  a  logarithmic  scale.  Also,  recall  that 
at  zero  pressure-gradient,  h+  =  230,  so  this  region  where  the  terms  are  significant 
corresponds  to  y+  <  20.  This  means  that  the  significant  transport  of  spanwise 
vorticity  is  occurring  primarily  in  the  viscous  layer  adjacent  to  the  surface.  (2)  In 
general,  diffusion  and  turbulent  transport  seem  to  be  the  two  terms  that  maintain 
the  vorticity  equilibrium  at  zero  pressure  gradient.  The  term  £>2  tends  to  increase 
the  negative  vorticity  near  the  surface  while  T2  tends  to  spread  out  the  vorticity 
away  from  the  surface.  (3)  At  a  =  18°,  the  diffusion  term  is  increased  in  the  down- 
wash  and  crossflow  regions.  This  effect  is  even  more  pronounced  once  the  adverse 
pressure  gradient  begins.  Figures  8.36a  through  8.38a  show  that  the  diffusion  term 
decreases  substantially  under  the  influence  of  the  adverse  pressure  gradient.  In  con¬ 
trast,  figures  8.36b  through  8.38b  show  that  the  diffusion  is  increasing  the  negative 
vorticity  in  the  downwash  and  crossflow  regions.  This  enhanced  viscous  diffusion  of 
vorticity  deep  within  the  boundary  layer  may  be  the  primary  mechanism  by  which 
the  vortex  generators  modify  the  transport  of  spanwise  vorticity.  (4)  While  T2  is 
positive  for  0.06  <  y/h  <  1.0  at  zero  pressure  gradient,  T2  becomes  nearly  zero 
or  negative  for  the  entire  inner  region,  in  response  to  the  adverse  pressure  gradi¬ 
ent.  This  indicates  that  while  T2  tends  to  spread  out  vorticity  in  opposition  to  the 
viscous  diffusion  during  zero  pressure  gradient,  T\ j  actually  assists  in  keeping  vor¬ 
ticity  near  the  surface  during  a  separating  flow.  While  the  vortex  generators  may 
modify  T2,  the  effect  is  quite  small.  (5)  Finally,  the  convection  term  is  insignificant 
at  a  =  0°  and  smaller  than  the  other  terms  at  q  =  18°.  Figures  8.43b  through 
8.46b  show  consistent  trends.  In  the  downwash  and  crossflow  regions,  and  negative 
value  of  C2  indicates  positive  vorticity  is  being  convected  into  the  surface.  This, 
therefore,  must  be  a  detrimental  secondary  effect  of  the  vortex  generators,  and  not 
the  primary  mechanism. 
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8.6  Transport  of  Momentum 

The  transport  equations  for  streamwise  and  spanwise  vorticity  have  been  exam¬ 
ined  to  increase  an  understanding  of  the  mechanisms  by  which  vortex  generators 
modify  the  behavior  of  a  separating  boundary  layer.  The  last  transport  equation 
that  will  be  examined  for  this  purpose  is  the  momentum  equation.  One  explana¬ 
tion  that  is  often  given  to  explain  how  vortex  generators  work  is  that  they  enhance 
mixing  and  energize  the  sluggish  flow  at  the  surface  with  energetic  fluid  from  the 
freestream.  In  terms  of  the  momentum  equation,  this  states  that  vortex  generators 
increase  the  convection  of  u  momentum  in  the  — y  direction.  The  results  that  follow 
will  demonstrate  that  this  is  true. 

The  Governing  Equation.  The  momentum  equation  was  stated  in  general 
form  as  equation  2.8.  When  applied  to  the  u  component  of  momentum  with  normal 
boundary  layer  assumptions,  it  can  be  expressed  as  follows: 

(u)  +  (uu)  ,x+(uv),y=  ^-(p),x+u(u),yy-(u'vl),y.  (8.19) 

The  intent  here  is  to  explore  only  the  second  convection  term,  (uv),y.  Often,  if 
continuity  has  been  subtracted  from  the  momentum  equation,  this  term  appears 
as  (v)(u), y.  Following  procedures  outlined  previously  for  taking  derivatives,  this 
convection  term  was  calculated  in  both  forms.  The  latter  form  tended  more  to 
map  out  the  velocity  gradient  than  the  actual  net  convection,  so  it  is  not  presented. 
The  first  representation  of  the  momentum  convection  is  presented  without  and  with 
vortex  generators,  for  four  different  phases,  in  figures  8.47  through  8.50.  To  further 
clarify  some  aspects  of  the  convection  term,  the  local  convection  of  momentum  was 
normalized  by  the  local  value  of  momentum.  The  analagous  information  is  shown 
in  this  form  in  figures  8.51  through  8.54. 

Observations.  Figures  8.47  through  8.54  do  verify  that  vortex  generators  help 
to  energize  the  inner  regions  of  the  boundary  layer,  and  to  suppress  the  loss  of 
momentum  from  the  outer  regions  of  the  boundary  layer.  Before  making  more 
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specific  observations,  it  is  worth  emphasizing  that  this  term  indicates  the  loss  of 
momentum  from  a  fluid  element  in  the  boundary  layer,  thus  a  positive  number 
represents  a  detrimental  effect  and  a  negative  number  represents  a  positive  effect. 
Also,  as  in  other  cases,  figures  8.47a  and  8.51a  indicate  zero  convection  because 
these  values  of  (v)  were  used  as  reference  values.  Now,  some  specific  observations 
are:  (1)  Figure  8.51b  shows  that  at  tUo/l  =  0  there  is  a  net  convection  of  momentum 
into  the  boundary  layer  in  the  downwash  and  crossflow  regions,  and  there  is  a  net 
convection  out  of  the  boundary  layer  in  the  upwash  regions.  The  only  trend  that 
is  at  all  surprising  is  the  small  region  of  momentum  loss  on  the  center  line.  This 
is  perhaps  due  to  a  stagnation-type  effect  as  the  downwash  parts  along  the  center 
line.  (2)  In  figures  8.52a  through  8.54a  there  is  loss  of  momentum  from  the  entire 
boundary  layer  which  is  quite  consistent  across  the  entire  span.  The  convection 
term  is  negative  very  near  the  surface  because  there  must  be  a  continual  flux  of 
momentum  into  the  wall  to  overcome  the  wall  shear  stress.  The  negative  convection 
term  appears  quite  large  because  it  is  normalized  by  the  very  small  amount  of 
momentum  that  remains  near  the  surface  as  separation  occurs.  (3)  Figure  8.52b 
shows  a  region  of  momuntum  flux  into  the  boundary  layer  at  zj D  —  0.125  and 
zfD  —  0.25.  This  region,  which  has  been  shaded  for  emphasis,  represents  the 
overlap  region  of  the  boundary  layer  (40  <  y+  <  200)  and  may  be  the  key  region 
that  is  energized  by  the  vortex  generator' .  Note  also  that  there  is  no  net  flux  on  the 
center  line,  where  the  downwash  is  not  nearly  as  strong  (reference  for  comparison 
figure  7.11).  Also,  as  expected,  there  is  a  flux  of  momentum  out  of  the  boundary 
layer  in  the  upwash  regions  of  the  vortex.  (4)  Figures  8.53b  and  8.54b  show  that  as 
the  adverse  pressure  gradient  continues,  there  is  momentum  loss  across  the  entire 
span,  but  the  average  value  is  considerably  less  at  0  <  zfD  <  0.25  than  the  average 
momentum  loss  shown  in  figures  8.53a  and  8.54a. 
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8.7  Changes  in  the  Structure  of  the  Turbulent  Boundary  Layer 

As  the  final  discussion  in  this  chapter,  this  section  will  address  how  the  struc¬ 
ture  of  the  turbulence  is  modified  by  the  combination  of  vortex  generators  and  an 
unsteady  adverse  pressure  gradient.  The  numerous  studies  of  streamwise  vortices 
embedded  in  a  steady  zero  pressure  gradient  turbulent  boundary  layer,  which  were 
summarized  in  section  3.1,  generally  all  indicated  that  the  turbulence  structure  was 
modified  significantly,  and  that  complex  models  will  be  required  to  predict  realistic 
values  of  the  Reynolds  stresses.  In  a  similar  manner,  the  unsteady  experiments  de¬ 
scribed  in  section  3.2  generally  found  significant  effects  of  the  unsteadiness  on  the 
turbulence  structure.  The  results  of  this  section  will  show  that  the  production  of 
turbulent  kinetic  energy  which  was  complicated  by  the  embedded  vorticity,  becomes 
even  more  confused  in  the  presence  of  the  unsteady  adverse  pressure  gradient. 

Chapter  7  already  presented  and  discussed  the  distribution  of  Reynolds  stresses 
in  the  boundary  layer,  and  how  they  are  modified  by  the  vortex  generators  and  an 
unsteady  adverse  pressure  gradient.  Due  to  the  cruder  data  grid  required  because  of 
the  long  sampling  times  associated  with  the  unsteady  cycle,  this  discussion  cannnot 
attempt  to  map  out  the  spanwise  variation  in  structural  parameters  in  a  manner 
analogous  to  Shabaka,  Mehta,  and  Bradshaw  (1985)  or  the  recent  work  of  Pauley 
and  Eaton  (1988).  This  section  will  address,  however,  the  nature  of  the  production 
of  turbulent  kinetic  energy  and  nature  of  the  anisotropy  of  the  Reynolds  stress 
tensor. 

Production  of  Turbulent  Kinetic  Energy.  Turbulent  kinetic  energy  is  de¬ 
fined  as  one-half  the  sum  cf  the  three  normal  Reynolds  stresses.  In  this  discussion 
k  will  denote  the  turbulent  kinetic  energy.  That  is 

k  =  *(<">  +  ^V'V^  +  =  \q2'  (8-20) 
Equation  2.9  defined  the  production  of  turbulent  kinetic  energy  as 

Pk  =  -(u'u'-Xu;,,). 


158 


When  this  term  is  summed  over  i  and  j,  there  are  nine  terms  that  contribute  to 
production: 

Pk  =  -(u'u')(u,x)  -  (uV)(u,y)  -  (u'v/)(u,z) 

-{u'v')(v,x)  -  {v'v')(v,y)  -  ( v'w'){v,z ) 

-(u'w')(w,x)  -  (v'w')(w,y)  -  (w'w'){w,z).  (8.21) 

Six  of  the  nine  terms  in  equation  8.21  could  be  estimated  by  taking  a  product  of 
the  appropriate  stress  and  ayorz  derivative  of  the  appropriate  velocity.  As  stated 
previously,  no  data  was  taken  that  would  provide  a  good  estimate  of  x  derivatives. 
In  addition,  due  to  higher  apparent  scatter  in  the  (w'w1)  measurements,  the  last 
term  was  also  excluded  from  the  estimate  of  turbulent  kinetic  energy  production. 
Estimates  of  production  based  on  the  five  remaining  terms  are  presented  in  figures 
8.55  through  8.57,  for  both  a  =  0°  and  a  =  18°,  for  six  discrete  points  in  time.  In 
the  production  calculation  the  second  term,  —(u'v')(u,y),  was  the  dominant  term 
over  most  of  the  region  in  most  of  the  cases.  This  is  consistent  with  the  observations 
of  Brereton  and  Reynolds  (1987)  for  unsteady  flows  without  vortex  generators. 

In  figures  8.55  through  8.57,  for  the  cases  where  a  =  0°,  the  magnitude  and 
location  of  most  of  the  production  are  consistent  with  the  observations  of  Brereton 
and  Reynolds.  For  comparative  production  terms  scaled  by  uT  and  y+,  the  appro¬ 
priate  scaling  factor  is  roughly  60.  This  means  that  the  peak  production  scaled 
by  t//ur  *s  between  0.15  and  0.2.  This  is  slightly  less  than,  but  of  the  same  rough 
magnitude  as  the  peak  production  value  of  0.25  seen  by  Laufer  (1954)  for  his  pipe 
flow  experiments.  Across  the  span,  as  the  adverse  pressure  gradient  evolves,  the 
peak  production  descreases  substantially.  After  tU$ //  =  2,  when  there  is  wall  region 
reverse  flow  over  most  of  the  span,  there  is  a  small  distribution  of  production  across 
the  shear  layer,  but  no  concentration  of  production  near  the  surface.  Since  the 
adverse  pressure  gradient  has  significantly  reduced  the  velocity  gradients,  or  strain 
rates,  that  are  responsible  for  production,  the  consequence  has  been  a  substantial 
drop  in  the  total  production  in  the  boundary  layer. 
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With  vortex  generators  on,  Pauley  and  Eaton  observed  “a  complex  interaction 
between  the  secondary  flow  field  induced  by  the  embedded  vortices  and  the  turbu¬ 
lence  structure  of  the  boundary  layer.”  The  same  would  appear  to  be  true  of  the 
Reynolds  stress  distributions  shown  in  figures  7.15  through  7.32  and  of  the  turbulent 
kinetic  energy  production  shown  in  figures  8.55  through  8.57.  Pauley  and  Eaton 
saw  a  high  concentration  of  turbulent  kinetic  energy  in  the  core  region  of  their  up¬ 
stream  case,  and  associated  this  with  the  velocity  deficit  in  the  core.  While  this 
experiment  did  not  show  much  velocity  deficit  in  the  vortex  core  (reference  figure 
8.5),  figure  7.15  does  indicate  a  large  concentration  of  ( u'u ')  in  the  core  region.  In 
contrast,  figure  7.18  indicates  there  is  not  a  large  concentration  of  ( v'v ')  in  the  core. 
With  respect  to  production  of  turbulent  kinetic  energy,  Pauley  and  Eaton  saw  an 
increase  in  the  vortex  core  region,  and  an  even  larger  increase  in  the  upwash  region. 
Figure  8.55  seems  to  indicate  an  increase  of  about  the  same  magnitude  in  both  the 
vortex  core  region  and  in  the  upwash  region.  After  the  adverse  pressure  gradient 
is  initiated,  production  decreases  in  the  downwash  region,  but  since  the  boundary 
layer  remains  attached  and  a  reasonable  strain  rate  is  maintained,  there  is  more 
production  here  than  in  the  separating  boundary  layer.  In  the  upwash  region,  there 
is  an  apparent  large  production,  despite  relatively  weak  velocity  gradients,  because 
of  the  large  magnitude  of  the  normal  stresses  in  the  chaotic  region  of  the  separated 
upwash. 

To  summarize  this  discussion  about  turbulent  kinetic  production,  vortex  genera¬ 
tors  definitely  alter  the  distribution  of  kinetic  energy  and  the  locations  of  maximum 
production.  In  an  unsteady  adverse  pressure  gradient,  the  modifications  observed 
for  steady  flow  are  greatly  amplified.  In  general,  turbulent  kinetic  energy  produc¬ 
tion  appears  to  respond  to  changes  that  occur  in  the  boundary  layer,  rather  than 
be  a  controlling  factor  in  the  changes  that  occur. 

Anisotropy  of  the  Reynolds  Stress  Tensor.  While  the  coarse  data  grid  and 
time  contraints  on  sampling  unsteady  events  made  a  detailed  investigation  of  the 
turbulence  structure  difficult  for  this  experiment,  it  was  possible  to  characterize 
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the  general  nature  of  the  turbulence.  The  general  nature  of  a  boundary  layer  is 
two-component  turbulence.  As  an  adverse  pressure  gradient  is  introduced,  does 
the  structure  remain  two-component,  or  does  it  tend  towards  one-component,  or 
isotropy?  This  section  will  briefly  address  this  question,  after  an  overview  of  essen¬ 
tial  definitions. 

First,  the  Reynolds  stress  anisotropy  tensor  is  a  useful  means  to  express  the 
variation  from  isotropic  turbulence.  It  is  defined  as 


Bn  -  Whi 


(8.22) 


In  general,  one  can  find  the  principle  axes  for  this  tensor,  along  which  only  the 
diagonal  terms  are  non-zero.  If  the  turbulence  is  isotropic,  then  all  the  terms  of 
the  tensor  are  zero.  The  manner  in  which  the  terms  of  6tJ  deviate  from  zero 
indicates  the  nature  of  the  turbulence  structure. 

A  compact  means  to  characterize  the  anisotropy  of  the  turbulence  structure  is  to 
map  b{j  onto  an  anisotropy  invariant  map.  While  the  foundation  for  this  technique 
lies  in  linear  algebra,  the  essence  of  the  technique  is  as  follows.  In  its  principle 
coordinates,  the  entire  tensor  can  be  expressed  in  terms  of  two  independent 

parameters,  which  will  be  named  II  and  III.  (II  and  III  are  the  invariants  of 
the  b{j  tensor.)  Therefore,  a  particular  state  of  turbulence,  characterized  by  a 
particular  set  of  values  of  the  6,-y  tensor,  corresponds  to  a  specific  location  on  the 
plane  formed  by  II  and  III.  Furthermore,  there  is  a  limited  region  or  space  on  this 
plane  that  is  physically  realizable.  The  limits  of  realizability  are  shown  in  figure  8.58. 
They  are  defined  by  isotropic  turbulence  at  the  origin,  a  straight  line  that  defines 
two-component  turbulence,  a  single  point  that  defines  one-component  turbulence, 
the  limiting  trajectory  from  isotropy  to  two-component  turbulence  (axisymmetric 
contraction),  and  the  limiting  trajectory  from  isotropy  to  one-component  turbulence 
(axisymmetric  expansion). 
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The  bij  tensor  was  computed  from  the  measured  values  of  R^j  at  the  center  line 
position,  z/D  =  0.  Of  the  various  spanwise  positions,  this  one  was  chosen  for  inves¬ 
tigation  because  it  is  a  region  where  the  vortex  generators  clearly  have  a  beneficial 
effect  in  preventing  separation.  While  the  (w'w1)  data  had  some  unrealistic  scatter, 
a  careful  review  of  the  data  at  z/d.  =  0  revealed  three  y/h  locations  in  the  middle 
portion  of  the  profile  where  something  had  apparently  malfunctioned.  After  these 
several  points  were  removed,  the  piecewise  least-squares  fit  of  ( w'w ')  looked  very 
reasonable.  The  resulting  values  of  6^  were  then  used  to  map  each  point  onto  the 
invariant  maps  shown  in  figures  8.59a  and  8.59b.  In  these  two  figures,  the  different 
symbols  indicate  different  points  in  time.  Two  observations  can  be  made  from  this 
figure.  (1)  The  state  of  turbulence  seems  to  tend  from  two-component  turbulence 
towards  isotropy.  To  see  this  tendency  more  clearly,  trajectories  at  individual  ylh 
locations  are  required.  (2)  There  is  no  major  difference  in  the  appearance  of  the 
two  maps  at  a  —  0°  and  a  =  18°. 

Figures  8.60  through  8.62  were  prepared  in  an  attempt  to  visualize  the  trajectories 
through  time  of  individual  points  in  the  boundary  layer.  To  keep  the  trajectories 
uncluttered,  a  representative  but  limited  number  of  y/h  positions  were  selected. 
Generally,  successive  arrowheads  correspond  to  the  dimensionless  times  tUg/l  =  0.4, 
tUo/l  =  2,  tUo/l  =  4,  and  tUo/l  =  8.  Occasionally  intermediate  arrowheads  were 
removed  when  they  obscured  some  aspect  of  the  trajectories.  Figure  8.60  shows 
the  innermost  region  of  the  boundary  layer  where  there  was  a  ve’-y  pronounced 
trajectory.  These  points  very  near  the  surface,  where  the  turbulence  is  originally 
two-component,  move  very  strongly  to  the  left  in  the  general  direction  of  isotropy. 
The  same  trend,  but  less  dramatic,  is  seen  for  the  inner  region  with  a  —  18°.  The 
trend  may  be  less  dramatic  because  the  vortex  generators  bring  some  more  isotropic 
turbulent  eddies  closer  to  the  surface.  Figure  8.61  shows  trajectories  for  points  in 
the  outer  portion  of  the  boundary  layer.  These  points  originate  near  isotropy  and 
tend  to  be  driven  to  the  left,  but  not  nearly  as  strongly.  Finally,  figure  8.62  shows 
points  on  the  outer  edge  of  the  boundary  layer.  At  a  =  0°  these  points  start 
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near  isotropy  and  remain  there,  but  at  a  =  18°,  some  of  the  points  tend  to  be  less 
isotropic  and  their  trajectories  are  more  chaotic.  This  may  be  because  the  enhanced 
mixing  from  the  vortex  generators  brings  out  some  packets  of  fluid  from  the  surface 
that  are  less  isotropic. 

To  conclude,  the  invariant  maps  indicate  several  things  about  the  turbulent  struc¬ 
ture.  First,  at  this  spanwise  location  where  the  vortex  generators  clearly  delay  sep¬ 
aration,  there  is  no  major  change  in  the  anisotropy  of  the  turbulent  structure  as  a 
result  of  the  embedded  vortices.  Second,  under  the  influence  of  a  strong  adverse 
pressure  gradient,  points  near  the  surface  are  driven  from  a  state  of  two-component 
turbulence  to  a  state  much  nearer  isotropy. 

8.8  Closure 

This  chapter  has  used  many  different  analytical  tools  to  investigate  the  nature  of 
embedded  streamwise  vortices  in  an  unsteady  adverse  pressure  gradient.  By  far,  the 
tool  that  yielded  the  most  information  was  the  time  response  analysis  of  secion  8.2. 
The  results  presented  in  figures  8.11  through  8.16  should  be  of  use  to  engineers  in 
the  design  of  flow  control  systems.  In  addition,  these  results  showed  that  the  vortex 
generators  do  not  modify  the  immediate  response  of  the  boundary  layer,  but  they 
do  modify  the  magnitude  of  the  slower  convective  response.  The  results  presented 
in  the  other  sections  are  useful  for  comparisons  to  other  studies  in  steady  flow,  and 
they  add  additional  insight  into  the  basic  conclusions  drawn  in  section  8.2. 
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Figure  8.1a  Steady  mean  profiles  in  wall  coordinates  for  a  =  0° 


Figure  8.1b  Steady  mean  profiles  in  wall  coordinates  for  a  =  18° 
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Figure  8.2a  Spanwise  variation  in  Cf  for  steady  flow 
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Figure  8.2b  Spanwise  variation  in  ur  for  steady  flow 
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Figure  8.6b  Contou 
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Figure  8.7a  Contours  of  {u)/Uq  at  Wq/1  =  4  for  a  —  0° 


Figure  8.7b  Contours  of  ( u)/Uq  at  tUg/l  =  4  for  a  —  18° 
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Figure  8.8b  Conto 
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Figure  8.9b  Contours  of  (u)/Uq  at 
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re  8.10a  Time  response  of  integral  parameters  at  z/D  =  0 


Figure  8.10b  Time  response  of  shape  factor  at  z/D  =  0 
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Figure  8.11  Mathematical  model  lor  time  response  of  velocity 


tUg/l,  a-C 

Figure  8.12a  History  of  $  for  a  =  0°  at  z/D  =  0 
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Figure  8.12b  History  of  $  for  a  =  18°  at  z/D  =  0 
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Figure  8.14d  Initial  delay  in  response  of  (u) 
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Figure  8.14e  Duration  of  initial  response  of  (u) 
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Figure  8.14f  Convective  time  constant  for  (u) 
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Figure  8.16a  Immediate  response  rate  of  (u'u1) 
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Figure  8.1Cb  Initial  delay  in  response  of  (u,u/) 


Skin  Friction  Coefficient,  C 


0  0.1  0.2  0.3  0.4  0.5 

z/D  (a=0) 

Figure  8.17a  Spanwise  variation  in  Cf  for  a  =  0° 


Figure  8.17b  Spanwise  variation  in  Cf  for  a  =  18° 
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Figure  8.18a  Spanwise  variation  in  ur  for  a  =  0° 
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Figure  8.19a  (u)  profiles  in  wall  coordinates  for  a  =  0°  at  tUg/l  =  0 


Figure  8.19b  (a)  profiles  in  wall  coordinates  for  a  =  18°  at  tUg/l  =  0 
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Figure  8.27b  Streamwise  vorticity,  {ux)h/Un,  at  tUo/l  =  0.4 
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Figure  8.30a  Streamwise  vorticity,  (wx)h/Uo,  at  tU0/l 
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Figure  8.30b  Streamwise  vorticity,  (ux)h/Uo,  at  tU^jl  =  8 
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Figure  8.32  History  of  streamwise  circulation 
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Figure  8.36a  Diffusion  of  spanwise  vorticity  for  a  =  0°  at  UJq/1  =  2 
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Figure  8.36b  Diffusion  of  spanwise  vorticity  for  a  =  18°  at  Wq/1  —  2 
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■e  8.37a  Diffusion  of  spanwise  vorticity  for  a  =  0°  at  tU^/l  =  4 


Figure  8.37b  Diffusion  of  spanwise  vorticity  for  a  =  18°  at  tUo/l  =  4 


202 


1  /Reh<a2coz/3y2>  (h3/Ua 


(aVv-ya^Kh/u/  f  o*(uW9y2)(wu0) 


y/h  (a=0,  tU0/l=0) 

e  8.39a  Turbulent  transport  of  spanwise  vorticity  for  a  =  0°  at  Wq/I  =  0 
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Figure  8.39b  Turbulent  transport  of  sp&nwise  vorticity  for  a  =  18°  at  Wq/I  =  0 
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Figure  8.40a  Turbulent  transport  of  spanwise  vorticity  for  a  =  0°  at  tU0/l  =  2 
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Figure  8.40b  Turbulent  transport  of  spanwise  vorticity  for  a  =  18°  at  Wq/I  =  2 
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Figure  8.41b  Turbulent  transport  of  spanwise  vorticity  for  a  =  18°  at  tU^/l  —  4 
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Figure  8.42a  Turbulent  transport  of  spanwise  vorticity  for  a  =  0°  at  tUo/l  = 
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Figure  8.42b  Turbulent  transport  of  spanwise  vorticity  for  a  =  18°  at  tUo/l 
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re  8.43a  Convection  of  spanwise  vorticity  by  ( v }  for  a  =  0°  at  tU$ [l  =  0 
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Figure  8.43b  Convection  of  spanwise  vorticity  by  (v)  for  a  =  18°  at  tUo/l  =  0 
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ire  8.44a  Convection  of  spanwise  vorticity  by  ( v )  for  a  =  0°  at  tUo/l  =  2 


Figure  8.44b  Convection  of  spanwise  vorticity  by  { v )  for  a  —  18°  at  tU$jl  =  2 
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ure  8.46a  Convection  of  spanwise  vorticity  by  (v)  for  a  =  0°  at  tUg/l  =  8 
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Figure  8.46b  Convection  of  spanwise  vorticity  by  (r)  for  a  =  18°  at  tUo/l  =  8 
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Figure  8.47a  Convection  of  momentum  by  ( v )  for  a  —  0°  at  tUo/l  =  0 
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Figure  8.47b  Convection  of  momentum  by  ( v )  for  a  =  18°  at  Wq/1  =  0 


212 


213 


214 


3<uv)/3y(h/U 


y/h  (a=0,  tU0/l=8) 

Figure  8.50a  Convection  of  momentum  by  (v)  for  a  =  0°  at  tUo/l  =  8 
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Figure  8.50b  Convection  of  momentum  by  (t;)  for  a  =  18°  at  tUo/l  =  8 
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Figure  8.53a  Normalized  convection  of  momentum  by  ( v )  for  a  =  0° 
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Figure  8.55  Turbulent  kinetic  energy  production  at  tUo/l  =  0  and  tUo/l  —  0.4 
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Figure  8.58  Limits  of  realizability  on  an  invariant  map  of  6, 
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Figure  8.61a  Trajectories  of  6ty  for  outer  boundary  layer  with  a  =  0° 
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Figure  8.61b  Trajectories  of  6^  for  outer  boundary  layer  with  a  =  18° 
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The  principle  is  most  important,  not  the  detail. 
—  Theodore  von  Karman,  1954 
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9.  Conclusions  and  Recommendations  for  Future  Work 


9.1  Conclusions 

This  program  has  fulfilled  the  objectives  which  were  identified  in  chapter  1. Chap¬ 
ters  4  and  5  described  the  facility  and  the  three-component  LDA  s>stem  that  was 
developed  and  tested  specifically  for  use  in  this  program.  Chapter  7  documented  in 
detail  the  profiles  of  each  velocity  component  and  all  six  components  of  the  Reynolds 
stress  tensor.  Finally,  chapter  8  described  the  time  response  of  the  boundary  layer 
and  provided  some  insight  into  the  mechanisms  that  make  vortex  generators  work. 
In  particular,  figures  8.11  through  8.16  summarize  the  most  significant  observations. 

Major  Conclusions.  The  major  conclusions  resulting  from  this  program  are  as 
follows: 

1.  An  array  of  weak  counter-rotating  vortices  markedly  alters  the  nature  of  the 
boundary  layer,  both  in  steady,  zero  pressure  gradient  flows  and  in  unsteady, 
adverse  pressure  gradient  flows. 

2.  In  steady,  zero  pressure  gradient  flows,  this  experiment  confirmed  the  findings  of 
earlier  experiments.  Key  effects  of  the  vortex  generators  on  the  steady  boundary 
layer  are: 

a.  A  decrease  in  the  interval  thicknesses  in  the  downwash  region  and  an  increase 
in  the  upwash  region. 

b.  An  increase  in  Cj  over  most  of  the  span,  but  a  decrease  in  Cj  in  the  upwash 
region  of  the  vortex  pair. 

c.  a  substantial  modification  of  the  wake  region  with  only  minor  distortion  of  the 
‘inner  region’  of  the  nominally  2-D  turbulent  boundary  layer. 

3.  For  unsteady  flow,  the  time  response  of  the  freestream  velocity  demonstrates  that 
convection  is  the  mechanism  by  which  vortex  generators  modify  the  response  of 
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the  boundary  layer.  Modeling  the  response  with  two  time  scales  proved  to  be  a 
useful  technique,  and  led  to  the  following  observations: 

a.  Vortex  generators  do  not  affect  the  basic  response  rates  of  the  boundary  layer, 
but  they  do  affect  the  magnitude  of  the  convective  portion  of  the  response. 

b.  There  is  an  initial  response  throughout  the  boundary  layer  which  is  unaffected 
by  the  presence  of  vortex  generators.  The  duration  of  this  initial  response  is 
the  same  throughout  the  boundary  layer  and  appears  to  reflect  the  control 
valve  response  time,  but  the  rate  and  magnitude  of  the  initial  change  are  less 
in  the  inner  region  of  the  boundary  layer. 

c.  There  is  also  a  slow  or  convective  response  throughout  the  boundary  layer 
which  is  substantially  modified  by  the  presence  of  the  vortex  generators. 

d.  Because  the  magnitude  of  the  convective  response  is  substantially  reduced  in 
the  downwash  region,  a  far  greater  percentage  of  its  response  occurs  during 
the  initial  response.  The  apparent  effect  of  vortex  generators  in  the  downwash 
region,  therefore,  is  that  they  substantially  decrease  the  time  required  for  the 
boundary  layer  to  adjust  to  a  step  change  in  the  pressure  gradient. 

4.  The  response  of  {v)  appears  to  be  governed  by  the  same  mechanisms  that  govern 
the  initial  response  of  (u),  and  in  the  outer  region  of  the  boundary  layer,  the 
response  rates  are  the  same  for  (u)  and  (v).  However,  near  the  surface,  the 
magnitude  and  rate  of  the  response  of  (v)  approach  zero  in  the  same  manner  as 
(v)  approaches  zero. 

5.  The  turbulence  in  the  inner  region  of  the  boundary  layer  begins  to  respond  at  the 
same  time  as  the  initial  response  of  (u) ,  but  responds  at  a  much  slower  rate.  In 
contrast,  the  turbulence  in  the  outer  portion  of  the  boundary  layer  responds  at 
the  same  rate  as  the  growth  of  the  boundary  layer,  which  occurs  on  the  convective 
time  scales  of  the  freestream  velocity. 

Additional  Observations.  The  following  are  additional  observations  that  add 
insight  into  the  behavior  of  vortex  generators: 
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1.  Under  unsteady  conditions,  the  displacement  thickness  increased  by  only  a  factor 
of  two  in  the  downwash  of  the  vortex  generators,  compared  to  a  five-fold  increase 
when  the  vortex  generators  were  at  an  angle  of  attack  of  0°. 

2.  Under  unsteady  conditions,  the  strength  of  the  the  embedded  streamwise  vortex, 
as  characterized  by  the  circulation,  holds  steady  for  about  two  convective  time 
periods  then  begins  to  decrease.  After  about  three  convective  time  periods,  the 
vorticity  from  the  embedded  vortex  loses  its  identity,  and  the  flow  appears  to 
become  dominated  by  streamwise  vortices  of  the  opposite  sign. 

3.  Plots  of  the  spanwise  vorticity  reveal  the  nature  of  a  boundary  layer  responding 
to  an  adverse  pressure  gradient.  With  a  zero  pressure  gradient,  the  bound¬ 
ary  layer  has  a  high  concentration  of  negative  spanwise  vorticity  very  near  the 
wall.  After  the  adverse  pressure  gradient  is  initiated,  a  portion  of  this  negative 
vorticity  diffuses  into  the  wall  due  to  the  pressure  gradient,  and  the  remainder 
redistributes  throughout  the  boundary  layer.  The  vortex  generator  downwash 
preserves  a  higher  concentration  of  negative  vorticity  near  the  wall.  This  results 
in  a  healthy  velocity  profile  and  precludes  the  reverse  flow  that  would  otherwise 
occur.  Estimates  of  the  terms  for  transport  of  spanwise  vorticity  indicate  that 
the  majority  of  spanwise  vorticity  transport  occurs  in  the  viscous  region  of  the 
boundary  layer,  where  large  vorticity  gradients  exist,  and  that  the  viscous  and 
turbulent  diffusion  terms  dominate, 

4.  Momentum  transport  by  { v )  is  a  significant  mechanism  by  which  vortex  gener¬ 
ators  modify  the  nature  of  a  separating  boundary  layer.  During  zero  pressure 
gradient  conditions,  the  downwash  convects  momentum  towards  the  surface.  It 
continues  to  carry  momentum  into  the  region  40  <  y+  <  200  for  a  brief  period 
after  the  adverse  pressure  gradient  is  initiated.  As  the  pressure  gradient  contin¬ 
ues,  (v)  begins  to  convect  momentum  away  from  the  entire  boundary  layer,  but 
the  net  efflux  is  significantly  less  in  the  downwash  regions  of  the  vortex  generator 
than  for  the  case  with  vortex  generators  at  an  angle  of  attack  of  0°. 
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5.  The  regions  of  high  turbulent  kinetic  energy  are  modified  by  the  presence  of  the 
streamwise  vortices  in  the  manners  observed  by  previous  steady  experiments. 
Under  an  adverse  pressure  gradient,  the  turbulence  is  amplified  and  moves  out¬ 
ward  with  the  separating  boundary  layer. 

9.2  Recommendations  for  Future  Work 

This  program  has  left  several  questions  unanswered  which  should  be  pursued  in 
future  programs.  It  has  also  identified  some  of  the  strengths  and  limitations  of 
this  facility  for  experiments  like  this.  In  particular,  the  strengths  of  the  facility  are 
its  computerized  flow  control  and  data  acquisition  which  make  it  ideal  for  varying 
many  parameters  and  flow  configurations.  Water  is  also  ideal  for  flow  visualization. 
In  contrast,  unsteady  experiments  which  require  mapping  the  flow  field  in  two  or 
three  dimensions,  are  less  ideal  because  of  the  very  long  sampling  times  required  to 
obtain  good  converged  turbulence  statistics.  The  recommendations  that  follow  are 
based  on  these  considerations. 

Flow  Visualization.  Some  very  strong  spanwise  secondary  flow  develops  during 
the  adverse  pressure  gradient,  particularly  in  the  upwash  region  (see  figures  7.13 
and  7.14).  Also,  the  contour  plots  of  streamwise  vorticity  indicate  that  vortices 
of  opposite  sign  are  formed  (see  figures  8.27  through  8.31),  perhaps  because  of 
the  spanwise  gradient  of  the  streamwise  velocity.  Flow  visualization  would  clarify 
whether  this  is  a  significant  flow  phenomenon  or  a  consequence  of  irregularities  in 
the  facility  which  are  amplified  by  the  adverse  pressure  gradient.  In  either  case, 
this  is  an  important  question  that  should  be  resolved.  A  new  test  plate  should  be 
constructed  that  allows  the  introduction  of  ‘tufts’  of  dye  directly  into  the  regions 
where  the  large  spanwise  flows  are  occuring. 

Parametric  Study  of  Active  Control.  The  time  response  information  from 
this  study  provides  information  that  should  be  useful  in  the  design  of  an  active 
flow  control  system.  However,  a  control  engineer  would  actually  like  to  know  how 
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these  response  parameters  vary  with  the  configuration  and  dimensions  of  the  vor¬ 
tex  generators,  and  how  these  parameters  change  when  the  vortex  generators  are 
moved  in  conjunction  with  the  unsteady  pressure  gradient.  This  would  require  two 
additions  or  modifications  to  the  existing  test  surface:  (a)  fabricate  test  surfaces 
for  several  alternate  configurations,  including  larger  blades  and  co-rotating  pairs; 
and  (b)  replace  the  belt-drive  mechanism  with  a  bevel  gear  mechanism  which  is 
driven  by  one  of  the  valve-control  motors.  The  data  acquisition  should  emphasize 
time  response  of  the  freestream  velocity,  with  very  little  or  no  emphasis  on  three- 
component  LDA  measurements.  While  phase  averaging  over  300-500  cycles  would 
still  be  recommended  for  good  convergence  of  the  velocity  data,  from  hindsight,  a 
4-second  cycle  would  be  adequate  for  the  time  response  information  presented  in 
this  program.  This  type  of  study  would  provide  a  wealth  of  information  for  design 
of  active  flow  control  system. 

9.3  Closure 

This  program  was  intended  to  be  a  first  look  at  embedded  streamwise  vorticity 
in  unsteady  flows.  Since  it  was  a  first  step,  the  approach  included  some  very  sophis¬ 
ticated  techniques,  and  some  rather  simple  techniques.  Because  of  the  complexity 
of  a  three-dimensional  unsteady  flow,  the  simple  time  response  analysis  actually 
provided  the  most  useful  information.  If  more  three-dimensional  LDA  measure¬ 
ments  jf  turbulent  quantities  are  desired  in  a  flow  such  as  this,  then  quality  data 
will  require  many  (several  thousand)  cycles  per  phase  average  on  a  denser  spatial 
data  grid  than  used  in  this  program.  The  most  useful  information,  however,  will 
probably  come  from  the  simpler  techniques. 
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Appendix  A:  Uncertainty  Analysis 


A.l  Measured  Quantities 

Methodology.  This  analysis  is  based  on  two  key  ingredients  of  an  uncertainty 
analysis,  as  described  by  Kline  and  McClintock  (1953).  They  are  (1)  a  clearly  stated 
confidence  level  for  which  the  uncertainty  interval  applies,  and  (2)  an  assumption  of 
linearly  independent  factors  that  contribute  to  the  uncertainty.  With  respect  to  the 
first  ingredient,  all  the  uncertainty  ranges  estimated  in  this  report  will  be  stated  at 
20:1  odds,  which  is  equivalent  to  a  95%  confidence  level.  With  respect  to  the  second 
item,  for  n  linearly  independent  sources  of  uncertainty,  the  total  uncertainty  can  be 
visualized  as  the  length  of  a  vector  in  n  dimensional  space  formed  by  n  orthogonal 
vectors  placed  head  to  tail.  The  length  of  the  total  uncertainty  vector  is  then  just 
the  square-root  of  the  sum  of  the  squares  of  the  individual  uncertainties. 

Moffat  (1988)  simplifies  the  categorization  of  uncertainties  as  bias  (fixed)  and 
random.  He  points  out  that  some  of  the  fixed  uncertainties  actually  result  from 
slowly  varying  effects,  but  for  simplicity  here  such  uncertainties  will  be  considered 
with  the  bias.  The  analysis  that  follows  addresses  uncertainties  in  the  three  phase- 
averaged  velocities,  (u),  ( v ),  and  ( w ),  the  three  normal  stress  terms,  ( u'u '),  ( v'v ’), 
and  (u/u/),  and  the  three  shear  stress  terms  ( u'v '),  ( u'w '),  and  (t/u/).  For  each 
mean  velocity,  the  uncertainty  will  be  calculated  for  two  cases,  a  representative 
peak  value  and  a  representative  low  but  significant  plateau  value.  These  values 
are  presented  in  table  A.l.  This  table  also  presents  the  cumulative  results  of  the 
discussions  that  follow. 

Bias  Uncertainties.  The  largest  fixed  uncertainties  in  LDA  measurements 
of  velocity  in  this  facility  are  (1)  those  in  the  calibration  resulting  from  uncer¬ 
tainty  in  the  beam  intersection  angle,  (2)  background  electronic  noise,  (3)  cross- 
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contamination  of  the  three  velocity  component  signals,  and  (4)  frequency-voltage 
conversion  errors  in  the  downmixing  and  tracker/counter  instrumentation. 

Calibration  Error.  One  significant  potential  source  of  bias  in  LDA  measurements 
in  this  facility  comes  from  determination  of  the  beam  intersection  angle.  For  the 
w  calibration,  holes  are  burned  in  construction  paper  at  several  heights  above  the 
surface.  While  the  distances  are  small,  the  height  change  is  known  quite  accu¬ 
rately  from  the  stepping  motor  controller,  and  the  distance  between  burn  spots  is 
determined  on  a  micro-positioning  platform.  Determination  of  the  angle  for  u  and 
v  involves  projecting  the  beams  about  2m  across  the  lab  onto  a  wall.  Here  the 
distances  are  greater  but  the  true  center  of  each  beam  is  harder  to  establish.  The 
calibration  constant  is  proportional  to  the  wavelength  of  laser  light  and  inversely 
proportional  to  the  sine  of  the  half-angle  of  intersection.  For  u  and  v,  this  angle  is 
on  the  order  of  2.8°,  while  it  is  about  7.6°  for  w. 


The  best  indictors  of  repeatibility  and  potential  bias  errors  were  obtained  during 
the  laboratory  validation  phase.  When  the  two  independent  LDA  systems  were 
operated  side  by  side,  repeated  tests  indicated  that  the  two  systems  agreed  with 
each  other  to  within  about  1%  of  the  freestream  velocity.  Therefore,  0.01  was  used 
as  the  relative  calibration  bias  estimator  for  u  and  v.  Because  of  the  larger  crossing 
angle  for  w,  its  sensitivity  to  the  same  amount  of  error  is  less,  about  40%  of  u, 
so  0.004  was  used  as  the  relative  calibration  bias  estimator  for  w.  An  error  in  the 
calibration  slope  affects  the  fluctuating  quantities  in  a  similar  manner,  except  that 
each  must  be  treated  as  a  product  term.  That  is 


A  (u  u)  .  Au' 

— L- =  2u  - . 

Aslope  Aslope 


(A.l) 


Electronic  Background  Noise.  Section  5.1  described  the  facility  validation  exper¬ 
iments  to  differentiate  between  electronic  background  noise  and  the  true  freestream 
turbulence  of  the  facility.  The  conclusions  of  that  study  were  that  the  down-mixing 
circuitry  added  an  apparent  turbulence  of  about  1%  of  the  0.7m/ sec  freestream 
velocity.  In  terms  of  the  variance  and  the  parameters  of  this  experiment,  this  is 
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equivalent  to  0.0001t/g.  This  value  was  used  for  terms  containing  u'  and  v1.  Note 
that  the  electronic  noise  is  subtracted  from  the  mean  values  and  therefore  does  not 
affect  them. 

Cross  Contamination.  In  a  normal  turbulent  boundary  layer,  u  is  about  30  times 
larger  than  v.  This  means  that  a  small  contamination  of  u  can  totally  obliterate  v. 
Fortunately,  during  the  adverse  pressure  gradient,  the  v  component  grows  signifi¬ 
cantly  while  u  decreases,  so  the  contamination  also  decreases.  While  the  alignment 
procedure  was  geared  toward  minimizing  the  cross-contamination,  the  optics  can¬ 
not  be  rotated  much  better  than  about  0.5°.  This  means  there  is  a  potential  u 
contamination  in  v  of  about  O.Olu,  which  shows  up  as  a  significant  bias  in  small 
values  of  v.  The  same  potential  for  contamination  exists  with  the  w  measurements, 
so  the  same  estimate  for  bias  was  used  for  w. 

Remaining  Bias  Uncertainties.  The  remaining  bias  uncertainties  that  have  not 
specifically  been  discussed  yet  are  those  due  to  biases  in  the  instrumentation,  such 
as  the  repeatability  of  the  Bragg  cell  frequency,  the  ability  of  trackers  and  counters 
to  detect  the  doppler  frequency,  the  analog  to  digital  conversion  process,  and  the 
voltage  calibration  slopes  and  offsets.  The  manufacturer  of  the  LDA  equipment 
gives  an  overall  system  accuracy  of  0.4%  of  full  scale  range,  which  corresponds  to 
0.02  uncertainty  on  u.  Calibration  data  indicates  that  the  uncertainty  in  the  analog 
to  digital  conversion  process  is  on  the  order  of  0.05%.  While  these  factors  were  not 
discussed  until  now,  they  were  implicitly  included  in  the  calibration  error  considered 
above.  That  was  a  system  to  system  comparison,  where  the  indication  was  that  the 
uncertainty  was  typically  1%  instead  of  the  2%  mentioned  above.  In  addition,  these 
bias  errors  will  be  shown  to  be  small  in  comparison  to  the  random  uncertainty  that 
will  be  considered  next. 

Random  Uncertainties.  Analysis  of  the  random  uncertainty  is  based  on  an 
assumed  statistical  distribution,  based  on  past  experience,  of  the  variable  being 
measured.  The  phase-averaged  mean  values,  the  normal  stresses,  and  the  shear 
stresses  each  require  separate  treatment. 
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Mean  Values.  The  velocity  usually  has  a  normal  distribution,  therefore  Student’s 
t  distribution  describes  how  to  estimate  the  uncertainty  from  the  measured  estimate 
of  the  variance.  In  this  case,  for  phase-averaged  values  (u),  ( v ),  and  (tu),  the 
corresponding  measured  variance  estimates  are  simply  (u'u1),  (v'v'),  and  ( w'w '). 
At  20  :  1  odds,  the  uncertainty  interval  is 


±(„_!(.975)-^.  (-4-2) 

In  equation  A.2,  s  is  the  square-root  of  the  measured  variance  and  lor  phase-averages 
performed  over  500  cycles,  n  =  500.  The  t  distribution  has  a  value  of  1.98  for  n  = 
120,  then  it  approaches  1.96  as  n  approaches  infinity.  These  uncertainty  calculations 
used  the  value  1.98.  In  addition,  the  peak  normal  stress  values  were  not  used  for 
this  estimate  because  the  peak  usually  applies  over  a  very  small  portion  of  the  data 
space  and  projects  an  unrealistically  large  uncertainty  for  most  of  the  data.  After 
viewing  the  data  presented  in  chapter  7,  a  more  realistic  value  for  the  normalized 
variance  estimate  appeared  to  0.005  for  all  three  normal  stresses.  Since  this  is  about 
one-third  the  peak  value  of  (u'u1),  there  are  regions  where  the  uncertainty  is  larger 
than  that  stated  by  this  analysis. 

Normal  Stresses.  Consistent  with  the  assumption  above  that  the  variance  of  the 
velocity  has  a  normal  distribution,  it  follows  that  the  variance  has  a  x2  distribution. 
Under  these  conditions,  Rice  (1988)  defines  the  confidence  interval  of  the  variance 
as 


no ‘ 


*n-l(a/2) 


, 


no * 


(4.3) 


where  d2  in  this  case  is  the  estimate  of  the  variance,  (u'u'),  ( v'v '),  or  (w'w').  Since 
X2  is  not  a  symmetrical  distribution,  the  two  values,  X2g9(-975)  and  x2gg(-025)  must 
each  be  calculated.  For  large  values  of  n  for  which  x2  is  not  tabulated,  Rice  gives 
a  useful  approximation 

Xp  -  \(zP  ±  -  l)2  (A.4) 

In  equation  A.4,  zp  is  1.96  for  the  standard  20:1  odds  used  in  this  analysis. 
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Shear  Stresses.  The  confidence  interval  for  the  shear  stress  is  a  little  more  com¬ 
plicated.  While  one  can  treat  the  cross-product  terms  as  the  product  of  two  x2 
distributed  variables,  a  more  correct  approach  that  yields  somewhat  tighter  inter¬ 
vals  is  given  by  Neter,  Wasserman,  and  Kutner  (1985).  They  sight  a  transformation 
by  R.  A.  Fisher  of  the  correlation  coefficient,  ri2, 


Unlike  the  shear  stresses,  this  transformed  variable  is  approximately  normal,  with 
a  variance 

»v> = -H-  i-4-6) 

n  —  6 

The  confidence  interval  for  z'  at  the  95%  confidence  level  is  therefore 

±1.96o,(2/).  (A. 7) 

For  estimation  purposes,  since  the  u'v1  correlation  coefficient  is  on  the  order  of  0.4 
across  half  of  the  boundary  layer,  and  since  the  correlation  coefficients  of  u'w'  and 
v'w'  approach  this  value  in  the  regions  where  vortex  generator  action  is  significant, 
0.4  is  a  representative  value  over  most  of  the  regions  where  the  shear  stress  is 
significant.  This  value  was  transformed  into  a  value  of  z'  to  get  a  confidence  interval 
in  z' ,  and  then  transformed  back  to  obtain  a  confidence  interval  of  the  correlation 
coefficient.  Normalizing  this  interval  by  the  original  value  of  0.4  yielded  a  relative 
uncertainty  estimate  for  the  shear  stress  terms  of  0.36. 

Summary.  Table  A.l  presents  the  significant  numbers  used  in  and  obtained 
from  this  uncertainty  analysis.  Figure  A.l  helps  to  visualize  the  total  relative  un¬ 
certainties  in  bar  graph  form.  Figure  A.l  and  table  A.l  indicate  that  the  smallest 
uncertainties  occur  in  the  larger  values  of  the  mean  velocities  but  the  largest  uncer¬ 
tainties  occur  in  the  very  small  values  of  v  and  w.  This  uncertainty  is  due  mostly 
to  the  potential  contamination  by  u.  Fortunately,  these  small  values  mostly  oc¬ 
cur  under  zero  pressure-gradient  conditions  without  vortex  generators,  where  not 
much  interesting  is  happening  anyway.  Errors  in  the  normal  and  shear  stresses  were 
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about  the  same  for  small  and  large  values,  hence  only  a  single  representative  value 
is  shown.  With  the  exception  of  small  values  of  v  and  w,  the  dominant  uncertainty 
is  random  in  nature. 

A. 2  Calculated  Quantities 

In  general,  the  goal  of  this  program  was  to  estimate  the  general  magnitude  of 
various  terms  that  describe  the  transport  of  momentum  and  vorticity.  In  contrast 
to  many  engineering  studies  where  regressions  yield  engineering  coefficients  with 
a  small  relative  uncertainty,  this  study  required  multiple  derivatives  of  the  data, 
a  dangerous  process  which  results  in  large  uncertainty  intervals.  This  section  will 
attempt  to  place  some  uncertainty  intervals  on  the  computations  presented  in  chap¬ 
ter  8.  While  the  philosophy  regarding  20:1  odds  will  carry  over  from  section  A.l, 
the  attempt  at  statistical  rigor  will  not.  The  following  paragraphs  summarize  the 
general  assumptions  and  judgements  that  went  into  the  estimates  of  computational 
uncertainty. 

The  determination  of  uT  and  Cj  were  based  on  curve-fits  to  the  (u)  profiles. 
Curve-fits  generally  have  a  lower  uncertainty  than  that  of  the  individual  points 
on  the  curve.  As  stated  in  chapter  8,  comparison  of  the  different  methods  for 
determination  of  ur  indicated  a  repeatability  to  within  about  2%  to  3%.  This  is 
the  basis  for  the  uncertainty  estimate  for  ur.  Since  Cj  is  computed  from  u^,  its 
uncertainty  is  double  that  of  uT. 

All  the  remaining  computations  involved  derivatives  of  the  data.  In  all  cases, 
a  piece-wise  least-squares  parabolic  fit  of  the  data  was  obtained,  and  then  a  cubic 
spline  was  used  to  estimate  the  slope.  The  curve  fit  should  generally  be  better  than 
the  uncertainty  in  the  individual  points  on  the  curve,  but  because  of  the  variety 
of  curves  fit,  no  improvement  was  assumed  here.  To  test  the  uncertainty  in  the 
spline  determination  of  derivatives,  the  spline  was  tested  against  a  rather  nasty 
non-polynomial  curve  representative  of  the  shear  stress  profiles.  The  spline  showed 
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a  worst-case  deviation  of  about  10%  from  the  analytical  slope,  so  10%  was  added 
to  the  uncertainty  of  each  velocity  first-derivative.  For  subsequent  derivatives  of 
the  velocity  and  for  both  derivatives  of  (uV),  various  manipulations  of  the  data 
seemed  to  indicate  that  it  was  more  appropriate  to  double  the  uncertainty  with 
each  operation.  Finally,  computational  products  were  treated  according  to  the 
standard  practice  square-root  of  the  sum  of  the  squares.  The  resulting  estimates 
are  summarized  in  table  A. 2. 
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Table  A.l  Estimates  of  Measurement  Uncertainty 


Calibration  Contam./Bkd.  Random  Total 


Value 

Abs. 

Rel. 

Abs. 

Rel. 

Abs. 

Rel. 

Abs. 

Rel. 

(U) 

1.0E-1 

1.0E-3 

1.0E-2 

0.0 

0.0 

6.2E-3 

6.2E-2 

6.3E-3 

6.3E-2 

(U) 

1.0 

1.0E-2 

1.0E-2 

0.0 

0.0 

1.1E-2 

1.1E-2 

1.5E-2 

1.5E-2 

(V) 

3.0E-2 

3.0E-4 

1.0E-2 

1.0E-2 

3.3E-1 

6.2E-3 

2.1E-1 

1.2E-2 

3.9E-1 

< V ) 

1.5E-1 

1.5E-3 

1.0E-2 

1.0E-2 

6.7E-2 

1.1E-2 

7.2E-2 

1.5E-2 

9.8E-2 

(W) 

3.0E-2 

1.2E-4 

4.0E-3 

1.0E-2 

3.3E-1 

6.2E-3 

2.1E-1 

1.2E-2 

3.9E-1 

(W) 

5.0E-2 

2.0E-4 

4.0E-3 

1.0E-2 

2.0E-1 

1.1E-2 

2.1E-1 

1.5E-2 

2.9E-1 

<u'u') 

1.5E-2 

3.0E-4 

2.0E-2 

5.0E-5 

3.3E-3 

3.7E-3 

2.5E-1 

3.7E-3 

2.5E-1 

<t/V> 

1.5E-2 

3.0E-4 

2.CE-2 

5.0E-5 

3.3E-3 

3.7E-3 

2.5E-1 

3.7E-3 

2.5E-1 

(l o'w') 

1.5E-2 

3.0E-4 

2.0E-2 

5.0E-5 

3.3E-3 

3.7E-3 

2.5E-1 

3.7E-3 

2.5E-1 

{u'v') 

2.5E-3 

5.0E-5 

2.0E-2 

5.0E-5 

2.0E-2 

9.1E-4 

3.6E-1 

y.lE-4 

3.6E-1 

( u'w ') 

5.0E-3 

1.0E-4 

2.0E-2 

5.0E-5 

1.0E-2 

1.8E-3 

3.6E-1 

1.8E-3 

3.6E-1 

{ v'w ') 

2.5E-3 

5.0E-5 

2.0E-2 

5.0E-5 

2.0E-2 

9.1E-4 

3.6E-1 

9.1E-4 

3.6E-1 

Table  A. 2  Estimates  of  Computational  Uncertainty 


Quantity 

Uncertainty 

«r 

3% 

Cf 

5% 

(wi) 

25%  -  50% 

(^z) 

10% 

y 

22%  -  45% 

(w2),yy 

40% 

(u'v'),yy 

100% 

(Uv),y 

20%  -  40% 
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Relative  Uncertainty  In  Measurements 


U-low  V-iow  W-low  <u’u'>  <wW>  <uV> 

U-high  V-high  W-high  <W>  <uV»  <vW> 


Figure  A.l  Relative  uncertainty  in  measurements 
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